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4 § Der Starkeffekt des Wasserstoffs bei kleinen Feldstärken 


| 


Inhaltsiibersicht 


ed Die Starkeffektaufspaltung bei Feldstärken unterhalb von 5000 Vem~! 
wurde an Umklappvorgängen in elektrischen Drehfeldern spektrographisch 
untersucht. Ein quadratisches Verschwinden der Feldaufspaltung in diesem 
Bereich, wie es Steubing und Junge direkten spektrographischen Messungen 
242 entnommen hatten, wurde nicht bestätigt. Für H, erreicht die Aufspaltung 
bereits bei 1500 Vem~! den von der Theorie geforderten Betrag. 


an den Linien H, und H,, des Wasserstoffs bei Feldstärken zwischen 0 und 
20000 Vem-! vorgelegt, aus denen die Autoren auf quadratisches Absinken 
der Aufspaltung bei Feldstärken < 4500 Vem und völliges Verschwinden 
unterhalb 2000 Vem-! schlossen. Nach Rechnungen von Lüders?), in denen 
Feinstruktur- und Lambshifteinfluß berücksichtigt wurden, ist ein derartiger 


«in Verlauf der erwarten. Bald nach Überschreiten einer 
gun Feldstirke Fy x 10° ‘x ist die 
bzüge Feinstrukturkonstante, a, der Radius der ersten Bohrschen Bahn Fu n 
we: die Hauptquantenzahl), bei der der Starkeffekt der Feinstruktur in die Fein- 


ıBer, struktur des Starkeffektes übergeht, sollte die Aufspaltung den normalen feld- 

stärkeproportionalen Verlauf der Schrédingertheorie zeigen (Fy liegt für 
ek: # #, bei 30 Vem-!). Gegen die von Steubing und Junge angewandte Methode 
sind zur Bestimmung der Starkeffektaufspaltung bei Feldstärken < 5000 Vem-! 

lassen sich Einwände erheben. Der Aufspaltungsfaktor wird von den Autoren 
Ver- aus der mikroskopisch gemessenen Differenz der Linienbreiten mit und ohne 
t ist. Feld erschlossen. Der aperturbedingte Dopplereffekt bei der benutzten An- 
rden. ordnung ist aber mindestens so groß wie die Feldaufspaltung bei 5000 Vem~. 


ittel- Die Differenz der Linienbreiten kann deshalb auch bei linearer Aufspaltung 
hren. den gemessenen Gang zeigen, wenn als das, was die Autoren unter „Linien- 
usge- breite‘‘ verstehen, ein genügend ungünstiges Maß vorgegeben ist (etwa die 


Halbwertsbreite der Linie). Es scheint nicht überflüssig, nochmais auf anderem 


) W. Steubing u. W. Junge, Ann. Physik 5, 108 (1950). 


Aun. Physik. 6. Folge, Bd.17 


Von Junge?) wurden Messungen Starkeffektes 


*) G. Lüders, Ann. Physik 8, an: (1951), siehe dort weitere Literaturangaben. 
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Wege eine Untersuchung des Starkeffektes bei kleinen Feldstärken vorzu- 
nehmen. Jede direkte spektrographische Bestimmung ist von vornherein wenig 
aussichtsreich. Einflüsse von Dopplereffekt, Intensitätsdissymmetrie und Ap- 
paratefunktion lassen kaum eindeutige Aussagen über die Aufspaltung zu. 

Für Mikrowellenuntersuchungen versagen die üblichen Nachweisverfahren, 
wegen der extrem kurzen Lebensdauer des metastabilen 2 S,,.-Terms im elek- 
trischen Feld (für E = 10 Vem! beträgt die Halbwertsbreite nur noch 10% 
sec®)); optischer Nachweis ist sehr unsicher, da es nicht gelingt, genügend 
unterschiedlichen Besetzungszustand benachbarter Starkniveaus zu realisieren. 

Es wurde deshalb im Anschluß an eine Arbeit von Ritschl®) eine Dreh- 
feldmethode benutzt, die es grundsätzlich gestattet, die Starkeffektaufspal- 
tung an H,, und H, zwischen einigen hundert Vem-! und 5000 Vem-! nach- 
zuprüfen. 


2. Untersuchungsmethode 
Schickt man ein H-atom (in angeregtem Zustand mit Hauptquantenzahl 


3 
n und Dipolmoment —ena,s in Feldrichtung; s = Starkquantenzahl: 


—(n—1)<s<m—1)) mit einer Geschwindigkeit v durch ein statisches 
elektrisches Feld, das längs einer Strecke d = v - At seine Richtung umkehrt, 
so behält das Atom nach dem Adiabatensatz seine Orientierung zum Feld 


bei, solange = wesentlich kleiner ist als die Termdifferenz &, benachbarter 


. . . + v 
Starkniveaus und es wird seine Orientierung zum Feld ändern, wenn T>% 


Wie weiter unten gezeigt wird, ist die Wahrscheinlichkeit, das Atom auch 


noch nach Durchlaufen des Drehfeldes mit einem Dipolmoment J eng: 


in Feldrichtung anzutreffen, fiir 72 ©, praktisch null und wächst für 


T< w, sehr schroff auf eins an. Die Messung dieser Wahrscheinlichkeit 


= 'z bestimmt den Betrag der Feldstärke, bei der die Aufspaltung 
Br 3 der Terme einen bestimmten Wert überschreitet. 

dy > Wir wählen folgende Versuchsanordnung: Ein H-Kanal- 

T strahl passiert ein starkes elektrisches Feld (zur Unter- 

2 suchung von H,, ~ 0,7—0,9-108 Vem-!). Dort wird ein Teil 

hen seiner Starkniveaus durch Präionisierung entvölkert, der 

Reet Kanalstrahl in einen unsymmetrischen Besetzungszustand 

isa überführt (von zwei Starktermen wird der energetisch 

#Ve -V. tiefere wesentlich stärker entvölkert). Die restlichen H- 

— atome durchlaufen danach das Umkehrfeld eines Kreuz 
kondensators (Abb. 1) längs der Symmetrieebene. 

Es wird die Starkeffektaufspaltung in dem homogenen 

“Vo *V. Feld vor und hinter der Kreuzungsstelle spektrographiseh 

aufgenommen und jeweils das Intensitätsverhältnis der rot- 

Abb.1. verschobenen Starkkomponenten zu den entsprechenden 


Feldanordnung violettverschobenen gemessen. 


3) G. Lüders, Z. Naturforschg. 5a, 608 (1950). 
4) R. Ritschl, Physik. 2. 40, -_ sees sowie Naturwiss. 38, 207 (1961). 
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3. Experimentelle Anordnung 


Zur Erzeugung eines möglichst lichtstarken, scharf gebündelten Kanal- 
strahls bei stabilen Entladungsbedingungen diente ein Ganzmetallentladungs- 
rohr, dessen Anodenzylinder in geeignetem Abstand von der Kathode (für 
P= 10! Torr, d ~ 2 mm) mit einer Ionenblende versehen war. (In ähnlicher 
Form von Hailer®) zur Erzeugung intensiver Protonenstrahlen entwickelt®).) 
Die Belastbarkeit des Entladungsrohres war begrenzt durch die Oberflächen- 
festigkeit der Kathode, die es gerade noch gestattete, das Präionisierungsfeld 
im Expositionsraum aufrechtzuerhalten. (Für 8-12 kV Brennspannung 
wurden Ströme von 8—14 mA erreicht.) Als Kathodenmaterial bewährte 
ich am besten hochlegierter, gehärteter Schnellarbeitsstrahl. Entladungs- 
raum und Feldraum waren getrennt durch einen Kanal von 0,25 x 2 mm 
Querschnitt und 30mm Länge. Hierdurch konnte im Feldraum mittels 
einer Pumpserie von Öldiffusionspumpen Leybold Mod.Q und Mod. P 
bequem ein Vakuum besser als 10? Torr aufrechterhalten werden. 

Die Platten des Kreuzfeldes waren durch Distanzstücke aus Quarz gegen- 


einander isoliert und durch Stellschrauben auf I mm genau justierbar. 


100 
Der Abstand betrug im Ionisierungsfeld 0,08—0,12 mm, im Kreuzfeld 1 mm. 
Feldspannung und Entladungsspannung wurden geglättet durch Konden- 
satoren von 1,5 uF bzw. 0,54F und auf 1% konstant gehalten, gemessen 
mittels Mikroamperemeter und Hochohm-Präzisionswiderstand. 

Die Untersuchungen wurden durchgeführt mit einem 1- und 2-Prismen- 
Spektrographen von 500 mm Brennweite und Öffnungsverhältnis f:5, der 
in der Institutswerkstatt von Herrn Mechanikermeister Herrmann gebaut 
wurde. Für Feldstärken oberhalb 75000 Vem-! konnte einwandfreie Fein- 
zerlegung der Starkkomponenten erreicht werden. Die Expositionszeiten 
lagen zwischen 4 und 6 Stunden. Für die Aufnahmen wurden Agfa-Raman- 
platten verwendet. Die Spektren wurden mit dem Schnellphotometer von 
Zeiß ausphotometriert. 


4. Berechnung der Übergangswahrscheinlichkeiten 

Es wird zunächst das Problem exakt gerechnet für den Fall, daß die Stark- 
effekteigenfunktionen unabhängig sind vom Betrag der Feldstärke, also für 
experimentelle Verhältnisse, bei denen die minimale Feldstärke des Dreh- 
feldes hinreichend groß ist gegen die Übergangsfeldstärke Fy (siehe 1), die 
Feldaufspaltung jedoch klein gegen den Abstand benachbarter Terme un- 
gleicher Hauptquantenzahl. Die Berechnung der Übergangswahrscheinlich- 
keiten reduziert sich dann auf die Lösung der Bewegungsgleichungen: 


(1) 
mit en und A= -(MxP-—-P M+- ‘) 
ı 2 mM x x 
‘ie’ C. Hailer, Wiss. Veröff. a. d. Siemenskonzern 17, 115 (1938), sowie Diss. Berlin 
*) Der Vorteil einer solchen Anordnung zur Erzeugung lichtstarker Kanalstrahlen 
wurde kürzlich auch von Gebauer [Z. Physik 185, 349 (1953)] quantitativ nachgewiesen. 


?) Vgl. W. Pauli, Z. Physik 36, 336 (1926) (unser A geht aus dem dort definierten 
durch Multiplikation mit n hervor). 
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(a, =: 2 ; n = Hauptquantenzahl, P = Bahnimpuls des Elektrons, M = 
0 

Bahndrehimpuls in Einheiten von hi). 

Die allgemeine Lösung von (1): y(t) = T - y(t,) wird vermittelt durch die 
sechsgliedrige Transformationsgruppe G, deren Tangentenraum durch die 
Operatoren A,, A,, A,, M,, M,, M, aufgespannt wird: 

T (G) = ei4z® gidy%y gids 9: gi eiMydy gi Mz 
Das gleiche Problem fiir ein magnetisches Drehfeld wird bereits durch die 


Transformationen der Drehgruppe allein gelöst®), wie man unmittelbar aus 
den bekannten Vertauschungsrelationen’): 


M,M,—M,M,=M,; M,A,—A,M,=A,; A, 4, 4; = My; 


(2) 
| M,4,— A; M,=0 
(j, k, l zyklisch) entnimmt®). Für ein ebenes Feld: 
E,=Fcod, E,=Fsind, E, = 0 (3) 


reduziert sich die zu (1) gehörige Transformationsgruppe auf eine dreigliedrige 
Untergruppe @ €G mit T’(@) = ei 42% Im Sonderfall 
F=0 und d = wz = const kann man dann (1) direkt (d.h. ohne explizit 
eine bestimmte Darstellung zugrunde zu legen) integrieren; nur dann 
lassen sich die Bewegungsgleichungen — durch die unitäre Transformation 
y=eMP mit tgO = — in eine Sehrödingergleichung 
mit zeitunabhängiger Hamiltonfunktion: 


F2 + h? wR A, (1') 


überführen, die man aus (1) erhält unter Benutzung der Relationen: 


A, ef Mi? — Ag cosd + (—1)Pixı A, sind 


= Ap cosB + (—1)Piazt M,sind 


cos + (—1)Psxx A, sin 
MO My =M cos + (—1)Pixz M,sin d 
_ 
i ABs für j,AK,L x, y,z durch Permutation 
iKL > 
((5) geht aus (2) hervor mit Hilfe der Beziehung U-!/(X) U = f(U1X U), 


die bekanntlich für jede, durch Potenzreihe darstellbare Ka, f (X) eines 
Operators X gilt.) Es wird dann 
Der allgemeine Fall (3): E, = F cos, E, =F sind, E, = 0 muß für 
eine spezielle Darstellung von @ gelöst, und 7’ (A, A ML, t) aus dem 


- 
8) E. Majorana, Nuovo Cim. IX, 43 (1932). J. Schwinger, Physic. Rev. öl, 
648 (1937). 


®) Vgl. S. Lie, Kontinuierliche Gruppen, Leipzig 1893, sowie Chevalley, Theory of 
Lie Groups, Princeton 1961. 
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‚ollständigen System von Lösungen der Anfangswertprobleme Y(t,) = Wo 
ntnommen werden. Experimentell interessieren nur Anfangsbedingungen, 
hei denen y, Eigenfunktion von EA ist, da sich vor Durchlaufen des Dreh- 
des jedes Atom in einem bestimmten Starkniveau befindet. 

Es sei also wo eee von A, (und von M,, zur eindeutigen Fest- 
kgung 0. B.d. A.): A, Yo = 8 Yq; M, Yo = = M Yo; dann ist die zu (1), (3) ge- 
hörige Lösung: 

w(t) = T'(A,, Ay, t) po = Av (6) 
irivialerweise Eigenfunktion von A; = T’ A, T’-! und M; = T’ M, T’ 
nit Ag p(t) = s p(t), Mey) = my(t). 

Aus (6) folgt wegen (5), daB A: =xM,+ BA,+y4, und M; = 
14, +6 M, +Y M,, skalare Produkte der beiden kommutativen ,,Vek- 
tor“operatoren A = an Ay A,} und M={4,, M,, M,} mit dem gleichen 
Einheitsvektor e; = {a, ß, y} sind. 

Da A, und M, und mithin auch A; und M: ein vollständiges System ver- 
tauschbarer Observabler im Unterraum der Hauptquantenzahl » bilden, 
folgt, daß y(t) eindeutig festgelegt ist durch den Einheitsvektor e;, die Inte- 
gation von (1), (3) zurückgeführt auf eine Berechnung der Richtung von 
el), wenn es(t) = {0,0, 1}, also der Transformation, die A, bzw. M, 
in de: bzw. Me; überführt. Mit Hilfe von (5) berechnet man: 
B 


T’ =e mit © = are cosy, arc tg (7) 

Wählen wir eine Darstellung, in der A, diagonal ist, so verursacht ® nur 
einen physikalisch uninteressanten Phasenfaktor und es genügt, den Winkel © 
m bereehnen. (Die Form (7) für 7’ ist für die Rechnung bequem, da die 
Matrixelemente von eiM=® durch einfache Transformation aus den bekannten 
Darstellungen der Drehgruppe erhalten werden.) 

Der niedrigste Grad einer Darstellung von @ ist vier (n = 2)?). Mit den 
bekannten Matrixelementen von A, und A,’) geht (1) über in das System 


simultaner Differe | a 
dc. a } 
dt 2 o™2 2 y~“3 177 
Cy = Ce = Cs = Cros = (8a) 
Die Transformation 


= mit U. 11 = = 1; U =U, an 


Ux, = 0 sonst 


» y 
M = 
1 die 
die 
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führt das System (8) über in a) 
db 

t y2 durch 
Wurz 
jurch 

und 
(8) ist dadurch fiir den Fall des zweidimensionalen Feldes auf ein System von } 
drei simultanen Differentialgleichungen zurückgeführt. Es sei: | 
1 1 
= 0,0 0 = — (A & = ; b, = — (ße — (9 
by = (0 By + bs =f, By; (Ba — (9) 


dann wird (8’) identisch erfüllt, wenn «,, 6, bzw. x, 8, Lösungen der beiden 
entkoppelten Differentialgleichungssysteme: 


= dx ; ; auf d 


S635 F INNE 


sind. Man überzeugt sich leicht, daß die Kombination aller unabhängigen 
Lösungen von (8”), in (9) eingesetzt, ein Fundamentalsystem von (8’) ergibt. 

Da wegen der eingangs erwähnten Anfangsbedingung Phasenfaktoren , 
uninteressant sind, genügt es, zur Bestimmung von © gemäß (7) die Gin. (8") jibe 
für die Anfangsbedingung «,(t)) = 0, 8, (tp) = 1 zu lösen. Aus (7) folgt dann 
mit Hilfe der Beziehung |«,|? + |ß,j? = 1 aus der zu n = 2 gehörenden Dar- 
stellung von #29; 


|B.(t)| = |eos (|. (10) 


Eine Darstellung von e‘M=® mit M, und A, diagonal kann aus der bekannten 
Darstellung der Drehgruppe durch unitäre Transformation erhalten werden 
an Hand der Entwicklung von parabolischen Eigenfunktionen des Wasser- 
stoffs nach den Eigenfunktionen bei Separation in Polarkoordinaten: 


(n, m, p) = S(n, m, p. I) (n, m, |) Abb. 2 
(s =n—2p—|m|— 1). 


Die Werte der S(n, m, p,l) sind von te elit angegeben worden: 


Ei 
stante 
S(n, m, p, l) (m* + p)! Der 


(n — p— 1)! (L— m*)! (1 + m*)!)}, wickle 
i 
Hierin 
für 


s 


10) V. Rojanski, Physic. Rev. 38, 1 (1929). 
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a) Umklappwahrscheinlichkeiten beim Durehlaufen des Kreuzkondensators 


Das Potential der im Experiment benutzten Feldanordnung (Abb. 1) wird 
dureh konforme Abbildung berechnet. 

Eine Parameterebene mit den aus Abb. 2a ersichtlichen Windungspunkten 
(Wurzel und Logarithmus sollen jeweils ihren Hauptwert annehmen) wird 


iach 

€ 


(s — 1) (# + 8— 4c] _ 

[im (le h, ih, (11a) 

| ‘1 auf das Innere des Kreuzkondensators (£-Ebene) und 

1/ 
durch 
(11b) 
= 

auf den (zweiblättrigen) Streifen der W-ebene zwischen U = + Vy und 

U=— V, abgebildet (vgl. Abb. 2b), wenn die Potentialdifferenz zwischen 

je zwei gegenüberliegenden Platten 2 V, beträgt. 

> U 

t» 7-0 


; 


Abb. 2. a) Konforme Abbildung des Kre vor ldes auf die Parameterebene ¢; b) Konforme 
Abbildung der Aquipotentiallinien auf die t-Ebene 


Ein H-Atom des Kanalstrahls durchläuft den Kreuzkondensator in kon- 
stantem Abstand Z, von der Y-achse (Abb. 2a) mit einer Geschwindigkeit v. 
Der Feldverlauf interessiert nur in der Umgebung von & = 0. Reihenent- 
wieklung von (11a) und (11 b) liefert 


V= 
Hierin ist 5 (€) < > wobei = 0, f 1, also 
2 i 2 ır 


für || S $A, zu vernachlässigen, wenn h, ~ hy. 
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2 t 2 (h? u 1(h? h2) ’ E, (1) 
(wenn das Atom zur Zeitt=0am Ort Y=0, Z=Z, ist). 
Mit (12) und den Abkürzungen 
r= zet K = iO (1?) 


(K ist das Verhältnis von Starkeffektaufspaltung im Minimum: «, ant 
B 


zur maximalen Drehfrequenz des Feldes: wp =F 


Die Substition «a =e 4 &; p= e* n führt (13) über in 
dr 
€ >)Ae und die Eliminierung von »: 
& 
d?& dé K 
wird gelést durch Integraldarstellung 
= =¢ fo (w) e®* dw 
H(t) =" wo (w) dw, 


wenn g(w) der Differentialgleichung iw + (w+k-+ = 0 genügt, 


und der Integrationsweg so gewählt wird. daß ’ wo(w).e”"—=0. Mit m= 


= 7 hat der Integrand in (19) 
bei wy = it einen Sattel; Richtung stark- 
ster Neigung ist w— wy) =s(1-+ 7%) mit 
reellem s. Zur Erfüllung der Randbedingung 
geniigt es, als Integrationsweg eine Gerade 
dieser Richtung durch w, zu wählen mit den 
Grenzen s = + co. Wegen der Singularität 
bei w = 0 muß für 7 > 0 der Integrations- 
weg die negativ reelle Achse ausschließen, 
Abb. 3. (Der Logarithmus nehme seinen 
Hauptwert an.) 

Abb. 3. Integrationswege; a) für Die Lösung interessiert für 
t>0; b) für r<o0 und + Co. des Inte: 
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granden (15) nach Potenzen von 2 und Integration ergibt unter Vernach- 
lässigung von Gliedern, die im lim t — + co verschwinden: 

2g(1 +i)Yr KT i= 

fir T> —oo: lim &(t) = 0; lim (t) = = (— 1) 


also wegen (10): g = — 


fir r> + oo: lim (t) = —— Yr(l+i)e (mi Ä 
@ (co) S 
(co)| —— | =e * =e (16) 
Aus (16) werden mit Für 27 ist B=1.5KV/em 
Hilfe von (*) die Matrix- RU re 
elemente von eiMz9(&) be- 
rechnet. Wir verzichten 
hier darauf, sämtliche be- | . = 


nutzten Matrixelemente | 


aufzuführen. Für das Ex- 
periment kommen haupt- 
sichlich die Elemente 


ei fiir die AuBen- 
-8,m,+ 8m 
komponenten in Frage, 
de die Zustandswahr- 


scheinlichkeit für diese Abb.4. Zustandswahrscheinlichkeit |a; (co)?| im Ab- 
Niveaus nach Durch- stand Zvonder Symmetrieebene, wenn |a; —(oo) |? = 1, 
laufen des Drehfeldes % (-») = 0 fiir k #j (für experimentelle Daten der 
angeben. Man erhält z. B. Aufn. 1) und 5) von Tabelle 2) 


n wg 


ws 2(n-1) 
(1 —e? =) . Den Verlauf als Funktion des Ab- 
standes Z des Atoms von der Symmetrieebene zeigt Abb. 4. (Die analogen 
Zustandswahrscheinlichkeiten für die übrigen Niveaus haben — bis auf gering- 


fügige Abweichungen bei = = 1 — ähnlichen Verlauf.) 
R 


b) Umklappwahrscheinlichkeiten bei feldstärkeabhängigem Dipolmoment 


Es bleibt abschließend zu diskutieren, welche Abweichungen von den 
zuvor erhaltenen Ergebnissen auftreten, wenn die Eigenfunktionen nicht 
mehr feldstärkeunabhängig sind. 

Die rein experimentelle Begrenzung in der Wahl der Ausgangsbedingungen 
bringt es mit sich, daß aus den Messungen zwar nur für Feldstärken, bei denen 
Ws/®z >1 ist, die Feldaufspaltung erschlossen werden kann, daß aber zu- 
gleich für 0 < wg/®g < 1 die Wahrscheinlichkeit dafür, daß sich ein 
Term nach Durchlaufen des Drehfeldes entvölkert hat, nicht wesentlich von 
eins abweichen darf. Wir geben im folgenden die Voraussetzungen an, unter 
denen dies auch dann der Fall ist, wenn infolge irgendeines Einflusses (H’) 
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das ne elektrische Feld noch nicht einen Starkeffekt mit ve „maxi- 
3 
malen“ Dipolmoment —ena, hervorbringt. (Die Feinstrukturaufspaltung 


stellt z. B. einen solchen Einfluß dar, ist aber für das Experiment ohne wesent- 
liche Bedeutung, da der Bereich Fin < 30 V/em nur etwa 1% des gesamten 
Kanalstrahllichtes ausmacht.) 

Mit einem Zusatzglied H’ in der Hamiltonfunktion lauten die Bewe. 
gungsgleichungen (H’ kann eine beliebige skalare, durch Potenzreihe dar. 
stellbare Funktion von r, M und o sein): 


a = [e F (z cos 8 + ysind) + H’]y. 


Hieraus ergibt sich durch die Transformation y = e‘Y** y, mit J, = M, +o, 
und Entwicklung nach Eigenfunktionen von e Fz + H’ = E: y = Say: 


(y, = y; (F)) für die Amplituden der einzelnen Zustände 


— =Ea+Y\ — ink (vi. 


f Bar 
mit a;=aje ° und (12) bzw. (12’) und “ae i folgt dann: 
da; ik T dy’, 
— ir Vice (vi; dF B)| 
(17) 
if 
.e 


oe. Fiir das Experiment interessieren im gr nur Lösungen a, (1) 


von (17) für > + 00, wenn a,(- oo) =1 und a,(- x) = 0 für alle 
K + j. 

Die Kenntnis der zugehörigen relativen Zustandswahrscheinlichkeiten 
lag (00)? mit X + j ist nur für im allgemeinen geringfügige Korrektions- 
glieder erforderlich. Aus der Struktur von (17) geht hervor, daß Lösungen 
a, (co) mit wachsender Hauptquantenzahl n (da jeweils n? Terme miteinander 
gekoppelt sind) immer unempfindlicher werden gegen den Betrag der Ma- 
trixelemente, sofern diese wenigstens in einem Teil des Zeitintervalles 


To 


-n<t<mjmit | - N Werte von der Größenordnung eins 
0 

nehmen, und daß Rn die a;(co) dann ausschließlich durch den Exponen- 
tialfaktor in (17) bestimmt sind, d. h. durch den Abstand benachbarter Niveaus 
©,x zwischen — T)< t< r, im Verhältnis zur maximalen Rotationsfrequenz 
des Feldes wp. Es sind demgemäß sechs Fälle zu unterscheiden: 

1. Die Eigenfunktionen sind feldunabhängig und J." im Unterraum der 
n irreduzibel, 
und 


oder 
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4 


3, Dipolmoment und damit Eigenfunktion sind für << F,") feldabhängig 
und _ für F (0) = B nicht notwendig irreduzibel, 
a) 1 und F (q) <F, 
b) t)< 1 und F (z,) 2 F, 
c) >1 und Fi )<F, 
d) >1 und F (1) = F,. 


Fall 1 a) und b) haben wir bereits exakt gelöst und festgestellt, daß 
‘a; (00) |? = 0 für 7,2 1 und |a, (co) |? > 1 für t)< 1 und zwar erfolgte der 
Übergang von null auf eins um so schroffer, je größer die Hauptquantenzahl 
war. Für 2 d) hat (17) genauso wie bei 1 a) den Charakter einer ,, Diffusions- 
gleiehung‘‘, da für |t|< r, keine eigentlichen Interferenzeffekte auftreten und 
für T > To: z —0. Wegen F(t,) x F, sind die H,, in (17) mindestens für 
einen Bereich At x 1 von der Größenordnung eins und eine numerische Inte- 
gration ergibt, daß dann für n = 5 der Term nach Durchlaufen des Kreuz- 
kondensators praktisch entvölkert ist (|a, (00) |? S 0,05). Die Fälle 2a), b), c) 
lassen ohne spezielle Annahmen über die H,x keine quantitative Aussage zu, 
als die, daß |a, (00) |? zwischen 0 und 1 liegt. 

Es lassen sich also Aufschlüsse über Abweichungen vom theoretisch zu 
erwartenden Gang der Starkeffektaufspaltung immer dann erhalten, wenn 
der Kanalstrahl den Kreuzkondensator in einem Bereich durchläuft, in dem 
theoretisch Fall 1a) in 1b) übergeht und für Abweichungen der Form 
® Kbeob. < ©; Ktheor. der Fall 2 d) noch über den Bereich 0 < Z< Z, (Abb. 3) 
hinaus realisiert ist und zugleich die Fälle 2 a)—c) für 0< Z< Z, durch geeignet 
gewählte Versuchsbedingungen ausgeschaltet sind. (Da nach Steubing und 
Junge F, S 5000 Vem-!, ist letztere Bedingung bei Kanalstrahlgeschwin- 
digkeiten von ~ 108 cms, einer Feldspannung von 10 kV und einem Platten- 
abstand von ~ 1 mm stets erfüllt.) 

Jede Abweichung dieser Art vermindert dann die Wahrscheinlichkeit 
'a, (00) |%, oder hat. wenn wir das mathematisch äquivalente Problem: Alle 
Terme gleichmäßig besetzt, nur der j-te entvölkert — betrachten, eine stärkere 
Wiederauffüllung des vorher entvölkerten Starkniveaus zur Folge. BER En 


Auswertung der Meßergebnisse 


Bei den erhaltenen Aufnahmen sind jeweils die rotverschobenen Außenkom- 
ponenten der H, (N = —20, —15, —10) durch Präionisierung geschwächt. Aus 
den Intensitätsverhältnissen der entsprechenden Komponenten werden die rela- 
tiven Besetzungswahrscheinlichkeiten Wf vor, und Wj nach Durchlaufen 
des Kreuzkondensators bestimmt. 

Der Kanalstrahl ist durch Spalte auf einen Bereich —Z,< Z < Z, begrenzt. 


Gemäß (7) wird dann : 
dZ|Z,. (19) 


Z, 
0 


U) Bei F = F, möge das Dipolmoment feldunabhängig werden. 
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Die Indizes sollen sich auf eine zum Feld orientierte Darstellung der y; und = 
y; beziehen, Z, ist Normierungsfaktor und definiert durch wg (Z,)/®p (Z,) =1 wi a 
(Abb. 4). Mit Hilfe der Unitarität von eM=9 kann (19) in die bequemere gi 
mittel 
Wi = WF + [75,4 2 le A 
gebracht werden. Hierin ist y gleich dem Verhältnis des Bereiches, in dem sicher 
O;K< Mp ist, zur Gesamtbreite des Kanalstrahls. Gesch 


18. 15, (10), bestimmt und in Abb. 5 als Funktion von Z,/Z, aufgetragen, 
Zum Vergleich ist der Verlauf von y (Z,/Z,) eingezeichnet den b). der sich 

bei einer Starkeffekt- 
aufspaltung entsprechend | —— 
den Meßergebnissen von 
Steubing und Jung | —— 
ergäbe. (Hierbei wurde 
zur Bestimmung der 
der Einfachheit halber 
mit feldunabhängigen AH,, 
in (17) gerechnet.) Ta- k 
belle 1 gibt die wichtig- 
sten /',, für n = 5 wieder 


Aus jeder Aufnahme wurde y unabhängig aus den Komponenten AN = der S 


hewir 


o) theoret Kurve 
¥ nach Junge 


5 (das Integral in (19) ist 
von 0 bis co genommen). 
Abb. 5. Integrale Umklappwahrscheinlichkeit y als Funk- Die eingeklammerten 


tion von Z,/Z, (Meßpunkte der Aufn. 1)—5) von Tabelle 2) Werte beziehen sich auf 
die zu Kurve b (Abb. 5) 

gehörigen Matrixelemente. Für k=s,m mit m +0 ist unter I’, stets --— 
m’)3 + zu verstehen, je nachdem ob s’=s oder § jy | 
s’ + s. (Siehe Tabelle 1). In Tabelle 2 sind Meßergebnisse und die zuge- 
hörigen experimentellen Daten 

aufgeführt. Spalte eins und Jig 0 


zwei geben die Intensitätsve-r 5 | 0 
hältnisse der rotverschobenen J” © 
zu den entsprechenden violett- 

verschobenen Komponenten für 
AN=+1, +15, +10 vor | 


und hinter dem Kreuzkonden- 
sator an, Spalte drei und vier § jg ( 
die hierfür hinter dem Kreuz- #15 = ¢ 


kondensator theoretisch zu er- 0 | ¢ 
wartenden Werte: Die nach 4) 
der Schrödingertheorie be- 8 | 0 


rechneten in der 3. Spalte, die 
Abb. 6. Starkeffektaufspaltunng der H, im Kreuz- nach den Ergebnissen von 
feld der Abb. 1. Kanalstrahlrichtung von oben Steubing und Junge (Sp. 4), 
nach unten. Die rotverschobenen Außenkompo- Spalte fünf h d sieben 
nenten sind vor der Feldumkehr deutlich ge- te 1b | ( 
schwächt den experimentellen und die Bw ( 
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theoretischen y-Werte für beide Fälle. In den übrigen Spalten die experimen- 
tellen Daten: Kanalstrahlbreite 2 Z,, Feldspannung 2 Vo, effektive Entladungs- 
spannung U (U beträgt, nach Geschwindigkeitsmessungen von Hailer l. c. 
mittels der Parabelmethode, 70% der angelegten Anodenspannung U’,: 
» = eonst- V Cet = const | "0,7 U,) und Minimumsfeldstärke B(Z,). 
alle Aufnahmen war h, = 0,6 mm, h, = 0,5 mm. (Siehe Tabelle 2.) 

Der experimentelle Fehier setzt sich zusammen aus den Einflüssen der Un- 
sicherheit in Z, (< 10%), der photometrischen Unsicherheit (< 10%) und der 
Geschwindigkeitsstreuung (2%). Stöße im Expositionsraum können wegen 
der Symmetrie der Anordnung nur eine Vergrößerung der Wiederauffüllung 
bewirken, also eine zu geringe Aufspaltung vortäuschen. 


Für 


Tabelle 1 


4;0 4;0 3; 1 3;1 
Ss m 
4;0 1,28 (1,21) 0,25 (0,4) 0,26 (0,17) 0,14 (0,14) 
-3; 1 0,13 (0,08) 0,07 (0,07) 1,34 (1,24) 0,25 (0,38) 
2;0 0,10 (0,07) 0,07 (0,06) 0,26 (0,17) 0,14 (0,12) 
2;2 0,04 (0,03) 0,03 (0,02) 0,19 (0,13) 0,12 (0,11) 
2;0 2;0 2;2 2;2 
m 
4;0 0,10 (0,07) 0,07 (0,06) 0,07 (0,06) 0,05 (0,05) 
—38; 1 0,13 (0,09) 0,07 (0,06) 0,19 (0,13) 0,12 (0,11) 
~2;0 1,50 (1,33) 0,19 (0,31) 0,15 (0,09) 0,09 (0,09) 
—2;2 0,08 (0,05) 0,05 (0,04) 1,39 (1,35) 0,23 (0,31) 
Tabelle 2 
WW w“ W* | dto. | dto. y ‚St Zo V, (Vert B (Z,) 
ya re  theor.) St. | beob. theor.| [kV] [kV] [kV em] 
we | wi 
18 | 0,61 | 0,70 | 0,70 0,80 | 0,255 0,26 | 0,42 0,49 10 60 1,30 
Ib 0,64 | 0,70 | 0,72 | 0,86 | 0,230 
10 | 0,9 
2 
18 | 0,24 | 0,48 | 0,42 | 0,60 | 0,30 | 0,25 | 0,41 0,49 5,0 7,0 : 1,65 
16 0,32 | 0,47 | 0,48 | 0,65 | 0,24 
10 | 0,43 | 0,68 | 0,66 0,88 | 0,26 
3) 
18 | 0,15 | 0,32 | 0,32 | 0,58 | 0,25 | 0,255) 0,42 049 #49 | 7,3 1,60 
15 0,20 | 0,39 | 0,37 | 0,53 | 0,26 
10 | 0,44 
4) 
Ik | 0,18 | 0,58 | 0,53 | 0,94 | 0,42 | 0,39 | 0,64 0,32 |49 | 7,3 1,60 
. 0,30 | 0,60 0,60 | 0,88 | 0,39 | SEEN 
9) | 
18 0,64 | 0,88 | 0,85 | 0,94 | 0,43 | 0,40 | 0,65 0,32 40 | 6,0 1,30 
15 | 0,68 | 0,80 0,84 | 0,93 | 0,36 
10 | 0,91 
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—_e dieser Fehlergrenzen bestätigen die in Tabelle 2 ui 
Ergebnisse die theoretisch zu erwartende Starkeffektaufspaltung für Feld- 
stärken zwischen 1,3 kV/em und 2 kV/em. Die mittlere Streuung der ge. 
messenen y-Werte beträgt +8%. Eine systematische Abweichung ist nicht 
festzustellen. 


Herrn Prof. Dr. Ritschl danke ich herzlich für die Anregung der Arbeit 
und das ständige Interesse, das er ihr gewidmet hat. 


Berlin, I. Physikalisches Institut der Humboldt-Universität. 


Bei der Redaktion eingegangen am 7. Juli 1955. 5 
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Beitrag zur Leitfähigkeitstheorie der Metalle 
bei tiefen Temperaturen 


Es werden Transporterse heinungen der Metallelektronen bei Anwesenheit 
von elektrischem Feld und Temperaturgradienten untersucht. Die ver- 
wendete Näherung für die gestörte Verteilung erlaubt eine anschauliche 
Deutung des Vorganges und liefert im Rahmen des einfachen Modells sinnvolle 
Grenzen für die linearen Gesetze. 


Von Koppe!) stammt eine Methode zur Berechnung der elektrischen 
Leitfähigkeit von Metallen, für die der Begriff der deformierten Fermi- 
Oberfläche charakteristisch ist. Das äußere elektrische Feld € treibt die 
Fermi-Kugel im f-Raum als Ganzes aus ihrer Ruhelage, wogegen die Wechsel- 
wirkung der Elektronen mit dem idealen Gitter ein diffusionsartiges Zurück- 
gleiten der auf der Oberfläche befindlichen Elektronen bewirkt. Die ge- 
störte stationäre Verteilung wird von Koppe durch eine „deformierte y 
Fermi-Kugel“ 

ge 


T, 5) =. 


+1 


I. Einleitung 


approximiert, welche sich von der Gleichgewichtsverteilung /, nur darin 
unterscheidet, daß an Stelle der Fermienergie ¢, eine von der Richtung im 
-Raum abhängige Funktion ¢ (9, y, T, ©) auftritt, die mit verschwindendem 
elektrischen Feld in ¢, übergeht. Bedingung der Stationarität ist: In einem 
zu einer beliebig vorgegebenen Richtung d, @ gehörenden infinitesimalen 
Kegel muß die Änderung der Teilchenzahl durch Diffusion (Stöße) gerade 
diejenige durch Drift (d. h. Einfluß des Feldes) kompensieren. Wegen der Ent- 
artung und auf Grund der gemachten Approximation resultiert hieraus eine 
Differentialgleichung für ¢, die bei Koppe an die Stelle der Blochschen 
Integralgleichung tritt. Sehr vorteilhaft an dem Ansatz für die gestörte Ver- 
teilungsfunktion ist, daß sie in jeder Richtung des f-Raumes ‚‚fermisch‘‘ von 
1 auf 0 abfällt, was ohne Schwierigkeit die Integration über |f| ermöglicht. 

Es erscheint nun wünschenswert, diese Methode so zu erweitern, daß ihr 
auch das allgemeine Transportproblem zugänglich wird, bei dem neben dem 


1) H. Koppe, Z. Naturforschg. 7a, 17 u. 156 (1952). ET 
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elektrischen Feld © noch ein Temperaturgradient grad T vorgegeben ist, 
Hierzu ist es erforderlich, näher auf die Struktur der Übergangszone ein- 
zugehen, innerhalb welcher die Verteilung von 1 auf 0 abfällt; die Wärme- 
leitung ist bekanntlich ein Effekt zweiter Ordnung in bezug auf die Tem- 
peratur. Will man obige Approximation verfeinern, ohne die damit verbundenen 
Vorteile aufzugeben, so besteht der einfachste Schritt in dieser Richtung 
zweifellos darin, neben £ auch die Dicke der Übergangszone richtungsabhängig 
zu machen, d.h. an Stelle von x T eine Funktion n einzuführen und 


f (f, T,&, grad T) = 


e”7 +1 


zu setzen. £ (0, p, T, &, grad T) und 7 (9, gp, T, E, grad T) sind — 
und müssen bei verschwindendem elektrischen Feld und Temperaturgradienten 
in ¢, und x T übergehen. Da nun zwei unbekannte Funktionen vorhanden 
sind, braucht man zwei Bestimmungsgleichungen. Sie ergeben sich aus der 
Boltzmanngleichung 


(D deutet auf Drift, S auf Stöße), indem man folgert: 


Kegel Kegel 


Hierbei wird über einen infinitesimalen Kegel in Richtung f integriert. Man 
stellt also im Kegel neben der Teilchen- auch die Energiebilanz auf; zweifels- 
ohne ist dies weniger als die Boltzmanngleichung allein, aber man braucht 
nicht mehr, da die grobe Form der Verteilung schon feststeht. 


Die folgende Darstellung stützt sich auf die einfachsten Annahmen, u. a.: 
Ideales Gitter, e = Schallquantengleichgewicht und Blochsehe 
Näherung für die Übergangswahrscheinlichkeit. Die hiermit verbundene 
Problematik soll nicht diskutiert werden, wohl aber wollen wir zum Schluß 
das geschilderte Verfahren dem üblichen gegenüberstellen. 


II. Anderung der Verteilung durch Stöße 


Wir beschäftigen uns zunächst mit der Frage, welche Änderung eine 
irgendwie aus dem Gleichgewicht gebrachte Verteilung f der Form (1) et 
fährt, wenn nur die Wechselwirkung mit dem Gitter berücksichtigt wird. 
Ausgangspunkt ist dabei 
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mit der Blochschen Übergangswahrscheinlichkeit für Metalle?): 


= 2 Na) + 1 —sgn [elt + m) — 
8 (le(t 10) — e(t)| — w). 


{und tw: Ausbreitungsvektor des Elektrons bzw. des en Schall- 
quants, «(f) Eagie des Elektrons, u, = Schallgeschwindigkeit, N(w) = 


4n 02 
—1 = 1jer Schallquanten im Gleichgewicht; ®, 
=const, hierbei ist C die Wechselwirkungskonstante in Bethescher 
Form, o die Dichte des Metalls. Die starke Entartung drückt sich Ze 
n<{ in f aus. Beschränken wir uns auf tiefe Temperaturen, also T<O 
ind nur „‚kleine‘‘ Schallquanten angeregt. Hieraus und aus dem in (3) ade 
den Pauliprinzip resultiert, daß Übergänge nur in unmittelbarer Umgebung 
der Fermi-Oberfläche e(f) —¢ =0 stattfinden können und dabei |tv|<|f| 
gilt. Daher kommt man mit der linearen Näherung e(f + w) —c(f) = 
w-gradre aus. ZweckmaBig ist die Einführung eines lokalen Koordinaten- 
systems f, dessen Ursprung der Endpunkt von f (Polarachse) ist, hier sind die 


Komponenten von w, = wsin cos w, = wsin dsin Q; Ww, = wcos ®. 
Wenn grads e wie im Falle freier Elektronen die Richtung von f hat, bedeutet 


I<d< Absorption, =< 9< 2 Emission eines Schallquants, und es 


gilt 
e(t + w) — e(f) = w- gradre = W, (k) = e’(k) weos®, €’ (k) 
Also wird die Übergangswahrscheinlichkeit: 
® 1 = 
= ie N(w) + 1—sgn [5 — #]}5 (cos (©) 
mit y = Da wir uns nur für f-Werte in der Umgebung der Fermi- 


Oberfläche zu interessieren brauchen, stellt y das Verhältnis von Schallge- 

schwindigkeit zur Geschwindigkeit an der Oberfläche dar; bei Metallen ist | 

y<1 und wird fortan als konstant angesehen. ‘tow Ich 
Das erste Integral in (3) wird nun unter Beachtung von (5): 


2 n)? (2 


+ (N(w) + 1) f dd sin #6 (cos + y) {1 + w)} 


ry 


*) Abgeleitet etwa bei Sommerfeld-Bethe®) S. 514, oder Wilson‘) S. 258. 
3) A. Sommerfeld u. H. Bethe, Hdb. XXIV. 2, 1933. 
*) A. H. Wilson, Theory of Metals, Cambridge 1953. 
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u | 
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Energie- und Impulssatz liefern die Bedingung cos ? = + y; die erlaubten jy 
liegen auf einem Kegel mit Offnungskosinus y, 
wir bezeichnen sie kurz mit tv, bzw. tw_.,, (vgl. 
Abb. 1). Damit wird aus (6): 


Wy 


dw w | N (w) f dp {1 — f (£ + w,)} 


5 


+ (N(w) + 1) fi dp {1 — f (f+ m3 (7a) 
0 


und ganz analog bekommt man fiir das zweite 
er aus (3): 


Abb. 1. Energie-Impulskegel 


f dw w? |(N (w) + 1 dp f (£ + w,) + N dg f 
(7b) 
Die mit t,, verknüpfte Impulsänderung führt von der Richtung #.» 
zur Richtung 9’, g’, d.h. 
Da |w.,|<|f|, können wir, sofern unphysikalische starke Schwankungen 
von ¢ und n in der Umgebung der Richtung #, g ausgeschlossen sind, nach 
g=LW,p) und entwickeln. Mit den 


Abkürzungen 
£ (8, p) Uy iw xT 
gibt das: 
1 
(8) 


Die p-Mittelung über den Kegel wirkt sich jetzt nur noch auf g und h aus, 
und es ist leicht zu bestätigen (App. I), daß gilt: 


(1-79) 22 


3 A; hdp = An. (9) 
2 2 


A = Laplaceoperator auf der Einheitskugel. 

Durch eine einfache Zwischenrechnung, bei der man (8) in (7) einsetzt, 
(9) verwendet und noch die Beziehung N(w) + 1 = — N(— w) benutzt 
findet man für (3): 


-@. = (35). I: 


®, (1 1 
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Fir (=, » =xT liegt die ungestörte Fermiverteilung vor. Wegen 
At = An = 0 fällt dann der zweite Term weg. Aber auch das erste Integral 
verschwindet, da die geschweifte Klammer für A = 1 Null wird, was man leicht 
einsieht. Daher ist (öf/öt)s = 0, wie es sein muß. 

Im folgenden betrachten wir nur Verteilungen, bei denen 7 nicht zu sehr 


von x T abweicht. 2 = liegt dann bei 1, und wir können die geschweifte 


Klammer im ersten Integral von (10) an der Stelle A = 1 entwickeln: 

Im zweiten Integral darf in dieser Näherung 72 = 1 gesetzt werden, was auf 

Vernachlässigung von Gliedern proportional (A — 1) g und (A — 1) A hinaus- 


läuft. Aus (10) wird dann: 
oT 


Pr (11) 

2 (27) (5) - | y| + y) (N (y) +1 (z)) 
ig [Ad + (2 + y) An]. 

Ohne Schw ierigkeit kann man jetzt die Teilchenzahl- und Energieänderung 

» einem Kegel in Richtung d, m des f-Raumes erhalten; dazu ist lediglich 


00 
(4), k? dk und | ek? dk auszurechnen. Wegen =‘ “ it dz = 
) 0 


“ak, und ein Blick auf (11) zeigt, daß wir Sak, ©, dk und £ dk 


brauchen. Hierbei ist die Kenntnis von e in dee von 


2 72 
lich. Wir nehmen ¢ = und finden mit 1: hat 
2m ay 
de =F Sdk = 


Der ee an pei Integrals in (11) hat die Eigenschaft G, (y, z) = 
~G,(—y, — 2). Beim zweiten Integral gilt dies nur für den An proportionalen 
Teil, daneben gibt es noch ein @, mit Gs bu 2) =@, (—y, — 2). Wegen ne <1 


entspricht “+ k< co praktisch -- z< co und wir haben 
fa dy G, (y, 2) = fa J =0 
zu berücksichtigen, "Mon bekommt mit (11) und (12): 
0 
m 
LY" 7) An | (18) 
22D) (xT\3 „m? (9 
(2 7)? G5) Fa Mo (F) 


(3) > K o( Ab. 
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Hierbei ist: 


co 


= de J dy |y| yf (2 + y) {N(y) + 


00 


= fa dylyl + wef Ny 


co 
= fd fayly (2+ y) {N (y) + 1—f 
-00 
Nach App. II besteht mit den bekannten Integralen J, (x) der Zusammen- 
hang: 


22? 
Ky (a) =— 255 (a); Lg (a) = (0) Tl); My (a) =J5 (0). (14 


Da A=xT/n voraussetzungsgemäß wenig von 1 abweicht, können wir 


sehreiben: 


*T\2 (ugh)? th 


rdw = -f dk 
und Energie 
| édw = de 
im Kegel der Öffnung dw: 
-(F)e [K AC L (5) An 
73 mus 
-(5), dw = dw [x(5 + (F (y — x 1)]. | 
Die K, L, M ergeben sich aus den Ky, Ly, M, durch Multiplikation mit — 
Dy m (1 — y*) (x T\5 
Bevor wir (15) diskutieren, wollen wir ohne Rec aia noch auf eine andere 
Méglichkeit hinweisen, um zu diesen Gleichungen zu gelangen: 


Kin Elektron sei bei f. Dann ist seine mittlere Verschiebung durch Stoß 
pro Zeiteinheit: 


u, (f) = fw {1 — (f + w)} Wee, 
Befindet sich bei f kein Elektron, ist also dort ein ,,Loch“‘, so ist dessen mittlere 
Verschiebung 
Uy (f) = f wwf rm) 
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Man kann auch sagen: — u, bedeutet die mittlere Verschiebung eines Elek- 
trons zu f hin, wenn dort ein „Loch‘ ist. Da f(f) bzw. 1—f(f) die Wahr- 
scheinlichkeit darstellt, daß sich bei f ein Elektron bzw. ein ‚‚Loch‘‘ befindet, 
erscheint es sinnvoll, 
(a) Sf) = FF (Hy — (O} 
als mittlere Strömungsdichte zu betrachten. 
jy = fy u, heißt dann, daß im Mittel ein 
von f nach F durch ein von f’ nach f gehendes 
Elektron ersetzt wird, so daß sich nichts ändert. 
Die Bestimmung des Vektors © hat natür- 
lich große Ähnlichkeit mit unserem bisherigen Vorgehen, und die Berechnung 
(14) des Teilchen- und Energieflusses in einem Kegel auf Grund der Gleichungen 


| Wir on > we > i 
do dive S dt = En do 
= div S = do — - graded 


Abb. 2. Mittlere Strömungs- 


men- dichte 


(integriert wird über das Volumen, bzw. die Oberfläche eines Kegels der 
Öffnung dw) gibt in der Tat ebenfalls die Formeln (15)5). 


| 
Jetzt können wir (15) anschaulich deuten. Die Teilchenzahl im Kegel ver- 
mindert sich hauptsächlich dadurch, daß die beweglichen Elektronen in der 
Ubergangszone in Richtung von — grad; ¢ (K (O/T) ist < 0) ausdem Kegel zu 
kleinerer Energie hin herausströmen. Eine mit dem Faktor »/¢ behaftete 
Korrektur, die grad: 7 also dem ,,Temperatur‘‘-gradienten auf der Fermi- 
fläche proportional ist, tritt hinzu. Natürlich ist mit den herausströmenden : 
Elektronen eine Energieverminderung im betrachteten Kegel verbunden, was 
sich offensichtlich im ersten Term von (15b) bemerkbar macht. Aber die ' 
Strömungsdichte © hat auch eine Komponente in Richtung f bzw. grad; e, 
P die auf die Teilchenänderung im Kegel keinen Einfluß hat, wohl aber auf die 
15a) Energieänderung, man erkennt das schon an (16). Diese Komponente gibt 
Anlaß zum zweiten Term in (15b). Die Energieverminderung ist proportional 
15b) n—xT'; die Wechselwirkung mit dem Gitter veranlaßt die Elektronen im 
Kegel, mit ihrer formalen ‚‚Temperatur‘‘ 4/x nach der Gittertemperatur hin- 
zustreben. 


Betrachten wir z. B. ein Elektronengas mit der Temperatur 7’ im Gleich- 
gewicht (7 = x T’), und denken wir uns die Wechselwirkung mit einem Gitter 
der Temperatur 7’< T’ eingeschaltet, so kühlt sich das Elektronengas auf 

dere die Gittertemperatur ab, in dem es gemäß (15b) sekundlich die Energie 


Stoß 32 mu, (7) T-T 
1—-yp 
5) Man muB sich aber vor Augen halten, daß die Strömungs- bzw. Diffusionsbeschrei- 


tiere bung des Stoßprozesses |10|< |f| zur Voraussetzung hat. Daher ist eine Gleichung 


ö 
(2), + div, S = 0 nur sinnvoll, wenn sie über ein Volumen V, (z. B. unsere Kegel) inte- 


griert wird, dessen Ausdehnung im f-Raum groß gegen die wahrscheinlichsten Stoß- 
änderungen iv ist. 
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abgibt. Da wir immer Schallquantengleichgewicht vorausgesetzt haben, 
wirkt das Gitter hierbei als Wärmereservoir. 

Allgemein gibt (15) quantitativ darüber Auskunft, wie eine gestörte Ver- 
teilung im Rahmen unseres einfachen Modells durch Gitterwechselwirkung in 
den Gleichgewichtszustand /, zuriickgleitet. 

Pewee 
III. Stationäre Verteilung 


Um die Änderung von % und é durch Drift zu bekommen, nn wir von 


ad T + - grad f 


aus; © und grad T mögen konstante, parallele Vektoren sein, die in der Folge 
die Polarachse z festlegen. Man hat: 


oo co oo 
| (2). kK? dk = (grad T). cos | k3 dk | grad, f) dk 
0 0 0 
(1%) 
=), ek?dk = —— (grad T), cos + | © (€ grade f) dk. 
0 0 


Da f sich in f-Richtung gemäß (1) fermisch verhält, können die in (17) vor- 
kommenden Integrale in bekannter Weise ausgerechnet werden, wobei man 


sich auf — | G(x) jf (x) dx =G(0) + =G ’’ (0) stützt. Das gibt (App. ID): 
-00 
en 
F ö a 


+ LAPS +7 | 2008 @. 
— 2m x (grad 7’), 
’ (22h) 


eingeführt. Die Driftänderung verschwindet mit Fy und F,. — E und 
grad T nicht in einer Richtung, oder ist auch noch ein Magnetfeld vorhanden, 
so werden die Formeln (18) entsprechend komplizierter und enthalten Ablei- 
en nach 9. Die Verteilung f ist stationär, wenn gilt: 


+ (2), = 0; (5), 4 
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Mit (15) und (18) kommen wir daher auf das Gleichungssystem: ee 
ar 
KAL- = OC) + # pe 24 — BS n)2 
(20a) 


rle+? „2 on 9 

Es ist unter der Bedingung [ n dw = n = const und: 
>0, n>»T 


zu lösen. 


IV. Näherungslösung für Z und 7 
Wenn £ und » nur wenig von den ungestörten Werten abweichen, also bei 
« 1, 1, (21) 
kann man leicht eine gute erste Näherung angeben. Man setzt zu diesem 
Zweck in (20) rechts fy und x T an Stelle von £ und », ferner in (20b) links noch 
=x {und werden nicht von abhängen, sodaß A = 5 sind 
Damit wird aus en 


sind dd sin fy sin’ dd (A + Py, 6 70 2 cos d 
(22a) 
K ind a | —x T) 
-[r (a+ Sox (1-F(% 208 0 
3 9 
at, 
Hierbei wurde = -5 verwendet. Aus (22a) folgt: 
So 


oder: 


= 


Da 5 und » keine Singularitäten haben dürfen, ist const = 0, also: 
F , F,). (23) 


SchlieBlich wird (22a) 
gleichung fir n—xT: 
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die Differentialgleichung der Legendreschen Polynome mit Inhamogeailll 
Für T — 0 geht A gegen co; etwa im Bereich T/O < 0,1 ist A > 2. Dort hat 
(24) die Lésung FT 


1 


— cos 3. 
Setzen wir sie in (23) ein, so resultiert, da c (F, F,) durch ¢ = {, fir F, = F, =0 
festgelegt ist: 


24% T 
3 


Diese Näherungslösungen können als Entwicklung nach Kugelfunktionen be- 
trachtet werden, die nach P, abbrechen; sie sind in Feld und Temperatur- 
gradient linear. 

Allerdings sieht man, daß ¢ und »; sehr unterschiedlich von F, und F, ab- 
hängen: In B hat F, gegenüber F, noch den Faktor x T/C, bei C ist es gerade 
umgekehrt. Das legt es nahe, überhaupt die mit x T/£, behafteten 
Terme gleich Null zu setzen. Der Fehler ist ~ x T/£,- 


& ant 


K) 
In dieser Näherung ist daher das elektrische Feld allein pe die Deformation 
¢(#) der Fermikugel verantwortlich, während der Temperaturgradient ledig- 
lich die Dicke der Ubergangsschicht n(d) 
bestimmt. Also ergibt sich von selbst, daß 
zur Behandlung von elektrischen Proble- 
men nur die Teilchenbilanzgleichung 
nötig ist. Ferner ist die Verknüpfung des 
n mit dem Temperaturgradienten in sehr 
anschaulicher Weise plausibel: Betrachten 
wir am Ort r mit der Temperatur T die 
Geschwindigkeitsverteilung. Dann werden 
Abb. 3. Fermi-Kugel am Ort r mit ‘ie Elektronen mit einer Geschwindigkeit 
Temperaturgradient in Richtung grad T aus Bereichen T’<T 
+ kommen, während bei Geschwindigkeiten 
a entgegengesetzt grad T gerade 7” >T ist: Die Fermikugel kehrt der 

one „Wärme ihre kalte Seite, der „Kälte ihre warme Seite‘ zu (Abb. 3). 
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Die k-Integration läßt sich wieder ausführen (App. III) und liefert: rae 


2 
Wenn wir = + und = x T(1 die Voraussetzung 
0 


(21) berücksichtigen, erhalten wir aus (28) in RE Nabsrungs 


1 
4mel, ¢ | | 


1 
= ja J de x 


i, =— 


Man erkennt, daß in einer Entwicklung von € und n nach Kugelfunktionen 
lediglich das Glied mit P, = x einen Beitrag gibt. Ist nur ein elektrisches 
Feld da, so folgt mit (27) und (19) aus (29) das Ohmsche Gesetz: 


Hat man dagegen einen Temperaturgradienten vorgegeben, so ist i = 0 mit 


gleichbedeutend. Dies ist die Bestimmungsgleichung für das Thomson- 
feld. Mit ihr wird der Wärmestrom 


3 (2 2)? 


) wurde dabei gleich 


Ein gemäß (29) noch vorhandener Faktor 1 -5(£ 


1 gesetzt. Nach (27) und (19) kommt: Be 
m x? TE 
— gy grad (32) 


Einsetzen der Ausdrücke für K, L und A liefert nach (30) die elektrische Leit- 


fähigkeit 
T\5 3 ht 
-=24(7)(5) 


und nach (32) die Wärmeleitfähigkeit 
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während aus (31) das Thomsonfeld 
Em ae grad 7 


resultiert. Für 7 — 0 ergibt sich das T5-Gesetz der elektrischen und das 
T-*-Gesetz der Wärmeleitfähigkeit, während das Thomsonfeld proportional 
T® grad T ist. Bei der reinen Wärmeleitung geschieht folgendes: Ein vorge. 
gebener Temperaturgradient beeinflußt nur die Übergangsschicht » der 
Fermikugel (Abb. 3). Dies hat neben den Wärmestrom auch einen elektrischen 
Strom zur Folge wie (29) zeigt. Das sich ausbildende Thomsonfeld führt 
aber zu einer zusätzlichen Deformation ¢, welche einen elektrischen Strom 
hervorruft, der jenen gerade kompensiert. 


> ise I. Vergleich mit dem üblichen Verfahren 


s in von Bloch und Wilson genau überein, wenn 1 — y? = 1 gesetzt wird. Dies 
kann man chne praktischen Fehler tun, denn z. B. bei £, = 4 eV ist y von der 
Größenordnung 10-4. 
Üblicherweise wird in der Leitfähigkeitstheorie für die gestörte Verteilungs- 
funktion der Ausdruck ae 
f = — cos f, (e) 
af 
he 
gilt. Damit folgt | j, (&) !<2% T und man findet als Gültigkeitsbereich des 
Ohmschen Gesetzes: 


verwendet. Allerdings ist das nur dann eine gute Näherung, wenn f,> 


Ihkol _- xT 
Ene (34) 


so 


Koppe!) hat darauf hingewiesen, daß diese temperaturabhängige Schranke 
lediglich eine Eigentümlichkeit der Lorentz-Sommerfeldschen Näherung 
ist. In der Tat resultiert bei uns aus (21) und (27) die Linearitätsbedingung: 


-K 
| 1 Me 


Es gilt daher das Ohmsche Gesetz, solange sich der Schwerpunkt mi /ne 
der gestörten Verteilung wenig gegen den Radius m v, der ungestörten Fermi- 
kugel im Impulsraum verschiebt, ein sehr plausibles Resultat ®). 


6) Die Beziehung (34) begrenzt von unserem Standpunkt aus nicht das lineare Ohm- 
sche Gesetz, sondern gibt an, bis zu welchen Stromdichten die gestörte Verteilung noch 
im Bereiche der Übergangszone x T der ungestörten Fermiverteilung f, von 1 auf 0 
\ abfällt. nev,x T/£, ist allerdings recht groß und hat fiir 6, ~ 5eV, 7’ = 1° K den Wert 
A x 10’ A/em?. Experimentell gefundene Abweichungen vom Ohmschen Gesetz bei 
je ~ 10% A/em?)®) sind wegen (35) sicherlich nicht im Rahmen der vereinfachten Theorie 
zu erklären. 

?) A. A. Ignatjeva u. C. C. Kalaschnikow, J. exp. teor. Fiz. 22, 385 (1952). 
8) E. S. Borowik, Doklady Akad. Nauk 91, 771 (1953). 


Die oben angegebenen Formeln fiir o und K stimmen mit den entsprechenden 


Wühlschl 


Auc 


nd (2 


Die mi 


daher 
portior 


wo le 


Wi 
verallg 
zu der 
die Ab 
sichtig 
die ge 
welche 
wonne 
Exper 
Richt 
Reche 
Bard 
ist 
Umkl 
nieht 
und 
könn: 


Ic 
und 


0 Annalen der Physik. 6. Folge. Band 17. 1956 we mz 
- ate 
| 
a 
9) 
Lond 
1 
G.K 
194 | 
i 
1 
1 


1 das 
‘ional 
orge- 
der 
schen 
führt 


trom 


nden 
Dies 


n der 
Ing3- 


(33) 


de 


und F;| < (A—2) und in (32) eingesetzt: 
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Auch der Fouriersche Ansatz kann gerechtfertigt werden, denn mit (21) 
3 «xT 


2m Cox T (x T)? 


n? (x T)?, 
C 


Die mittlere Energie eines Elektrons im Ferm igas ist ¢(7’) (0) #% 


daher ergibt sich als Bedingung für einen dem RER In. pro- 
portionalen Wärmestrom: 


9 
c< nAé v%, 
wo Je =e(T) —e(0) der Überschuß über die Nullpunktsenergie 260) ist. 


Wir diskutieren zum Abschluß kurz, inwieweit die hier geschilderte Methode 
verallgemeinerungsfähig ist und wo ihre Grenzen liegen. Zunächst ist daran 
zu denken, an der recht guten Näherung (1) eine Korrektur y (f) anzubringen, 
die Abweichungen vom fermischen Verhalten der gestörten Verteilung berück- 
siehtigt. Dies wurde bereits von Koppe angedeutet. Jedoch können auch 
die genaueren Lösungen der Boltzmannschen bzw. Blochschen Gleichung, 
welche mit Hilfe der Variationsmethode durch sukzessive Näherungen g>- 
wonnen werden’), nicht die schlechte quantitative Übereinstimmung mit den 
Experimenten!) bei der Wärmeleitung überwinden. Ein Versuch in dieser 
Riehtung, der sich an die Diffusionsnäherung anlehnte!!), beruhte auf einem 
Rechenfehler!2). Es erscheint erfolgversprechender, nach Ziman!?) vom 
Bardeenschen Ansatz für die Gitterwechselwirkung auszugehen. Allerdings 
ist unsere Methode auf kleine Impulsänderungen beschränkt, weshalb ihr 
Umklapp-Prozesse oder Transportvorgänge im nichtentarteten Gas zunächst 
nieht zugänglich sind. Dagegen ist das Verfahren auf starke Felder zugeschnitten 
und sollte darüber hinaus auch eine Schallquantenverteilung erfassen 
können, die nicht im Gleichgewicht ist. 


Ich danke Herrn Dr. Koppe sehr herzlich für Ratschläge und Diskussionen 
und bin Herrn Prof. Möglich für freundliche Unterstützung verbunden. 


Appendix I 
Im f-System . bei £ = 0, d. h. am Endpunkt des Vektors f: es 
at et 


®) M. Kohler, Z. Physik 125, 679 (1949). E. H. Sondheimer, Proc. To Soc. 
London (A) 208, 75 Dada. 

%) R. Berman u. D.K.C. MacDonald, Proc. Roy. Soc. London (A) 209, 368 (1951); 
2 White, Proc. Phys. Soc. 66, 844 (1953); P. G. Klemens, Proc. Phys. Soc. 67, 
194 (1954). 

1) M. Toda, J. phys. Soc. Japan 8, 339 (1953). 
18) A, Blatt, J. phys. Soc. Japan 9, 440 (1954). 
8) J. M. Ziman, Proc. Roy. Soc. London 226, 436 (1954). 
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= ). Daher lautet die Taylorentwicklung: und mit 


00 
f ee 
0 
Wegen w, = w/1—y? cosq, w, = w yi—y sing ergibt sich hieraus 
f gdy = ya (= ~ Transformation von f auf f führt zu (9), 
Man bringt die Integrale zunächst auf die Form 1 [a 
de 
Ko(a)= — Wo (y), 
2 Ber 
= — J dy| y| y® N(y) er (we(y) + y Wiy) 
Myx) J dy|y| (ye (Wy(y) + Waly) 


-& 


co 


mit We. dz z, 


Aus der Reihenentwicklung 4) ta 


je 2 [8/2] N 
= (— 12-27 (° (2 r)! cay yorı-ar 


Cy = 4, Cg, = (1— 21-27) (27), f(a) = Zetafunktion, folgt: 


1 
——, W, 


Damit kommt man auf den in (14) festgestellten Zusammenhang. Es ist 


J, (a) = dy y" N2(y) ev = 


(17) läßt sich mit E- grad, f = dive e qradyf div; 
— €-gradze vereinfachen. Wegen div; A = (k? — 5: 2 (sin?) 
(Wist ein Vektor, der nur eine Komponente 1 in Richtung der Blase :hse hat) 
folgt mit A iG: 


E grad; k2 dk = sin? fde 


2 sind eb 


Vgl. Wilson‘) S. 336. 


t 
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= 
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f E- grade = —™% - (sine Fe fa + 20080 f ef 

|sin® 

hi N of 3 of feast on 

!eraus Ferner wird k3 dk = -2(5) + ; =) 

u (9), 

und Sh = 7 oF 
Benutzt man 


2 
de = +- n(n — 1) eats de =a or 


wergibt sich (18) und (28). 


Berlin, Institut für theoretische Physik der Humboldt-Un Universität. 


‘Bei der Redaktion eingegangen am 15. Juni 1955. 


% 
hat) 


4 


“ 


Rotverschiebung und Bewegungsgleichungen 


Von A. Papapetrou 


Es wird eine unmittelbare Berechnung der relativistischen Rotverschiebung 
mit Hilfe der Feldgleichungen des kombinierten Gravitations- und elektro. 
magnetischen Feldes sowie der Verallgemeinerung der Dirac-Gleichung für 
einen Riemannschen Raum gegeben. 


Die WIEN Relativitätstheorie hat bekanntlich drei neue Effekte vor- 
ausgesagt: Die Periheldrehung, die Lichtablenkung und die Rotverschiebung. 
Die Formeln für die ersten zwei Effekte leitete Einstein aus dem urspriing- 
lich unabhängig von den Feldgleichungen postulierten geodätischen Bewe- 
gungsgesetz ab. Für die Berechnung des dritten Effektes war eine weiter 
Hypothese notwendig. Nach dieser sollen Atome, die sich an verschiedenen 
Orten innerhalb eines Gravitationsfeldes befinden, sich wie Normaluhren ver- 
halten, d. h. Spektrallinien emittieren, deren Eigenperioden — am Emissions- 
ort gemessen — vom Gravitationsfeld unabhängig sind. 

Später konnte man beweisen, daß das geodätische Bewegungsgesetz kein 
unabhängiges Axiom bildet. Wie es zuerst Einstein und Grommer zeigten)), 
sind in der allgemeinen Relativitätstheorie die Bewegungsgleichungen durch 
die Feldgleichungen festgelegt. Die Bewegungsgleichungen darf man also 
nicht willkürlich postulieren, sondern muß sie aus den Feldgleichungen ab- 
leiten. Die ausführliche Behandlung der Bewegung eines Probeteilchens hat 
dann tatsächlich ergeben, daß für diesen Fall aus den Feldgleichungen das 
geodätische Bewegungsgesetz folgt?). Dadurch wurden aber offenbar die 
Periheldrehung und die Lichtablenkung, die sich beide auf Bewegungen von 
Probeteilchen beziehen, unmittelbar auf die Feldgleichungen der allgemeinen 
Relativitätstheorie zurückgeführt. 

Eine ähnliche Zurückführung der Rotverschiebung auf die Feldgleichungen 
wurde bisher nicht gegeben. Dies stellt zweifellos keine praktisch wichtige 
Lücke dar. Das Verhalten der Normaluhr im Gravitationsfeld folgt nämlich 
unmittelbar aus dem Äquivalenzprinzip. Allerdings liegt dabei das Vorbild 
einer makroskopischen Uhr zugrunde, welches bei der Diskussion der Rot- 
verschiebung nur im übertragenen Sinne herangezogen werden kann. Um den 
Unterschied völlig klar hervortreten zu lassen, wollen wir an das Haupt- 


1) A. Einstein u. J. Grommer, S.-B. Berl. Akad. 2 (1927). 
2) L. Infeld u. A. Schild, Rev. mod. Physics 21, 408 (1949); A. Papapetrou, 
Proc. Roy. Soc. London (A) 209, 248 (1951). 
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merkmal einer (makroskopischen) Uhr denken: Es handelt sich im wesent- 
lichen um zwei materielle Punkte, die infolge einer Koppelung sich gegen- 
einander so bewegen, daß eine Reihe von Koinzidenzen resultiert. Bei Ab- 
wesenheit eines Gravitationsfeldes wiederholen sich diese Koinzidenzen nach 
gleichen Zeitintervallen, die der Periode der Uhr gleich sind. Ist nun diese 
Uhr in ein Gravitationsfeld gebracht und dort freigelassen, so wird man dann 
ihre inneren Bewegungen am zweckmäßigsten in dem mitbewegten Koor- 
dinatensystem betrachten. Da dies aber ein frei fallendes Koordinatensystem 
ist, wird in ihm das Gravitationsfeld bis auf Terme zweiter Ordnung aufge- 
hoben. Daher wird die in diesem Koordinatensystem gemessene Periode der 
Uhr, die mit ihrer Eigenperiode identisch ist, in der ersten Näherung von dem 
Gravitationsfeld unabhängig. 

Geht man nun zur Rotverschiebung zurück, so sieht man gleich, daß es 
sich hier um die Perioden bzw. die Wellenlängen der Spektrallinien handelt, 
die von einem Atom emittiert werden. Im Atom selbst gibt es aber keine 
Koinzidenzen, die sich mit diesen Perioden wiederholen, so daß bei der Be- 
handlung der Rotverschiebung eine zu enge Anlehnung an die der Normaluhr 
nieht gerechtfertigt ist. Trotzdem kann man die von einem im Gravitations- 
feld frei fallenden Atom emittierten Spektrallinien ähnlich wie die makrosko- 
pische Uhr mit Hilfe des Äquivalenzprinzips behandeln und auf diese Weise 
die Einsteinsche Formel für die Rotverschiebung erhalten. Man wird aber 
noch bedenken müssen, daß die z. B. am Sonnenspektrum beobachtbare Rot- 
verschiebung sich auf Atome bezieht, die nicht frei im Gravitationsfeld fallen, 
sondern die mittlere Beschleunigung Null haben. Schon aus diesem Grunde 
wäre eine direkte Berechnung der Rotverschiebung unter genauer Berück- 
siehtigung der Atomstruktur erwünscht. Wichtiger noch ist aber die Tat- 
sache, daß durch eine solche Rechnung die innere Kohärenz der Theorie sowie 
die gegenseitige Verträglichkeit der theoretischen Behandlung der verschiede- 
nen Teilprobleme, in welche sich die Rotverschiebung aufspalten läßt, ge- 
prüft wird 3), 

Die unmittelbare Berechnung der Rotverschiebung ist das Thema der 
vorliegenden Arbeit. Den Ausgangspunkt liefert die Bemerkung, daB es sich 
auch bei der Rotverschiebung im wesentlichen um ein Bewegungsproblem 
handelt: Ein Elektron bewegt sich um den Atomkern. Gegenüber den beiden 
ersten Effekten kommen aber jetzt zwei neuartige Umstände hinzu, welche 
das Problem bedeutend komplizierter machen. Zunächst hat man es mit 
einer Bewegung im kombinierten Gravitations- und elektromagnetischen Feld 
zu tun, und zwar in dem elektromagnetischen Feld eines im (makroskopischen) 
Gravitationsfeld der Sonne befindlichen Atomkerns. Ferner handelt es sich 
im vorliegenden Falle um eine quantisierte Bewegung, zu deren Unter- 
suchung also auch die Quantentheorie berücksichtigt werden muß. 

Die genaue Berechnung des kombinierten Gravitations- und elektro- 
magnetischen Feldes wäre nach Einstein nur bei Verwendung der einheit- 
lichen Feldtheorie möglich. Da die endgültige Formulierung dieser Theorie 


3) Ähnlich sind die Verhältnisse bei der Periheldrehung und der Lichtablenkung. 
Im wesentlicher folgt nämlich das geodätische Bewegungsgesetz unmittelbar aus dem 
Aquivalenzprinzip. Trotzdem bildet die Ableitung dieses Gesetzes aus den Feldgleichungen 
ein wichtiges und lehrreiches Ergebnis. 


Fl 
a 
E 
- 
4 
a 


.. Annalen der Physik. 6. Folge. Band 17. 1956 wie q 

noch nicht bekannt ist, kann man die Rotverschiebung nur provisorisch be- 
handeln, indem man diejenigen Feldgleichungen verwendet, welche sich ays 
der bekannten, additiv zusammengesetzten Lagrange-Funktion des kom- 
binierten Feldes ergeben. Andererseits läßt sich die Quantisierung der Be- 
wegung ohne weiteres berücksichtigen, wenn man die für einen Riemann. 
schen Raum verallgemeinerte Dirac-Gleichung benutzt*). Dabei hat man 
offenbar für das Gravitationsfeld die Schwarzschildsche Lösung einzu- 
setzen, da es sich um ein makroskopisches Feld handelt, auf welches der 
betrachtete Atomkern einen völlig unbedeutenden Einfluß hat. Es sei noch 
darauf hingewiesen, daß die Verwendung des Schwarzschildschen Gra- 
vitationsfeldes bei den folgenden Rechnungen in einer anderen Hinsicht eine 
entscheidende Rolle spielt. In diesem statischen, kugelsymmetrischen Feld 
ist nämlich eine Zeitskala physikalisch ausgezeichnet. Der Existenz dieser 
ausgezeichneten Zeitskala ist es zu verdanken, daß man die stationären Atom- 
zustände und die Energie eines Probeteilchens ohne weiteres definieren sowie 
die Energien von Probeteilchen an verschiedenen Raumpunkten unmittelbar 
miteinander vergleichen kann. 

Die Behandlung der verallgemeinerten Dirac-Gleichung wird besonders 
einfach, wenn man die isotrope Form des Schwarzschildschen Linien- 
elements verwendet) : 


ds? = — a (r) (dx? + dy? + dz?) + y (r) dt?; Bi 
m\4 m \2 
a =(1 ‘ y= (i -5,) (1 + =) (1,2) 


rechnung der Rotverschiebung in der Dirac-Gleichung die Größen x und y 
als konstant ansehen, da sie sich innerhalb des Atomvolumens äußerst wenig 
ändern. Durch die Vernachlässigung der Ortsabhängigkeit von « und y wird 
auch eine weitere, für die folgenden Betrachtungen notwendige Annahme in 
besonders einfacher Weise eingeführt. Es handelt sich um die Aufhebung der 
Gravitationsbeschleunigung. Konstante Werte von « und y bedeuten näm- 
lich eine verschwindende Gravitationsbeschleunigung. Daher dürfen wir 
bei der Rechnung den Atomkern als ruhend annehmen®). Diese Rechnung 
wird im letzten Abschnitt der vorliegenden Arbeit durchgeführt. Vorher 
wollen wir im nächsten Abschnitt eine einfachere Methode anwenden, die zu 
einer besonders anschaulichen Beschreibung der Rotverschiebung führt. 
Br Es handelt sich um die relativistische Behandlung des Bohrschen Atom- 
a) 2 modells nach Sommerfeld, die sich für den Fall, daß es ein Gravitationsfeld 
von der Form (1,1) gibt, unmittelbar verallgemeinern läßt. Wie wir sehen 


4) E. Schrödinger, S.-B. Berl. Akad. 105 (1932); W. Pauli, Ann. Physik 18, 
337 (1933); L. Infeld u. B. L. v. d. Waerden, S.-B. Berl. Akad. 380 (1933). 

5) Vgl.z.B. G. Möller, The Theory of Relativity, Oxford (Clarendon) 1952; § 123. 

6) In Wirklichkeit wird die Gravitationsbeschleunigung durch die Bildung eines 
Druckgradienten in der Sonnenatmosphäre ausgeglichen. Die Zusammenstöße zwischen 
den Atomen, auf welche der Druck nach der kinetischen Gastheorie zurückgeführt wird, 
sind für die so schwer zu berechnende Druckverschiebung der Spektrallinien verant- 
wortlich. — Es sei noch bemerkt, daß es einen — allerdings nur theoretisch interessanten — 
Fall gibt, in welchem « und y streng ortsunabhängig sind. Gäbe es nämlich innerhalb des 
das Schwarzschildsche Feld erzeugenden kugelförmigen Körpers eine konzentrische 
Hohlkugel, dann wären in dieser Kugel « und y streng konstant. Unsere Rechnung würde 
also für Atome, die sich innerhalb einer solehen Kugel befinden, streng gelten. 


werden 
Reehnur 


Die | 
und elek 
lautet 


wobei 
seine Gt 


und 


bedeute 
hewegel 
stimme! 
Die Ha 
kanonis 
lassen s 


wobei q 


bedeute 
pulse di 
die Ene 


Man be 
und de: 


nach F 
Es | 


u 
&. 
U 
ab 
Bs 
er 
Wie störungsrechnerische Betrachtung zeigt, darf ie Be- 
daß H 
Feld k 
Age RER 
E 
sich ab 
514 (19: 
Shiwa 
schlager 
Bf 8) 
we der Fal 
wir uns 
nischen 
6. Fal 
vir 7 
— 


4 A. Papapetrou: Rotverschiebung und Bewegungsgleichungen 217 


werden, führt auch diese einfache, einem spinlosen Elektron entsprechende 
Reehnung zu demselben Ergebnis wie die Rechnung mit der Dirac-Gleichung’). 


§ 2. Das Bohrsche Atommodell 


Die Bewegungsgleichung eines geladenen Probeteilchens im Gravitations- 
und elektromagnetischen Feld, die sich aus der additiven Feldtheorie ergibt, 


lautet 

c + mc = 0, (2,1) 
wobei m, und —e die Ruhmasse und Ladung des Teilchens, u“ = dzt/ds 
vine Geschwindigkeit, F,, = — Pu,» die elektromagnetische Feldstärke 
und 

bedeuten. Das uns interessierende Probeteilchen ist das sich um den Atomkern 
bewegende Elektron. Die Energie dieses Teilchens läßt sich eindeutig be- 
stimmen, wenn man die Bewegungsgleichung in kanonischer Form schreibt: 
Die Hamiltonfunktion ist dann mit der Teilchenenergie identisch. Die 
kanonischen Gleichungen, die dem Bewegungsgesetz (2,1) gleichwertig sind, 
lassen sich unmittelbar hinschreiben®). Setzt man 


P, = Myc Uy, + Pu 


Un = Gyro ur (2,3) 


wobei p, das elektromagnetische Potential und) 


bedeuten, so sind die drei ersten Komponenten von P, die kanonischen Im- 
pulse des Teilchens. Dagegen liefert ce P, die Hamiltonfunktion, also auch 
die Energie des Teilchens: 
E=H =cP,=mcu + eq. (2,4) 
Man bestimmt H als Funktion der kanonischen Impulse, der Koordinaten 
und der Zeit, wenn man die Gleichung j 


Pu) (P, ay 9») = 
nach P, auflöst. 


Es folgt aus der allgemeinen Beziehung wer eee 
( 


dH _ 0H 

a 
daß H in einem zeitunabhängigen Gravitations- und elektromagnetischen 
Feld konstant bleibt. Dieser Fall liegt bei der Bewegung des Elektrons um 


2,4a) 


’) Es sei noch auf zwei frühere Berechnungen der Rotverschiebung hingewiesen, die 
sich aber auf andere Gravitationstheorien beziehen: Moshinsky, Phys. Rev. 80, 
5l4 (1950), hat die Birkhoffsche Gravitationstheorie zugrunde gelegt; Belinfante u. 
Shiwart, Report of Research, Purdue University 1954, betrachten die von ihnen vorge- 
schlagene lineare Gravitationstheorie. 

8) Vgl. A. Papapetrou u. W. Urich, Z. Naturforschg. 10a, 109 (1955). Dort wurde 
der Fall des allgemeinen Pol-Dipolteilchens behandelt. Das einfache Polteilchen, für das 
wir uns hier interessieren, ergibt sich daraus, wenn man f = m, = const setzt. Die kano- 
nischen Bewegungsgleichungen des einfachen Polteilchens wurden schon früher von 
6. Falk, Z. Physik 182, 44 (1952), jedoch in vierdimensionaler Form, abgeleitet. 
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den Atomkern vor. Wir müssen zunächst die Energiewerte der stationären 
Bewegungszustände des Elektrons bestimmen. Daraus können wir dann dure 
Subtraktion die Energiewerte der vom Atom emittierten Photonen im Augen- 
blick ihrer Emission ableiten. Nach der Emission bewegen sich die Photonen 
vom Emissionsort weg. Die Photonenergie wird ebenfalls durch eine Formel 
der Form (2,4) gegeben, wobei aber diesmal e = 0 zu setzen ist, da es sich um 
ein neutrales Teilchen handelt. Da sich nun die Photonen in einem statischen 
Gravitationsfeld bewegen, bleibt ihre Energie nach (2,4a) während der Be. 
wegung konstant. Dieses Ergebnis zeigt, daß die Rotverschiebung nur ans 
einer Änderung der Energiewerte der stationären Atomzustände folgen kann. 
Dagegen läßt sie sich nicht mit Hilfe einer potentiellen Gravitationsenergie 
des Photons erklären?). 

Wir berechnen nun das elektrostatische Potential eines Protons, welche 
sich im Gravitationsfeld (1,1) befindet, unter der Annahme, daß in der für 
das Atom interessierenden Umgebung des Protons x und y konstant sind. 


Die Maxwellsche Gleichung lautet oie 
== 3". (2,5) 


FH = = — 9): 


®) Einen solchen Erklärungsversuch findet man bei W. H. Mc Crea, Phil. Mag. 4 
1010 (1954). Im Schwarzschildschen Gravitationsfeld kann man eine potentielle Energi 
des Probeteilchens definieren. Man definiert nämlich zunächst die kinetische Energie des 
Teilchens, indem man es durch einen passenden Stoß zur Ruhe gebracht denkt: 


=E-E, 


wobei E und £’ die Energie des Teilchens vor und nach dem Stoß bezeichnen. Wenn die 
Ruhmasse des Teilchens beim Stoß unverändert bleibt, wird Z’ durch 


E' = mc? ug = myc Is, 


gegeben, wobei u“ die Geschwindigkeit nach dem Stoß bedeutet. Es ist also uw’! = u?= 
w®= 0 und daher 1=g,,u*u” = gy, (u. Daraus folgt 


E'= m V 908 = M, c? Vy. 


Es ist aber in erster Näherung Yy = 1 — >= 1 — > Py (@ die Newtonsche Gra 
tionskonstante und M die Sonnenmasse) und daher 
Man kann also schreiben hy 
M m, E 
= E— m. 


Da die rechte Seite dieser Beziehung eine konstante ist, kann der Term —@ M my/r ähnlich 
wie in der Newtonschen Gravitationstheorie als potentielle Energie gedeutet werden. 
Es ist aber zu beachten, daß dieser Ausdruck der Ruhenergie des Teilchens proportional 
ist und daher beim Photon gleich Null wird. 
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Die Strom-Ladungs-Dichte 5“ hat im vorliegenden Falle nur eine nicht-ver- 
schwindende Komponente, nämlich 34 = o. Wir dürfen also gy, = 0 fir k = 
1,2,3 setzen, wodurch sich (2,5) auf folgende einzige Gleichung reduziert: 


(= — 9) 
pa nach (1,1) V-g =] asy ist, kann man diese Gleichung auch so schreiben: 


2, 


das elektrostatische Potential des Protons ohne Gravitationsfeld bedeutet. 
Wir setzen (2,6), (2,7) in (2,4) ein und erhalten für die Energie des Elektrons 
die Formel 

/ye 


Im folgenden wollen wir nur die kreisférmigen Bahnen betrachten. Es 
sei also 


x«=Reoswt, y=Rsinwt, 2=(. (2,9) 
Die Bewegungsgleichung (2,1) liefert dann eine einzige Bedingung”): 
2 ps — e? 
; c dt R?\- 
mit ut =. = (y—a ) (2,11) 


Rine weitere Beziehung wird durch die Quantenbedingung geliefert. Diese 
hat die einfache Form 
J, (2,12) 


dabei sind J, der Drehimpuls des Teilchens und / eine ganze Zahl!!). Der 
Ausdruck für J, ergibt sich eindeutig aus der Integration der Bewegungs- 
gleichung. Es gilt nämlich wegen der Kugelsymmetrie des Feldes um den 
Atomkern die Beziehung 


J, =—2«xP,+ y P, = const. (2,13) 
Aus (2,12) und (2,13) ergibt sich für die Bewegung (2,9): eb 
am, oR? u = 1h. (2,14) 


) Der zweite Term in (2,1) verschwindet infolge der Annahme &,y = const. 

1) Die rechte Seite von (2,12) muß als Phasenintegral eine invariante Größe sein 
und kann daher in der allgemeinen Relativitätstheorie nicht vom Gravitationsfeld ab- 
hängen. Nur in einer Theorie mit zwei Tensoren 9,,und ut vgl. z. B. N. Rosen, 
Phys. Rev.57, 147 (1940) — könnte man eine feldabhängige skalare Größe bilden und sie 
eventuell in die rechte Seite von (2, 12) einzuführen versuchen. 
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Und wenn wir dies in (2,8) einsetzen: at ree 


(2,15) 


ea 
£,=Vy 1 za) 
Diese Energiewerte unterscheiden sich von den entsprechenden Sommer. 
feldschen lediglich um den Faktor Vy . Daher werden auch die Energien 
der vom betrachteten Atom emittierten Photonen um den Faktor Vy kleiner. 
Da nun die Energie der Photonen während ihrer weiteren Bewegung nicht ge. 
ändert wird, so folgt, daß Atome, die sich an verschiedenen Orten innerhalb 
des Schwarzschildschen Gravitationsfeldes befinden, Photonen emittieren, 
deren Energie den an den entsprechenden Emissionsorten herrschenden 
Werten von Vy proportional ist. Dies entspricht aber gerade dem Einstein- 
schen Wert der Rotverschiebung. 


2 § 3. Die Diraesche Gleichung 


‘Der Rechnung wollen wir die von Infeld und v. d. Waerden gegebene 
Form der verallgemeinerten Dirac-Gleichung mit dem zweidimensionalen 
Spinraum zugrunde legen2). Der Darstellung von Infeld und v. d. Waerden 
folgend, werden wir in diesem Abschnitt die Tensoren des vierdimensionalen 
Raumes mit lateinischen, die des Spinraumes dagegen mit griechischen In- 
dizes bezeichnen. Die verallgemeinerte Dirac-Gleichung lautet: 


(*, —y? x 


imyc fr 
2) 
Ak @ 


ii ye = 0, riz 


A 
(ye + very’) + “ 


Der gemischte Tensor o*4 genügt der Beziehung 


wobei y, » der Fundamentaltensor des Spinraumes ist. Die zwei Spinvektoren Dat 
y* und x; stellen die (vierdimensionale) Wellenfunktion des Elektrons dar. f die 
Yar sind die Komponenten des affinen Zusammenhangs im Spinraum und 
sind aus folgenden Beziehungen zu berechnen: 


Ger — yur, 


1 A + yt gkiv — 
+H ¢ o 0, 


ri i. sind die Christoffelschen Symbole des Riemannschen Raumes. a 

Wir wählen das Koordinatensystem im Spinraum derart, daß m 


Ya=- Ya = 1 
‘Vel. Anmerkung *). 


Au let 
N ‘3 
‘vs 
Bei 
schi¢ 
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(2,15) 
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wird. Dann ergibt sich als die einfachste Lösung von (3,2) für die Metrik (1,1): 
; 1 : 
(ot4¥),. 


okiu (oki), für k= 1, 2, (3,5) 


(#%#), bedeutet eine Lösung im Minkowski-Raum. Im folgenden wollen 
wir die von Infeld und v.d. Waerden angegebene Lösung (in), ver- 
wenden : 


Die allgemeine Lösung von (3,3) läßt sich leicht bestimmen und lautet: 
1 


Bei der Metrik (1,1) sind die folgenden Christoffel-Symbole von Null ver- 
schieden : 


; 
un 


1 . 


‘vs 


= = =! = = 557° 
11 22 13 23 - en. 8) 

, 


Dabei bedeutet der Strich die Ableitung nach r. ,7) ergeben sich dann 


die Werte von 


Aus (3 


1 
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Es genügt nun, (3,9) in (3,1) einzuführen, um die Dirac-Gleichung in 
ihrer endgültigen Form zu erhalten. Es ist aber zweckmäßiger, daß man 
zunächst die Gleichungen (3,1) mit Hilfe vierdimensionaler Matrizen zusammen- 
faßt. Dies läßt sich ohne weiteres erreichen, wenn man folgende Matrizen 
einführt: 


0—1 0 0 0— 0 0 100 9g 
0 00 17°% 10 0 0 001 OF 
8 r 6 0 0 0009-1 
0100 0001 9 
% 0 0 1 07°%-11 0 0 OF? 0 
( w ; 
0 0 0 l ) l ( ) 0 0—y} —pB, 
und 
x2 


Die Dirac-Gleichung nim mt dann die Form an 000 

| 3 a, é iM, 

Multiplizieren wir diese Gleichung mit . Yy, so erhalten wir die üblichere 


Form 


mit 
he y 

H V Xo (= Vy — 

Benutzt man nun die Werte (3,9) fiir Yan so ergibt sich 2 


x’ 3 e 


Daraus folgt ee 


h 3 \ ay’ 
== VZ i | 
+ e€6,9,+ Vy Mo Gy. 


Sind U und V zwei beliebige Spinoren der Form (3,11), so folgt aus (3,15) 
die Beziehung 


he y/y BY x 
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4. Papapetrou: Rotverschiebung und Bewegungsgleichungen 


Danach wäre der Operator H nur im Falle x, y = const selbstadjungiert. 
Multiplizieren wir aber diese Beziehung mit a°:, so wird 


3 

£ (U* HV — V A* U*) =" E(Yya U* oa, 

1=1 

d.h. der Operator «’/:H ist selbstadjungiert. H ist der Energieoperator und 
E = [ H YdV (3,16) 

der Energiemittelwert für den durch VY beschriebenen Zustand. Der Faktor x": 


ist bei der Mittelwertbildung aller Operatoren nötig. Man sieht das unmittelbar, 
wenn man den Ausdruck für die Wahrscheinlichkeitsdichte betrachtet. Diese 


ist der dureh e dividierten vierten Komponente der Stromdichte 


: 


gleich. Beachtet man (3,5), (3,6) und (3,11) sowie den aus (1,1) folgenden 
Wert / g= Vy x3, so ergibt sich fiir die Wahrscheinlichkeitsdichte der 
Ausdruck x’: YW*,. Der Mittelwert eines Operators A wird also durch ee 


A = [ Pr 
gegeben. 
Wir kommen jetzt zur Berechnung der Rotverschiebung. Beachtet man, 
daß in (3,15) r den Abstand des Atoms vom Mittelpunkt der Sonne darstellt 
und daher von der Größenordnung des Sonnenradius sein wird, so ist es leicht 
einzusehen, daß in (3,15) die den Faktor x’ bzw. y’ enthaltenden Terme um 
einen Faktor von der Größenordnung Atomradius/Sonnenradius kleiner als 
die für die Rotverschiebung interessierenden Terme sind. Man darf also diese 
Terme ohne weiteres weglassen, d.h. für x und y die am Ort des Atomkerns 
herrschenden Werte einsetzen, wie wir es schon in der Einleitung bemerkten. 
Beachtet man noch, daß 9, = %, = Ya = 0 ist, so vereinfacht sich (3,15) zu 
/ 3 


ı=1 


Wir interessieren uns für die stationären Bewegungszustände des Elektrons, 
welche durch die Beziehung 


= yp (a, y, 2) - e-iEtih 


charakterisiert sind. Dann nimmt die Dirae-Gleichung die Form an ee 
(E—H)y= 

Führen wir in (3,17) den Ausdruck (2,6) für das elektrostatische Potential des 

Protons ein, so erhalten wir die Dirac-Gleichung des Wasserstoffatoms in 

der Form 


Diese Gleichung unterscheidet sich von der normalen (d.h. von der im 
Minkowski-Raum geltenden) Dirac-Gleichung lediglich um die konstanten 
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Faktoren yz und Vy vor den drei letzten Termen. Man kann sie auf die 


Normalform zurückführen, wenn man die neuen Konstanten 


(hc)’ = y2 hc, = me = Vy mg (3,19) 


einführt. Die Eigenwerte der normalen Dirae-Gleichung lauten Zu 
E, (n, j) =merlı + My 
o (”; 7) 0 (n +VP- (3,20) 


Man erhält daraus die Eigenwerte der Gleichung (3,18), wenn man fic, & 
und m, durch (Ac)’, (e2)’ und m, ersetzt. Beachtet man, daß nach (3,19) 
(e2)’/(h c)’ = e2/hc ist, so findet man für die Eigenwerte von (3,18) die Formel 


E (n, j) = Vy By (n, j)- (3,21) gleich 
In der hier betrachteten Näherung werden also alle Energiewerte des Atoms Spezi 


mit dem Faktor Vy multipliziert!3). Wie am Ende des vorigen Abschnitts 
gezeigt wurde, folgt daraus der Einsteinsche Wert der Rotverschiebung. 


Für n = 0 geht die Formel (3,21) mit j—1 in (2,15) über. 


Ei 

Berlin, Forschungsinstitut für Mathematik der Deutschen Akademie der fallbe 
Wissenschaften. on 
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Zur Theorie des Kathodenfalles einer Glimmentladung 


Von Johannes Wilhelm 


Inhaltsverzeichnis 


Ausgehend von gewissen Näherungsannahmen wird eine Differential- 
gleichung aufgestellt, die den örtlichen Verlauf des Potentials beschreibt. 
Spezielle Lösungen ergeben räumliche Abhängigkeiten des Feldes, wie sie 
durch experimentelle Untersuchungen gefunden wurden. => 


Eine kritische Betrachtung der physikalischen Gegebenheiten im Kathoden- 
fallbereich einer Glimmentladung!) zeigt die besondere Schwierigkeit, welche 
sieh bei der theoretischen Erfassung des Problems ergibt: 

Sie liegt nämlich im wesentlichen darin, daß es sich um ein Gebiet handelt, 
wo ein inhomogenes elektrisches Feld vorliegt, welches von relativ hohen Feld- 
stärkewerten in Nähe der Kathodenoberfläche bis zu verschwindend geringen 
in der Gegend des Glimmlichtes variiert. Dementsprechend ergibt sich, daß 
einige zur Beschreibung der Verhältnisse wichtige funktionale Zusammen- 
hänge, welche für den Fall homogener Felder abgeleitet worden sind, nur noch 
bedingte Gültigkeit besitzen. Das betrifft vor allen Dingen die Ausdrücke für 
die Abhängigkeit der Ionisation vom Feldverlauf sowie die Darstellung der 
Trägerbewegung im elektrischen und im Konzentrationsfeld. Die im allge- 
meinen für diese Größen existierenden Funktionen wurden nur in hinreichend 
homogenen Feldern bestätigt, so daß eine einfache Übertragung auf das 
Kathodenfallgebiet nicht ohne weiteres gerechtfertigt erscheint. 

Erste Beiträge zur Erfassung der Ionisation in inhomogenen Feldern 
sind durch die Untersuchungen von Morton?) und Johnson?) gegeben. 
Andererseits wurde bezüglich der Trägerbewegung bereits von R. Seeliger) 
anläßlich seiner Untersuchugen über die Dicke des Fallraumes darauf hinge- 
wiesen, daß im Kathodenfallbereich wahrscheinlich die theoretische Behand- 
lung eines Gemisches von Freifall und Driftbewegung erforderlich ist. Ein 
Weg hierzu wird durch die Untersuchungen von Höcker und Bez°) vor- 
gezeichnet, wenn man ihre Ansätze sinngemäß auf das Kathodenfallgebiet 
erweitert, in dem man inelastische Stöße in der Energiebilanz berücksichtigt. 


1) Fine ausführliche Diskussion hierzu wurde in Universität Greifswald, 

math. nat. R 5, 495 (1954/55) durchgeführt. OTe. © 
2) P. L. Morton, Physic. Rev. 70, 358 (1946), 
3) G.W. Johnson, Physic. Rev. 78, 284 (1948), ging 
4) R. Seeliger, Z. Naturforschg. 8a, 1, 74 (1953). Kes 


5) K.-H. Hicker u. W. Bez, Z. Naturforschg. 9a, 1, 64 (1954). re zei Er Er; = 
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Be bekannter Feldabhängigkeit der Wirkungsquerschnitte sowie der bei 
einem Stoß abgegebenen Bruchteile an Energie ist damit die Möglichkeit ge- 
geben, die wichtigsten Kenngrößen des Kathodenfallgebietes zu errechnen. 
Allerdings sind dabei kaum geschlossene Lösungen zu erwarten, vielmehr 
dürfte die Lösung in einer numerischen Rechnung bestehen. Vorliegende 
Arbeit beschränkt sich daher auf die Aufgabe, mit Hilfe plausibler Näherungs- 
annahmen eine Differentialgleichung für die örtliche Abhängigkeit des Po- 
tentials aufzustellen, wobei spezielle Lösungen dieser Gleichung Feldver- 
teilungen in Längsrichtung der Entladung liefern, welche mit den in der 
Literatur vorkommenden vereinbar sind. Nun wurde die Bestimmung des 
Feldverlaufs bisher nach drei experimentellen Methoden durchgeführt: 

I. durch Anwendung von Sonden, 2. mittels Querkathodenstrahl und 
3. unter Benutzung des Starkeffektes®). 

Eine Übersicht über die vorhandenen Messungen zeigt dabei, daß zwei 
unterschiedliche Feldverteilungskurven vorliegen: Einerseits vollständig grad- 
liniger Längsabfall über den ganzen Bereich des Kathodenfalles, andererseits 
aber durchhängende konvex zur Koordinatenachse hin gekrümmte Vertei- 
lungskurven (abgesehen von selten vorkommenden konkav gebogenen, wie 
sie z. B. von Warren a. a. O. bei Messungen im anomalen Bereich angegeben 
wurden), wobei bei letzteren noch zu unterscheiden ist, ob es in Kathoden- 
nähe zur Ausbildung eines Feldmaximums kommt, oder ob kein solches vor- 
handen ist. Zur Aufstellung einer Differentialgleichung für die Potential- 
verteilungsfunktion wird nun im folgenden die Trägerbewegung durch einen 
Drift- und Diffusionsstrom beschrieben, wobei die hierin auftretenden Beweg- 
lichkeits- und Diffusionskoeffizienten infolge der Inhomogenität des Feldes 
noch Funktionen des Ortes sind. Wegen der Unsicherheit in der Kenntnis 
dieser Funktionen verzichten wir von vornherein auf eine explizite Angabe 
dieser Größen. Damit lassen sich zunächst die Stromdichtegleichungen in der 
folgenden Form angeben: of 


j, =e D, grad n, + eb, n, &, ee: 


Nun beschränken wir uns bei den weiteren Betrachtungen auf ein ebenes 
Problem, legen also eine Längskoordinate z in die Mittellinie der Entladung 
und führen eine weitere Koordinate x in seitlicher Richtung ein. Spaltet man 
dann nach Komponenten auf, so erhält man: 


= —eD,,gradn,,„teb,,n 


v 


ön,, 
I+v,x e b;, n,, E,. 3 (1b) 


‘wid 


Ganz entsprechend ergibt sich fiir die Elektronen: a innate 


i. =e +eb,n,E,, ar (2a) 
lea = @ (2b) 


8) G. Mierdel, Handbuch d. Experimentalphys. Bd. 13. 3. Teil, S. 363ff. Neuere 
Arbeiten vor allem: W. H. Ernst, Helv. phys. Acta 8, 381 (1935); R. P. Stein, Physic. 
Rev. 89, 134 (1953); R. Warren, Physic. Rev. 98, 1650 (1955). 
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Hierzu kommt noch die Poissongleichung: 
AV =— — me) aba (3) 


die zusätzlichen Beziehungen: 


= Ve Ney = Vy M1 (5) 
mit } 

Ein Vergleich mit den Kathodenfalltheorien, wie sie in der Hauptsache 
von Rogowski und von Weizel, Rompe und Schön’) entwickelt wurden, 
zeigt den Unterschied gegenüber den bisherigen Grundannahmen: 

Es wird durch Hinzunahme einer zweiten Koordinate der seitlichen Aus- 
dehnung der Entladung Rechnung getragen; die Einwirkung der Inhomogeni- 
tät des Feldes auf Beweglichkeit und Diffusionskoeffizienten wird — wenn 
auch ohne definitive Festlegung — berücksichtigt, mehrere Ionenarten werden 
zugelassen und schließlich soll auf eine explizite Angabe der Ionisations- 
verhältnisse verzichtet werden. 

Ausgehend von den eben angegebenen Beziehungen gelangt man zur Auf- 
stellung einer Differentialgleichung für die Potentialverteilung, wenn man 
eine Reihe von Voraussetzungen vornimmt: 

Zunächst soll hinsichtlich der Kontinuitätsgleichung angenommen werden, 
daß innerhalb wenigstens eines kleinen Bereiches um die Rohrachse herum 
mit annähernd konstanter Entladungsstromdichte in bezug auf die Längs- 
riehtung gerechnet werden darf. Dann spaltet sich (4) innerhalb dieses Gebietes 
in die beiden Gleichungen 


auf, wobei (4b) durch Integration unter Berücksichtigung von j, (z, 0) = 0 
folgt. Weiterhin soll eine Änderung der Größen y, d.h. der Verhältnisse der 
verschiedenen Trägerdichten zueinander innerhalb des eben gekennzeichneten 
Bereiches um die Rohrachse herum in seitlicher Richtung vernachlässigt 
werden. Damit ergibt sich durch Kombination der Beziehungen (1b), (2b), 
(3) und (4b) sowie (5) 


AV =I'VAV (6) 


als partielle Differentialgleichung für die Potentialverteilungsfunktion V (z, x). 
Hierbei bedeuten die Punkte partielle Ableitungen nach der x-Koordinate 


und J” die Abkürzung für das Verhältnis Sara 


D, D; v V+ v 
Eine letzte anzunehmende Näherung soll nun bezüglich dieser eben neu 
eingeführten Größe darin bestehen, daß die Änderung der Beweglichkeiten und 
Diffusionskoeffizienten in seitlicher Richtung gegenüber der axialen vernach- 


7) W.Rogowski, Arch. Elektrotechn. 26, 643 (1932); W. Weizel, R.Rompe u 
M. Schön, Z. Physik 112, 339 (1939). 
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lässigt bzw. weniger gefordert, daß der Ausdruck J” näherungsweise von der 
seitlichen Richtung unabhängig angesehen wird. Spezielle Lösungen der La 
werden dann dadurch gewonnen, daß sowohl mit einem Produktansatz V (2, x) 
h(2) g(x) als auch mit V(zx) = h(z) + g(x) in die Ditferentialgleichung 
eingegangen wird. Benutzt man z. B. einen cdditiven Ansatz der genannten 
Form, so erhält man zunächst die Beziehung 

rg th" +9) 
Die dieser Gleichung genügenden Funktionen h(z) und g(x) müssen auch die 
einmal nach x diffenrenzierte Gleichung erfiillen. Bei entsprechender Division 
der von x abhängigen Koeffizienten kann es aber immer erreicht werden, 
daß jeweils eins der drei in (8) vorkommenden Glieder nur noch von der Koor- 
dinate z abhängig bzw. überhaupt eine Konstante ist, so daß dieses Glied bei 
der nachfolgenden Differentation aus der neu entstehenden Gleichung ganz 
herausfällt. Da hierbei von vornherein Lösungen der Art e = 0 und V = Va) 
sowie h’’ = 0 als physikalisch nicht brauchbar ausgeschlossen werden können, 
ergibt sich in jedem Falle, ganz gleich wie die Division durchgeführt wird, 
daß die Größe J” und auch h” Konstante sein müssen. Damit erhält man 


die folgenden Beziehungen 


hk” =c 


Die Verwendung eines Produktansatzes liefert dagegen; TE 
+ BR'+CrH+Dh=0, (10 
wobei die Abkiirzungen 


A=gg B=-g C=9g9 D=-g 


verwendet worden sind. Das fiihrt durch analoge, den obigen cutapenthallll 


wegen der der Mehrzahl der Glieder jetzt aber auch REN. Über- 
legungen auf die folgenden Differentialgleichungen: 


g —k*g —-kk*ggtkg=0 (AL) 
Th=kt. (116 


Eine vollständige Auswertung dieser Beziehungen würde Annahmen 
über die Größe der in den Gleichungen vorkommenden Konstanten c, «* 
sowie k und k* erforderlich machen, die aus entsprechenden Rand- bzw. An- 
fangsbedingungen zu entnehmen wären. Für unsere Zwecke genügt es, sich 
auf den z-abhängigen Anteil der Potentialverteilungsfunktion zu beschränken. 
Zunächst ist die durch einen additiven Ansatz erzeugte spezielle Lösung so 
gebaut, daß ihre zweite partielle Ableitung nach z konstant wird. Bei passender 
Bestimmung der auftretenden Konstanten aus den Rand- bzw. Anfangs- 
bedingungen ist also eine solche Lösung — wie ohne weiteres ersichtlich — 
in der Lage, linearen Feldabfall in Längsrichtung (wie er z. B. unter den 
Meßergebnissen von Stein a.a. O auftritt), zu beschreiben. Dementgegen 
eignet sich die weitere spezielle Lösung gut dazu, die durchhängenden Feld- 
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stirkekurven wiederzugeben. Hierzu ist von der ersten Beziehung des Glei- 
chungssystems (11) auszugehen. Negatives k führt auf die Überlagerung 
einer cosinus- und sinus-Funktion, eine Lösung, die wir im weiteren nicht 
diskutieren wollen. Dagegen liefert positives k = x? eine reelle Lösung der 
Form: 


h= hg | oj x z + Sin x2] 


wobei h, und A, Integrationskonstante bedeuten. Bei entsprechender Multi- 
plikation mit einer Konstanten und Festlegung des Potentialnullpunktes auf 
V(0, 0) = 0 ergibt sich hiermit für V (z, 0) die Beziehung: 

V (z, 0) = A [(Gojxz+ BSinxz— 1]. (12) 
Die hierin auftretenden-Konstanten A, B und x müssen aus Rand- bzw. An- 
fangsbedingungen festgelegt werden. Bezeichnet man mit E, die Größe 
B,=— V’ (z,0) und sieht man als Grenze Kathodenfall-Glimmlicht die- 
jenige Stelle an, wo die Längsfeldstärke verschwindet, so hat man folgendes 
Gleichungssystem zur Bestimmung der erforderlichen Größen: 


E, =—xAB, (13 a) 
BGojxd)=0, (13h) 
= A (Cof xd + BSinxd — 1]. 
Hieraus folgt: 
(14) 


Letztere Gleichung stellt einerseits bei vorgegebenem V,, E, und d eine 
transzendente Bestimmungsgleichung für x dar; andererseits kann sie aber 
auch als ein Ausdruck für den Kathodenfall V, benutzt werden, wenn man 


nur beachtet, daß x auch durch das Verhältnis x? = > ausgedriickt werden 


kann, d.h. durch Feldstärke und Krümmung derseiben an der Kathoden- 
oberfläche. 

Im folgenden wird ein Vergleich mit Meßwerten durchgeführt, die von 
Ernst a. a. O. bei Untersuchungen mit Glühsonden im normalen Kathoden- 
fallbereich einer Wasserstoff- bzw. Argonentladung gefunden worden sind. 
Da die Meßpunkte nur bis 0,11 em vor der Kathode angegeben werden und 
außerdem eine graphische Differentation entsprechende Ungenauigkeiten in 
sich birgt, wurde an Stelle von (13a) die Beziehung: 

V, =A [oj xz, + BSinxz, — 1) (13 a) 
verwendet wo z, = 0,11 und V, den zugehörigen Potentialwert bedeuten. Als 
V,-Werte wurden diejenigen gewählt, wo das Potential konstant zu werden 
beginnt. Die hiermit verbundene Ungenauigkeit in der Festlegung dieser 
Stellen mußte in Kauf genommen werden. Damit ergeben sich modifizierte 
Gleichungen zur Bestimmung der erforderlichen Konstanten. Die in den Ab- 
bildungen angegebenen ausgezogenen Potentialkurven stellen die nach Fest- 
legung der Integrationskonstanten errechneten Werte dar, die Punkte bzw. 
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Kreuze die entsprechenden Meßwerte. Als Parameter geht dabei der 
in dem Sinne ein, daß die in Kathodennähe höher gelegenen Kurven jeweils 
auch zu höheren Drucken gehören. Während sich die in Abb. 1 angegebenen 
Kurven auf Wasserstoff beziehen, handelt es sich bei Abb. 2 um Messu 

in Argonentladungen. Die in Wasserstoff benutzten Drucke betrugen dabei der 
Reihe nach p = 3,15; 2,60; 2,10 und 1,60; im Argon p = 1,25; 1,09; 0,78 und 
0,624 mm Hg. Es ergeben sich Abweichungen der Meßpunkte von den theore- 
tisch errechneten Kurven, die nieht mehr als 51/,% betragen. In den Ab. 


Qn 02 32 “Ss 06 OF OB 9 oz Qn 2 OF 065 06 
Abb. 1. Potentialverteilung in Wasserstoff Abb. 2. Potentialverteilung in Argon 
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Abb. 3. Feldverteilung in Wasserstoff Abb. 4. Feldverteilung in Argon 
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Druck # bildungen 3 und 4 wurden außerdem die aus den Rechnungen folgenden Ver- 
jeweik # teilungskurven für die Ableitung des Potentials in Längsrichtung angegeben. s 
ebenen Nachdem in dieser Weise ein Vergleich mit experimentellen Daten vorge- _ B 
sungen # nommen worden ist, sind noch einige Bemerkungen zu den Begleitumständen 4 
bei der zu machen, die mit den von uns gewählten speziellen Lösungen verknüpft sind. 
78 und Eine erste Bemerkung bezieht sich hierbei auf die Raumladungsverhältnisse, 


heor- # wozu natürlich die Poissongleichung herangezogen werden muß. Hier- 

n Ab. nach ergibt sich die Raumladung für die Lösungen (9) und (11) beziehungs- 

r Andererseits erhält man durch Ti der Gleichungen (9b) bzw. (11b) 


für die seitliche Abhängigkeit des Potentials die — 
SS. 


| EN A, = A, = — 2c, A, = —)2 [= — -¢ | 
— Io 
rgon B, B, = —k, = — 7 AC} 


k=( + kg) 


d.h. in die Lösung gehen als Integrationskonstante die Größen c bzw. k und 
%=9(0) ein. Dabei stellen aber c bzw. k im wesentlichen den seitlich ab- 
hängigen Anteil der Raumladung in der Rohrachse dar. Zur Integration 
müßte also dieser Anteil als Integrationskonstante vorgegeben werden. Damit 
kann folgendes gesagt werden: Wird im Falle der 1. speziellen Lösung an einer 
Stelle in der Rohrachse eine positive Raumladung vorgegeben, so besteht 
längs der Rohrachse überall dieselbe konstante positive Raumladung. Ähn- 
liches läßt sich auch hinsichtlich der 2. Lösung sagen, nur mit dem Unter- 
schied, daß, wegen der Abnahme des Betrages von h mit zunehmender Längs- 
koordinate, diese positive Raumladung dann ebenfalls auch mit Annäherung 
an das Glimmlicht immer mehr abnehmen wird. 

Während so über die Raumladungsverhältnisse Aussagen gemacht werden 
können, sind solche bezüglich der Trägerverteilung nicht ohne weiteres möglich. 
Hier ergibt sich ein Gemisch von Gleichungen bestehend aus Poisson- 
gleichung, der axialen Stromdichtebeziehung und der aus den Lösungssystemen 
FEN folgenden Gleichungen (9c) bzw. (11e). Eine Ausschöpfung dieser Beziehungen 

4 würde die Kenntnis der funktionalen Abhängigkeit der Beweglichkeiten und 
Diffusionskoeffizienten bzw. je nach der Anzahl der vorkommenden Ionen- 


J > 2 . 
= 
2 
4 


fikation der bestimmenden Differentialgleichungen. 


Differentialgleichung als wenig aussichtsreich. 


wertvolle Diskussionen. 


der Wissenschaften zu Berlin. 
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arten auch noch von Größen y zur Voraussetzung haben. Tatsächlich erweist 
sich die Ableitung von Ausdrücken für die zuerst genannten Größen in einem 
inhomogenen Feld als ein so kompliziertes Problem, daß es den Rahmen dieser 
Arbeit weit übersteigen würde. Es muß daher an dieser Stelle auf eine Nach- 
prüfung verzichtet werden, ob die sich hinsichtlich der Potential- und Raum. 
ladungsverteilung bewährende Lösung auch einer Probe bezüglich der Träger- 
dichteverteilung standhält; jedoch soll nieht versäumt werden, auf diesen noch 
offenen Punkt der vorliegenden theoretischen Untersuchung hinzuweisen. 
Schließlich ist letzten Endes noch eine Bemerkung zur Wahl des Koor. 
dinatensytems zu machen. Von Anfang an war den Betrachtungen ein ebenes 
Problem durch Einführung einer Längskoordinate z und einer seitlichen x 
zugrunde gelegt worden; selbstverständlich wird im allgemeinen eine wirk- 
liche Entladung durch Zylinderkoordinaten wesentlich besser erfaßt werden. 
Wie zu erwarten, hat die Einführung solcher Koordinaten z und r keine 
Änderung bezüglich der für die Längsriehtung gewonnenen Vorstellungen 
zur Folge; lediglich in radialer Richtung ergibt sich eine entsprechende Modi- 


Zusammenfassend kann damit folgendes gesagt werden: Ausgehend von 
gewissen Näherungsannahmen wird eine Differentialgleichung aufgestellt, 
welche die räumliche Abhängigkeit des Potentials beschreibt. Spezielle Lö- 
sungen dieser Gleichungen werden aufgesucht und in ihrer Anwendbarkeit 
hinsichtlich der Wiedergabe des Potential- und Raumladungsverlaufes über- 
prüft; eine Beschreibung der Trägerverteilung selbst scheitert jedoch im 
wesentlichen an der Unkenntnis der funktionalen Abhängigkeit der Beweg- 
lichkeiten und Diffusionskoeffizienten. Auch erscheint eine Aussage über die 
physikalischen Bedingungen hinsichtlich der Auswahl der speziellen Lö- 
sungen aus der Mannigfaltigkeit aller anderen bei der. Kompliziertheit der 


Herrn Professor Dr. R. Seeliger danke ich an dieser Stelle herzlich für 


Greifswald, Institut für Gasentladungsphysik der Deutschen Akadeni 


Bei de der Redaktion eingegangen a am 16. ‘Oktober 1955. 
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Magnetische Pulvermuster auf Kobalt 
bei höheren Temperaturen 


BR Vorläufige Mitteilung 


Von Wilfried Andrä 
Mit 3 Abbildungen 


Inhaltsübersicht 


Mit einer temperaturbeständigen Suspension wird die Veränderung der 
magnetischen Bereichsstruktur von Kobalt zwischen Zimmertemperatur und 
380° C beobachtet. 


Die Richtungsabhängigkeit der magnetischen Kristallanisotropie-Energie 

läßt sich bei Kobalt in der Form 

Ex = K, sin? O + K,sin? © (1) 
schreiben. Dabei ist mit © der Winkel zwischen der hexagonalen Achse und der 
Magnetisierungsrichtung bezeichnet; K, und K, sind Konstanten, deren 
Temperaturverlauf in Abb. 1 angegeben ist. 

Wie aus den Kurven zu entnehmen ist, wird bei Temperaturen oberhalb 
von etwa 250° C die Kristall- 
anisotropie-Energie negativ. 
Die magnetische Vorzugs- 
riehtung geht demnach von 
ihrer Lage parallel zur hexa- 
gonalen Achse in eine dazu 
senkrechte Richtung über. 
Es ist zu erwarten, daß in- 
folgedessen eine andere An- 
ordnung der magnetischen 
Elementarbezirke eintritt als 
bei Zimmertemperatur. 


0°erg cm’? 


Konstonten der Kristallanisotr. -Energie 


Diese Veränderung der 0 
Weissschen Bezirke wurde 4 
merst von Fryer und 
Fowler mit Hilfe des ma- 200-100 0 oo 200. #0 
gnetooptischen Kerr-Effek- Temperotur 


2 
tes beobachtet ). Zur gleic hen Abb. 1. Temperaturabhängigkeit der Konstanten K, 
Zeit und unabhängig von und K, in Formel (1), nach Bozorth') 


') R. M. Bozorth, Ferromagnetismus, D. van Nostrand Comp. New York 1951. 
*) E.M. Fryer u. Ch. A. Fowler, Physic. Rev. 98, 270 (1955). 
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den genannten Autoren versuchten wir, durch Entwicklung einer temperatur. 


beständigen Suspension den Einsatz der Bitter-Streifentechnik bei die. > 
sem Problem möglich zu 

machen. Es gelang, eine hr 

se 4 


Suspension aus Eisen 
oxyden in Paraffinöl her- 
zustellen, die bis zu Tem- J 
peraturen von etwa 400°C J wisse 
zu verwenden ist. 

Mit Hilfe der tempe- 
raturbeständigen Suspen- 
sion wurden unter anderem 
magnetische Pulvermuster 
auf Kobalt im Temperatur- 
gebiet von Zimmertempe- 
ratur bis 380° C beob- 

achtet und photographiert. 
_ olmm Die Abb. 2 und 3 zeigen 
Abb. 2. Magnetische Pulvermuster bei 20°C auf einem 2- B. die Veränderungen 
Kobaltkristall, dessen hexagonale Achse senkrecht zur der magnetischen Bezirks- 


Probenoberfläche steht. Die Eindrücke des Mikro- struktur auf einem Kobalt- 
härteprüfers dienen zur Markierung 


kristall bei Temperatur- 
erhöhung von 20°C auf 
380°C. Aus dem ‚,Stern- 
chenmuster** der Abb. 2 ist 
zu entnehmen, daß die 
Probenoberfläche parallel 
zur Basisfläche des Kristall- 
gitters liegt?). Die magne- 
tische Vorzugsrichtung 
steht bei Zimmertempe- 
ratur senkrecht zu dieser 
Ebene, woraus sich die 
Bildung des eigenartigen 
Musters erklärt. Die gleiche 
Stelle der Probe enthält 
bei 380° C offensichtlich 
o Tv 
Abb. 8. Dieselbe Stolle der Probenoberfläche wie bei °° — Bloch — vue 
Abb. 2, jedoch entstanden die Muster bei 380° C. und „‚Dolche‘, ein Zeichen 
Probe und Suspension (feinstes Eisenoxydpulver in dafür, daß eine ausgeprägte 
Paraffinöl) wurde zunächst getrennt auf die angegebene Vorzugsrichtung parallel 
Temperatur ‚erhitzt, dann vereinigt und nach einigen zur Oberfläche vorhanden 
Sekunden wieder getrennt. Der an der Probe haftende . 


Suspensionsrest wurde bei der gleichen Temperatur ist. 
eingetrocknet Eine _umfangreichere 


Veröffentlichung, in der 
über die Untersuchungen bei höheren Temperaturen ausführlicher berichtet 
werden soll, ist in Vorbereitung. 


olmm 


“ 
» 
C. Elmo Physic. Rev. 58, 757 (19 
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Herrn Prof. Dr. W. Schiitz bin ich für seine Anregung zu diesen Experi- 
menten zu besonderem Dank verpflichtet. Ebenso habe ich Herrn Dipl.-Phys. 
W. Stephan für die Überlassung des von ihm hergestellten Kobaltkristalles 


sehr zu danken. 


Jena, Institut für magnetische Werkstoffe der Deutschen Akademie der 


Wissenschaften zu Berlin. 
= 
Bei der Redaktion eingegangen am 5. November 1955. wi 
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Bemerkungen zur Arbeit: 
Grundzustand des Wasserstoffmolekils* 


Durch eine igheteit von H. Schulz stellte sich heraus, daß in der 
oben genannten Arbeit!) (im folgenden als I zitiert) bei der Berechnung des 
Normierungsintegrals der Eigenfunktion 


(tpt) =e ? Vl+ePy(m)e *" + Ps (m) (1) 
ein Fehler unterlaufen war. In den Nennern in den Gleichungen (I,11), (1,13), 
(I, S. 74 oben), (I, S. 75 oben) muß 


| 
- 1 1 1 1 4 


ersetzt werden. Die Korrektur gibt natiirlich nur kleine Anderungen fiir die 
Gesamtenergie des H,-Molekiils. Da aber seine Bildungsenergie als Differenz 
zweier fast gleicher Zahlen entsteht (Gesamtenergie minus ein Doppelrydberg), 
ist die Änderung für die Bildungsenergie beträchtlich. Die günstigsten Para- 
meterwerte für die Eigenfunktion (1) liegen in der Nähe von k = 1,7; c = 0,3; 

= 1,40 und führen zu einer Bindungsenergie von 3,59 eV. Infolgedessen ist 
also in (I, Abb. 9) der mit Gl. (10) gekennzeichnete Strich um 7 mm nach 
links und um 14 mm nach oben zu verschieben. Damit ist die in I gezogene 
Folgerung, daß bei gleicher Zahl von variierbaren Parametern in elliptischen 
Koordinaten mit den Kernen als Brennpunkten aufgebaute Eigenfunktionen 
den aus Wasserstoffatomeigenfunktionen mit Abschirmungsparametern auf- 
gebaute Eigenfunktion überlegen sind, falsch. Man muß schon wie James 
und Coolidge mit einer großen Zahl von Parametern arbeiten, um zu guten 
Ergebnissen zu gelangen. 


1 1 1 


Unter diesen Umständen hat es keinen Sinn, die Beziehungen in (I, Abschn. 
III) noch einmal auszuwerten. Die Eigenfunktionen yı, yy, (I, (16)) führen 
nach Berücksichtigung der Korrektur zu Bindungsenergien von 1,95 und 
1,90 eV. Das Kästchen in (I, Abb. 9) ist also um 6mm nach oben zu ver- 
schieben. 

Es ist vielleicht noch interessant, dieses letzte Ergebnis mit der Unter- 
suchung von H. Preuß?) zu vergleichen. In I wurde eine von v. Mohren- 
stein?) angegebene Eigenfunktion des Wasserstoffmoleküls unverändert über- 
nommen und die zugehörige Bindungsenergie mit Hilfe des exakten Hamil- 


') B. Kockel, Ann. Physik (6) 15, 64 (1954) 
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B. Kockel: Der Grundzustand des Wasserstoffmoleküls 


tonoperators bestimmt und ein um etwa den Faktor 2 zu kleines Ergebnis 
gefunden. Preuß übernimmt den Hamiltonoperator v. Mohrensteins 


nur im nur im 
A-Halbraum B-Halbraum 


(1,2 die beiden Elektronen; A, B die beiden Kerne) und bestimmt durch ein 
Variationsverfahren die zugehörige Bindungsenegrie. Er findet das Doppelte 
des experimentellen Wertes (9,2—9,8 eV je nach dem Ansatz zu variierender 
Funktionen) ?). 

4) Verändert man den Hamiltonoperator, so kann man durch ein Variationsver- 
fahren natürlich sowohl Werte der Bindungsenergie finden, die größer, als auch solche, 
die kleiner als der experimentelle Wert sind. 


Leipzig, Theoretisch-Physikalisches Institut der Universität. 
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Über der Thomas-Fermi- 


Von H. A. Buchdahl eu 


Es werden gewisse elementar-algebraische Funktionen als Approximationen 
der Lösung der Thomas-Fermi-Differentialgleichung betrachtet und ta- 
bellarisch dargestellt. 


“= . . . . ~ . 
bee In einer kürzlich erschienenen Arbeit!) hat Tietz die Funktion 


mit c= (1) 


als Näherungslösung der Thomas-Fermi-Differentialgleichung 


mit den Randbedingungen 
y(0) = 1, yo) = 0 (3a, b) 
vorgeschlagen. Die Funktion (1) hat er erhalten, indem er, analog einem = 


fahren von Brinkman?), in (2) auf der rechten Seite V x durch ( (62 y y) 
ersetzt hat. Daraus erfolgt eine abgeänderte Differentialgleichung, welche das 
singuläre Integral (1) besitzt. Die Konstante c ist dabei durch die aus (2) und (3) 


folgende Identität 
00 
J= | Vay"dx=1 #73 

6 


bestimmt. In der Tat ist aber nun Vx y nur in einem kleinen Teil des Gesamt- 
intervalls 0, oo näherungsweise konstant; und überhaupt ist das Brinkman- 
Tietzsche Verfahren nicht recht sinnvoll. Denn es wird angenommen, daß y 
eine gewisse Funktion f(x) von x ist, nur um dann mit Hilfe dieser Annahme 
für y eine von f(x) ganz verschiedene Funktion g(x) zu erhalten. Dem falschen 
asymptotischen Verhalten von f(x) entsprechend ist auch das asymptotische 
Verhalten von g(x) verkehrt. Die Tietzsche Lösung ergibt für x > oo 


(4) 


a y > const = c~. (5) 


Andererseits besitzt die Funktion J, (P) = x? y von Coulson und March’) 
ihrer Kurve II, Abb. 3 zufolge, bei P (= (y/x)':) x 0,03 ein steiles Maximum 
(= 2,5), und fällt sodann für kleinere P, d. h. größere x, rasch nach Null ab. 

Die vorliegende Arbeit hat den Zweck, den eben erwähnten Mangel beiseite 
zu bringen. Ich beschränke mich dabei wesentlich auf möglichst einfache 


5 T. Tietz, Ann. Physik (6), 15, 186 (1955). 

 %) H. Brinkman, Physica XX, 44 (1954). 

Sat *) C. A. Coulson u. N. H. March, Proc. Physic. Soc. (A) 68, 867 WR: 
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elementar-algebraische Näherungsfunktionen, 7(x), welche im Unendlichen 
das asymptotische Verhalten von y, nämlich?) 


y~ 144/23, (x + oo), (6) 


richtig wiedergeben®). Will man nun Tabellen von y nicht verwenden, so 
darf 7 nur eine unbekannte Konstante enthalten, welche man dann mit Hilfe 
von (4) bestimmt. Die für den gegenwärtigen Zweck geeignete einfachste 
Funktion ist wohl 

1 


7 
Be erryeur) mit a,b >0, und 144ab? = 1. (7) 


Sie leistet (3) und (6) Genüge, und hat weiterhin den Vorteil, daß sich J ele- 
mentar als Funktion von a und 6 auswerten läßt. Setzt man, wenn bja > | 
1 
ist, — = = 1+, so erhält man 
12a"! & 


J = 3(& + 1) [-E(5& + 6 & — 1) arctan + (15 + N]. (8) 
Die Gleichung J = 1 hat die Wurzel § = 0,25551. Also 


YO) (1 + 0,02966 x) (1 + 0,4839 x)? 


Nimmt man b/a < 1 an, und geht ganz ähnlich vor, so erhält man die zweite 
Möglichkeit 


1 

(1 + 1,0738. x) (1 + 0,08042 x)’ (10) 
Es ist von Interesse zu bemerken, daß man 7,,, (aber nicht Jz) !) als eine ein- 
fache Modifikation der Tietzschen Funktion ansehen kann, bei welcher der 
letzteren das nur für große x wesentliche Korrektionsglied (1 + 0,02966 x)-1 
zugefügt ist, und demgemäß die Konstante c = 0,5376 durch die angemessenere, 
etwas kleinere, Konstante 0,4839 ersetzt ist. Die Maximalwerte von J,(P) 
sind (bei dreistelliger Genauigkeit) 2,29 bei Verwendung von 7,1), bzw. 2,71 
bei 7. In Tabelle 1 werden die durch (9) bzw. (10) gegebenen Funk- 
tionswerte verglichen einerseits mit denen, die (1) entspringen (y7,), 
andererseits mit bekannten Werten y,, von y. Für 0 < «< 30 sind die letzte- 
ren aus der Tabelle von Bush und Caldwell®) entnommen, für x > 30 aus 
der Sommerfeldschen asymptotischen Formel”) berechnet. Die Rechnungen 
sind durchwegs auf drei Stellen genau. j,, ist offenbar eine bessere Näherung 
als J). Für x< 10 ist sie im ganzen nicht schlechter als y7,, aber für große 
x natürlich sehr viel besser. 

Die einfachste elementare Funktion, die (3) und (6) entspricht, und die 
zwei frei wählbare Konstanten enthält, ist wohl 

mit 144y = 1. (11) 

4) A. Sommerfeld, Z. Physik 78, 283 (1932). ae : 

5) Nur bei der letzten der unten betrachteten Approximationen ist diese Bedingung 
nicht erfüllt, aber nur in dem Sinn, daß der Koeffizient 144 durch eine etwas kleinere 
Zahl ersetzt ist. 

6) V. Bush u. S. H. Caldwell, Physic. Rev. 38, 1898 (1931). 
7) 8.2, GI. (3). 
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> Tabelle 1 


x Ynu Yri Ya) | 4 (2) Ya) 
0 1 1 1 1 1 
0,5 0,607 0,621 0,639 0,601 0,618 
1,0 0,425 0,423 0,441 0,413 0,425 
2,5 0,193 0,182 0,191 0,188 0,189 
5,0 0,0788 0,0735 0,0745 0,0868 0,0773 
7,5 0,0408 0,0395 0,0382 0,0430 0,0406 
Be 0,0244 0,0246 0,0226 0,0269 0,0244 
ae 0,0158 0,0168 0,0147 0,0171 0,0160 
0,0109 0,0122 0,0101 0,0122 0,0111 
<2 20,0 0,0058 0,00724 0,00551 0,00654 0,00602 
30,0 0,0022 0,00341 0,00220 0,00259 0,00239 
40,0 0,00108 0,00197 0,00110 0,00128 0,00122 
100.0 0.0,987 0.0,339 0.0,103 0.0,113 0,05108 


a und b bestimme ich jetzt aus den bekannten Werten von y (1) und y (10), 
und erhalte auf diese Weise 
J, = (1 + 1,126 x + 0,2178 x? + 0,006944 a?)-1 
1 
(1 + 0,03978 x) (1 + 0,1962 x) (1 + 0,8900 x) ' 


Die sechste Spalte der Tabelle 1 gibt aus (12) berechnete Funktionswerte. 
Diese Approximation ist also bedeutend besser als die Tietzsche. Sie gibt 
für /,(P) den befriedigenden Maxi- 
Tabelle 2 malwert 2,48, und zwar bei x = 13,4, 

d.h. bei P = 0,032. 
Fa), 92), nähern sich, wie von 
Anfang an verlangt wurde, der 


| Ya Ynu 


1 1 1 
0 0,771 0,755 wahren Funktion y fiir hinreichend 
q 0,50 0,617 0,607 große x beliebig nahe an. Man sieht 
— ae 0,507 0,502 nun aber, daß selbst für x = 100 7 
100 0,425 0,425 immer noch um etwa 10% zu groß 
1,75 0,273 0,276 ist. Fü 
2,50 0,190 0,193 ist. Für praktische Zwecke könn 
“S75 0,117 0,118 man also versuchen in (11) die Be- 
5,00 = oe dingung 144y = 1 fallen zu lassen, 
150 0.0408 0.0408 und a, ß, y aus drei bekannten 
I 8.75 0.0311 0.0310 Werten von y zu bestimmen. Man 
un‘ 0,0244 0,0244 kann erwarten, auf diese Weise trotz 
a 12,5 0,0158 0,0158 der Einfachkeit der Approximations- 
= funktion eine sehr gute Näherung 
30.0 0.00228 0.0022 zu erhalten, wenn man von allzu 
40,0 0,00111 0,00108 großen x absieht. Zieht man etwa 
100 0.0,968 0,0,987 y(1), y(7,5) und y(15) hinzu, so be- 
200 0,0,136 0,0,144 kommt man 
1 
7 = — 13 
Wa) (T + 0,05727 x) (1 + 0,1536 x) (1 + 0,9288 x) 4 3 (19) 
Fiir sehr groBe x hat man also “u 
ya) ~ 122,4 a8, (x > 00), ag = 
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‚daß man für hinreichend große x einen von % unabhängigen relativen 

4 fehler von etwa 15% begehen wird; in der Tat ist er selbst bei x = 100 nur 

— ungefähr 1%- Tabelle 2 zeigt, mit welcher G4) 
jie wahren Funktionswerte approximiert. 
Hobart (Australia), Physies Department, University of 
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Zur Theorie der in Metallen 


Von Karl-Heinz Schramm eve 
Mit 1 Abbi 4 
Mit 1 ildung 4 


Inhaltsübersicht 


Es wird zunächst eine Beziehung hergeleitet zwischen der Fortpflanzung- 
geschwindigkeit kleiner Störungen im System der Leitungselektronen eines 
Metalles einerseits und derjenigen im Ionensystem andererseits. Daraus er- 
gibt sich dann eine Formel für die Schallgeschwindigkeit in Metallen, in die 
die folgenden Größen eingehen: Planck sches Wirkungsquantum, Elektronen- 
masse, Loschmidtsche Zahl, Atomgewicht, Dichte, Anzahl der Leitungs 
elektronen pro Atom und Poissonsche Zahl der Querkontraktion; der 
Elastizitätsmodul kommt in der erhaltenden Formel nicht vor. Zum Schluß 
werden die theoretischen Ergebnisse mit einigen experimentellen Meßwerten 
verglichen. 


1. Einleitung 


_ Ein Metall besteht strukturell aus 2 Teilchensystemen: dem System der 
Metallionen und dem der Leitungselektronen. Als Analogon zur Ausbreitung 
elastischer Wellen im lonensystem (Schall!) soll die Fortpflanzung kleiner 
Störungen im Elektronensystem diskutiert werden. Dabei interessiert ins- 
besondere die Fortpflanzungsgeschwindigkeit solcher Störungen (,,elektronische 
Schallgeschwindigkeit‘‘ v) entsprechend der „klassischen Schallgeschwindig- 
keit‘‘ c im Ionensystem. 
Die Notwendigkeit zur E inführung einer solchen Geschwindigkeit v ergab 
“ sich bei der Ausarbeitung meiner statistischen Theorie der Supraleitung)). 
Ich ging dort von der — erstmalig von de Lar Kronig geäußerten — Vor- 
stellung aus, daß sich das System der Leitungselektronen in einem Supn- 
leiter zumindest näherungsweise wie ein Elektronengitter verhalte. Der 
Phasenübergang normalleitend—supraleitend entspricht in dieser Vorstellung 
etwa einer Kondensation des Elektronengases im Metall. Zur Berechnung 
der spezifischen Wärme eines solchen Elektronengitters benutzte ich: — it 
Anlehnung an die Vorgehensweise beim klassischen Kristall — die Debye 
sche Eigenschwingungsmethode unter angemessener Berücksichtigung der 
Tatsache, daß, durch die Anwendung der Fermi-Statistik bedingt, das 
schwingende Elektron nach einem ganz anderen Strahlungsgesetz strahlt ak 
nach dem Planckschen des klassischen (Boltzmann) Oszillators. An die 
Stelle der Schallgeschwindigkeit c des klassischen Kristalls tritt dabei zudem 
die erwähnte elektronische Schallgeschwindigkeit v. Diese wurde aber nieht 


4) K. -H. Schramm, Ann. Physik [6] ms 55 (1956). 
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wellenmechanisch berechnet, sondern nur auf die empirisch bestimmte Sprung- 
temperatur 7, des Supraleiters zurückgeführt. Er ergab sich: 
k 4 Ms 
k = Boltzmann-Konstante; h = Plancksches Wirkungsquantum; 
L=Loschmidtsche Zahl; o = Dichte des Metalles; 
A = Atomgewicht; f = Anzahl der Leitungselektronen pro Atom. 


(Die Gl. (1) wurde in der genannten Arbeit nicht explizit hingeschrieben. 
Man erhält sie durch Gleichsetzen der beiden dortigen Gln. (18b) und (26b/e), 
wenn man die Konstante A in (18b) aus (19) bestimmt, indem man dort für 
J (7) die nullte Näherung J (co) einsetzt.) 

Numerisch ergibt sich nach Gl. (1) z.B.: | oe N 


; 


Z Symb. A f o [g- em] | T, v [cm - | c 
13 Al 26,97 8 | 2,70 | 1,14 0,15 - 105 5,1 - 105 
30 | Zn 65,38 2 | 7,13 0,79 | 0,12- 105 3,7 - 105 
50 Sn 118,7 4 7,28 3,72 | 0,53 - 105 2,6 - 105 
80 Hg 200,6 2 | 13,595 4,17 0,73 - 105 2,7.105 
82 Po. | 207,2 a | 11,34 7,26 | 1,08 - 105 1,7 - 105 


Man erkennt, daß im allgemeinen v < c ist und nur bei den Elementen mit 
den höchsten Sprungtemperaturen » ~ c wird. Da diese Theorie der Supra- 
leitung bezüglich der berechneten Werte für die spezifische Wärme und die 
Schwellwertkurve H, = H,(T) des kritischen Magnetfeldes relativ gut mit 
der Erfahrung übereinstimmt, so wird man den berechneten Zahlenwerten für 
v auch eine gewisse Realität nicht absprechen können. Die Ausbreitungs- 
geschwindigkeit v „elastischer Wellen‘ im Elektronensystem ist also wesent- 
lich kleiner als die Lichtgeschwindigkeit mit der sich ihr elektromagnetisches 
Feld fortpflanzt. 

Im folgenden wird nun zunächst eine Formel für diese elektronische Schall- 
geschwindigkeit v im Supraleiter hergeleitet. Diese hat die Gestalt: 


v2 = 1840: A - (x? — c?) (A = Atomgewicht). (2a) 


Zur numerischen Berechnung von v ist Gl. (2a) jedoch nicht geeignet, da wegen 
v< .c, in der Klammer die sehr kleine Differenz zweier sehr großer (und nicht 
genau genug bestimmbarer!) Größen steht. Damit ergibt sich aber die Möglich- 
keit, die Fragestellung umzukehren und Gl. (2a) zur Berechnung von c zu be- 
nutzen, da wegen v<c mit genügender Genauigkeit: 
(2b) 
gilt. 
Der so berechnete Wert c der Schallgeschwindigkeit bezieht sich herlei- 
tungsgemäß zunächst nur auf den Supraleiter. Experimente?) haben jedoch 
ausgewiesen, daß die Schallgeschwindigkeit in einem Metall im ei ain 


) Squire, Physic. Rev. 95, 1126 (1954) 
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Zustand von derjenigen bei Zimmertemperatur nicht merklich versclitii 
ist. Gl. (2b) gilt also auch bei Zimmertemperaturen, d. h. im normalleitender 


2. Die Strömungsgleichung der Elektronen 


Wie sich thermodynamisch zeigen läßt!)3), bleiben sämtliche Zustand. 
funktionen des Elektronensystems beim Phasenübergang normalleitend- 
supraleitend eines Metalles (in Abwesenheit äußerer Magnetfelder!) unver. 
ändert. Dies bedeutet, daß auch für das Elektronengitter ein ,,hydrostati- 
scher Druck“ p existiert, der sich, wie für das Elektronengas, nach Som. 
merfeld*) zu: 


h = Plancksches Wirkungsquantum; m = Elektronenmasse; 
N = Anzahl der Leitungselektronen im Volumen V. 


(Genau genommen gibt Gl. (3a) nur den Nullpunktsdruck p, für T = 0° K; 
die temperaturabhängigen Terme sind jedoch um mehrere Zehnerpotenzen 
kleiner, so daß man sie hier vernachlässigen kann.) 


Bezeichnet man mit D die Massendichte des Elektronensystems, so formt 


sich Gl. (3a) um zu: ER 


Wegen der Existenz eines solchen allseitigen Druckes p im Elektronengitter 
kann man dieses „strömungsmechanisch‘ als „Flüssigkeit‘‘ schematisieren 


und die hydrodynamischen Grundgleichungen anwenden: re 

av 1 
a 


Dabei bedeutet: 


& die auf das Elektronensystem wirkende Volumenkraft, bezogen auf 
die Masseneinheit, 
b die Strömungsgeschwindigkeit der Elektronen. 


Durch Einführung des Ausdruckes (3b) für p erhält man aus Gl. (4a), 
wenn man noch das quadratische Glied (v, grad ») vernachlässigt: 


(=n) grad (D*/). (4e) 


Es soll nun die Fortpflanzung kleiner Dichteschwankungen o untersucht 
werden: 


D=D, (1+0; o<l. (5a) 


®) v. Laue: „Theorie der Supraleitung‘“, Springer Verlag, Berlin-Göttingen-Heidel- 
berg, 2. Aufl. 1949. 

*) Vgl. z.B.: A. Sommerfeld: „Thermodynamik und Statistik“, Dieterichsche 
Verlagsbuchhandlung, Wiesbaden 1949. 
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‘Wegen: div (Dv) =D-divv + (v, grad D) ~ D,-divo erhält man aus 
Gl. (4b) in Verbindung mit Gl. (5a): 

do 
Durch Bildung der Divergenz von Gl. (4c) und Einsetzen der Gln. (5) erhält 
man dann die Strömungsgleichung: 
40. (6a) 


div 6 = — m 

Soweit hält sich der Gang der Rechnung eng an den, den man zur Be- 
reehnung der Schallgeschwindigkeit in Flüssigkeiten normalerweise ein- 
schlägt. Während man dort aber für die Volumenkraft © = 0 einsetzen, 
d.h. die Flüssigkeit als kräftefrei betrachten konnte, ist das hier nicht möglich, 
da die Leitungselektronen stets unter dem Einfluß der elektrostatischen An- 
ziehungskräfte der Metailionen stehen. 


3. Die Berechnung der Volumenkraft © 


Das Leitungselektronensystem in einem Metall steht nach Gl. (3a) unter 
der Wirkung eines expandierenden Druckes p, der für Kupfer z. B. den Wert 
p= 3,8- 105 at annimmt. Dieser enorme Druck entspricht der elektrostati- 
schen Anziehung zwischen Ionen und Elektronen; die Ionen halten das Elek- 
tronensystem gewissermaßen zusammen. Steigt nun die Elektronendichte D 
und damit der Druck p entsprechend GI. (5a) über den normalen Wert an, 
so wird dieses Kräftegleichgewicht gestört und es resultiert eine expandierende 
bzw. komprimierende Volumenkraft & auf das Elektronensystem 

Zur Berechnung von ® kann folgende Überlegung dienen: Man denke 
sich die Dichte 9 des Ionensystems durch Einwirkung geeigneter äußerer 
Kräfte im gleichen Sinne verändert wie D: 


D = Dy: (1+a)70 =@: (1+ 0). 


Das Ionensystem soll also teils gestaucht, teils gedehnt vorgestellt werden.. 


Hierdurch wird das durch die Änderung von D nach Gl. (5a) gestérte Gleich- 
gewicht der auf das Elektronensystem wirkenden Kräfte wieder hergestellt. 
Datür wirkt auf das Ionensystem jetzt eine resultierende Volumenkraft $, 
die dieses in der geforderten Weise gespannt hält. Da sich § nach den Ge- 
setzen der Elastizitätstheorie berechnen läßt, ergibt sich damit eine Möglich- 
keit, & auf die elastischen Daten des Metalles zurückzuführen. 
Zur mathematischen Durchführung d’eses Gedankens soll zunächst 
Gl. (6a) auf die Strömung durch einen unendlich langen, kreiszylindrischen 
Draht spezialisiert werden. Alle auftretenden Größen hängen dann in 1. Nä- 
herung nur von einer Ortskoordinate, der Richtung x der Drahtachse ab: 
0G, @o, 1 3 ls ney 
(Vernachlässigt worden ist dabei die Strömung der Raumladungen radial 
auf die Leiteroberfläche zu.) 
Um die Kraft G(x) als Funktion der Elektronendichte D(x) zu bestim- 
men, denke man sich den Draht in lauter kleine Elemente AV der Länge 4x 
aufgeteilt, wobei Az so klein gewählt werden soll, daß sich D(x) innerhalb Ax 


(6b) 
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satiate nur ew andert. Sodann greife man 2 solcher Elemente 4 ail B 
zur näheren Betrachtung heraus: 


Nach Voraussetzung beträgt die mittlere Elektronendichte in B: 
Dz =—:[D(x) + D(a + Az)], 
(7a) 
und analog in A: 


Ist nun D, die Elektronendichte im Gleichgewichtszustand und D, >D, 
so wirkt auf die Elektronen in B = AV die mittlere expandierende Kraft 
D,:G,:AV von B nach A. Da G, in 1. Näherung nur vom Dichtegra- 
dienten @D/dx abhängt, so kann man ohne Beschränkung der Allgemeinheit 
annehmen, daß man es bei den beiden Drahtelementen A und B gerade mit 
denjenigen zu tun hat, deren mittlere Gesamtdichte D = 1/, (Da + Dy) =D, 
ist, d. h. es sollen in B gerade so viele Elektronen zu viel sein, wie in A zu 


wenig sind. Mit D(x) = D,- [1 + o(z)] gilt also: Be u 
Pr 


Entsprechend dem eingangs nie wird der Kraft D,-G,-AV in 
1. Näherung gerade durch die Kraft X, - AV das Gleichgewicht gehalten, die 
man aufwenden muß, um das Volumen von B um öAV zu stauchen und gleich- 
zeitig das von A um 6AV zu dehnen, und zwar so, daß in B jetzt statt der 


Ionendichte 9, die mittlere Dichte 

Op = 2%" [2 (?e) 

K, berechnet sich nach dem Heskeniien Gesetz zu: 
or 1 AV 4 


E = Elastizitätskonstante; «4 = Poissonsche Zahl der Querkontraktion. 
(Da die Kraft K,-AV gleichzeitig A dehnen und B stauchen soll, so entfällt 
die Hälfte davon !/,- K,- AV auf die reine Stauchung von B im Sinne der 
Anwendung des Hookeschen Gesetzes!). 

In den Gln. (8a) und (8b) bedeutet 64x die bei der Stauchung von B um 
öAV auftretende Verkleinerung von Ax und dr die gleichzeitig eintretende 


Vergr 
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Vergrößerung von r: 


dAV (AV +5AV)— AV [(r — Or)? - (Ax + — Ar] 


iv 7 AV Ag 
d.h. 
AV Az r 
Aus den Gln. (8) erhält man für K,: 
Ax 1— (9a) 
Nach Voraussetzung besteht zwischen o und AV entsprechend Gl. (7c) die 
Beziehung: 
AV = on (AV — =o(1 (av — 84V), 
(9e) 


Durch Einsetzen der Gin. (7) und (9c) erhält man aus Gl. (6b) für die 
elastische Störung o im Elektronensystem die Wellengleichung: 
Coa fi h? ( 3 E 1 | 
\3 m? \8-a-m D, 1-24‘ 
Hieraus liest man unmittelbar die elektronische Schallgeschwindigkeit v ab: 

1 M 3 Vis E 1 
Di: 


Die Massendichte D, des Elektronensytems läßt sich leicht auf die ,,normale 


(10a) 


(10b) 


Dichte o des Metalles zurückführen: — = — =? f ZEN, 
7 m V A 
L- f 


Beachtet man nun noch, daß] += =c die Schallgeschwindigkeit (longitudinaler 


Wellen in unendlich langen Drähten) bedeutet, so erhält man aus (10b) und 
(10e) die Gleichung: 
Diese Gl. (11a) entspricht der Gl (2a). Man erhält aus ihr für die Schall- 
geschwindigkeit c entsprechend Gl. (2b): EHE 


(h = Wirkungsquantum; m= Elektronenmasse; LL = Loschmidt-Zahl; 
A = Atomgewicht; o = Dichte; f = Anzahl der Leitungselektronen pro Atom; 
#=Poissonsche Zahl der Querkontraktion) 


i 
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bzw. numerisch: 


© = 4,0905 - 108 — 2 wem - see. 


Vergleiche zwischen Theorie und Experiment 


Leider läßt sich die entwickelte Theorie nur in wenigen Fällen mit der Er- 
 fahrung vergleichen. Insgesamt wurden 9 Beispiele durchgerechnet. Für alle 
anderen Metalle liegen entweder keine Meßwerte über die Schallgeschwindig- 
keit c bzw. die Poissonsche Zahl u vor, oder aber die Atomtheorie liefert 
keine brauchbaren Aussagen über die Anzahl f der Leitungselektronen pro 


Atom. 
RN Die Tabelle 2 bringt die Angaben zu einigen Elementen, die nachweislich 
auch wirklich supraleitend werden, auf die also die vorstehenden Überlegungen 
anwendbar sind. Man erkennt, daß die berechneten Werte relativ gut mit den Es w 
gemessenen übereinstimmen. EMK be 
Anders ist dies für die Elemente Cu, Ag und Au, die alle 3 der ersten Gruppe # jes Lich 
des periodischen Systems angehören und alle 3 daher nicht supraleitend werden. # ywuerstı 
dieser V 
Tabelle 25) Reaktio 
c [em - an PbS: 
| Symb. | elg- f berechnet | gemessen gradien' 
| Photo-F 
13 Al 26,97 2,72 3 0,34 5,18-10° | 5,1 - 106 von Gi 
30 Zn 65,38 7,12 2 0,2/0,3 3,30/2,72-105, 3,7- 105 
48 ‘d 112,4 8,64 2 0,30 1,85-105 | 2,3. 108 
50 | Sn 118,7 7,28 4 0,33 2,74-105 2,6 - 105 
82 Pb 207,2 11,34 4 0,45 1,08-10® 1,3. 105 
8 | Bi 209,0 9,80 5 0,33 2,28-105 1,8- 10° 
Tabelle 35) Neb 
ohne d 
Symb. 4 elg- em] f | berechnet | gemessen auftrete 
Görlie 
29 Cu 63,54 8,93 1 0,35 | 1,46. 105 3,9 - 105 Liehtso 
47 | Ag 107,9 10,51 1 0,38 | 0,89-105 | 2,7.10 allen 
79 Au 197,2 19,3 1 0,42 | 0,54 - 105 | 2,1. 108 . 
perime? 
erstma‘ 
u Alles in allem läßt sich wohl feststellen, daß eine quantitative Berechnung end: 
der Schallgeschwindigkeit c nach Gl. (11b) nicht möglich ist, daß Gl. (11b) ment 
aber sicherlich einen interessanten Fingerzeig fiir eine kiinftige atomistische Ladun; 
_ Determination der elastischen Konstanten vines Metalles abgibt. interkr 
5) Die Meßwerte wurden entnommen aus: F. Kohlrausch: ‚Praktische Physik“, tialwal 
Bd. 2, Teubner, Leipzig 1950. Sperrse 
der Sc 
Mülheim a.d. Ruhr, z. Z. Deutsche Versuchsanstalt für Luftfahrt. und el 


Bei der Redaktion ei 18. A 
i der Redaktion eingegangen am 18. August 1955. 
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Mit 10. Abbildungen 
en 
len Es wird über eine an polykristallinen PbS-Schichten auftretende Photo- 


EMK berichtet, welche nach Größe und Vorzeichen von der Einfallsrichtung 

pe Ff des Lichtes abhängt. Als Entstehungsursache dieser Photo-EMK werden von 

en. # Sauerstoffadsorption herrührende Oberflächenzustände angesehen. Unter 
dieser Voraussetzung werden Rückschlüsse auf den Schichtenaufbau und den 
Reaktionsmechanismus bei Belichtung gezogen. Ferner werden Untersuchungen 

in PbS-Schichten mitgeteilt, welche nach Formierung in einem Temperatur- 

n gradienten eine beträchtliche Unsymmetrie in der Winkelverteilung der 

—  Photo-EMK zeigen. Dieses Verhalten kann auf Grund einer Thermodiffusion 

von Gitterstörstellen wenden. 


1. Einleitung 


Neben der Widerstandsänderung bei Belichtung können an PbS-Schichten 
— AH ohne die Mitwirkung von Metallkontakten photoelektromotorische Kräfte 
auftreten. So sind von Sosnowski, Soole und Starkiewiez?) und von 
— #6örlich und Krohs?) beim Abtasten solcher Schichten mit einer feinen 
liehtsonde örtlich verschieden kleine Photospannungen wechselnden Vor- 
zeichens festgestellt worden. In dieser Tatsache hat man des öfteren eine ex- 
perimentelle Bestätigung der von Sosnowski, Starkiewicz und Simpson‘) 
erstmalig vertretenen Auffassung gesehen, wonach der hohe Dunkelwider- 
ng [stand und die Photoleitfähigkeit von PbS-Schichten nicht lediglich auf die 
im unbelichteten bzw. belichteten Zustande vorhandene Konzentration freier 
Ladungsträger zurückzuführen ist, sondern in erster Linie auf die Existenz 
nterkristalliner Sperrschichten, welche bei Belichtung die Höhe ihres Poten- 
tialwalles vermindern. Nach den genannten Autoren soll es sich bei diesen 
Sperrschichten um ,,p-n-junctions‘‘ handeln. Dies würde für die Struktur 

der Schichten bedeuten, daß etwa gleichzeitig vorhandene elektronegative 

und elektropositive Störstellen sich nicht völlig statistisch verteilen, sondern 


- 


1) Gekürzte Wiedergabe der Dissertation, G. Schwabe, Jena 1966. SEN, 
*) L. Sosnowski, B. W. Soole, J. Starkiewicz, Nature 160, 471 (1947). Scag i 
3) P. Görlich, A. Krohs, Jenaer Jahrbuch 1952, S. 155. ag 


‘) L. Sosnowski, J. Starkiewicz, O. Simson, Nature 159, 818 (1947). aie 
Ann. Physik. 6. Folge, Bd. 17 „22 17a da 
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daß es kleine benachbarte Bereiche gibt, in denen die beiden Störstellen- 
arten einzeln bei weitem überwiegen. Der mit solchen p-n-Ubergiingen yer. 
bundene Potentialwall kommt bekanntlich durch die Forderung zustand. 
daß die Fermigrenze auf beiden Seiten das gleiche Energieniveau einnehmen 
muß. Die n-Seite trägt dabei gegenüber der p-Seite ein positives Potential 
Bei Belichtung eines p-n-Überganges ändert sich auf jeder Seite die Konze- 
tration der in der Minderheit vorhandenen Ladungsträger relativ sehr viel 
stärker als die Konzentration der in der Mehrheit vorhandenen. Dies bedeutet 
aber, daß zur Einstellung des neuen, durch die Belichtung gegebenen Gleich 
gewichtszustandes Elektronen von der p-Seite auf die n-Seite übergehen und 
entsprechend Defektelektronen von der n-Seite auf die p-Seite. Die n-Seite 
nimmt also bei Belichtung gegenüber der p-Seite ein negatives Potential an. 
In der Natur dieser ,,p-n-junctions“ liegt es daher begründet, daß die an ihnen 
auftretenden Photospannungen orientierungsgebunden sind. Man wird deshalb 
bei vollständiger Ausleuchtung von Halbleiterschichten, welche aus mikro- 
kristallinen ,,p-n-junctions‘‘ aufgebaut sind, im allgemeinen keine Photo- 
spannung erwarten können, da sich die an den einzelnen Elementen entste- 
henden Photospannungen im Mittel kompensieren; denn man kann annehmen, 
daß sich diese p-n-Elemente hinsichtlich ihrer Orientierung statistisch über 
die Schicht verteilen. Starkiewize, Sosnowski und Simpson?) konnten 
aber zeigen, daß bei Belichtung der ganzen Halbleiterfläche dann eine resul- 
tierende Photo-EMK möglich ist, wenn die Schichten zuvor während der 
Formierung bei Temperaturen von 2—300° C mit einem Strome von einigen 
mA belastet wurden. Diese Tatsache kann man zwar so verstehen, daß wegen 
der Beweglichkeit der Störstellen bei höheren Temperaturen durch die Strom- 
formierung ausgerichtete ,,p-n-junctions‘‘ zustande kommen; man darf 
daraus jedoch nicht den Schluß ziehen, daß PbS-Schichten, welche ohne 
Stromdurchgang unter normalen Bedingungen hergestellt wurden, schon , ,p-n- 
junctions“ in einer zur Erklärung der Photoleitfähigkeit ausreichenden Zahl 
enthalten. Jedenfalls ist es denkbar, daß während der Stromformierung ein- 
heitlich ausgerichtete „p-n-junetions‘‘ sogar leichter aus einer völlig ungeord- 
neten Störstellenverteilung entstehen können. 


Die folgenden Untersuchungen werden darüber hinaus zeigen, daß sowohl 
bei partieller als auch bei vollständiger Ausleuchtung der Halbleiterfläche 
an PbS-Schichten photoelektromotorische Kräfte auftreten, welche nicht in 
interkristallinen Sperrschichten ihre Ursache haben können. Auf Grund dessen 
dürfen aber auch Photospannungen, welche bei örtlich begrenzter Belichtung 
entstehen, nicht ohne weiteres als Beweis für die Existenz von interkristallinen 
Sperrschichten angesehen werden, gleichgültig ob es sich dabei um „pr 
junctions“ oder um die von Schwartz®), Gibson’) und Smith*) diskutierten 
Schottky-Randschichten handelt, welche durch Sauerstoffadsorption an 
den Korngrenzen zustande kommen sollen. Trotzdem ist dieses Sperrschichten- 
modell besser mit den folgenden Versuchsergebnissen in Einklang zu bringen 
als die Annahme von statistisch verteilten ,,p-n-junctions‘ (siehe Abschnitt §). 


5) J. Starkiewicz, L. Sosnowski, O. Simpson, Nature 158, 28 (1946). aS” 
6) E. Schwartz, Proc. physic. Soc. (A) 62, 530 (1949). 


7) A, F. Gibson, Proc. physic. Soc. (B) 64, 603 (1951). 
8) R. A. Smith, Semicond. Materials Butterworth, scientif. public. London 1951. 
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G. Schwabe: Von der Lichtrichtung abhängige Photo-EMK 


2. Photo-EMK bei örtlich begrenzter Belichtung 


Ausgangspunkt für die folgenden Untersuchungen waren die Mitteilungen 
von Görlich und Krohs?). Glasplatten der Größe 10 x 20 mm? wurden 2 
nach bekanntem Verfahren?) auf chemischem Wege mit einer 1 bis 2 4 dicken 
Bleisulfidschicht verspiegelt und bei Temperaturen von 530—620° C in einem 
elektrischen Ofen unter Luftatmosphäre formiert. Dabei reagierte das Blei- 
sulfid mit dem Sauerstoff der Luft und die Schichten zeigten bei einem Dunkel- : 
widerstand von 0,5 bis 2 MQ eine Widerstandsverminderung bei Belichtung. 
Die Schichten wurden an den beiden kiirzeren Seiten mit aufgedampften 
Gold- oder Aluminiumelektroden versehen, so daß eine wirksame Halbleiter- 
fläche von 10 x 20 mm? erhalten blieb. 


Diese PbS-Zellen wurden an den Eingang eines gegengekoppelten Gleich- 
sttomverstärkers angeschlossen. Die Anordnung hatte in Verbindung mit 
einem Lichtmarkengalvanometer eine Empfindlichkeit von 1,2 x 10% Volt/ 7 
SkT. Beim Abtasten der Zellen mit einer Lichtsonde (d = 0,05—0,1 mm) a 
wurden ebenso wie von Görlich und Krohs je nach Lage des Lichtpunktes 7 
kleine Photospannungen (i10-4—10-% Volt) unterschiedlichen Vorzeichens 
beobachtet. Jedoch traten diese Photospannungen nur an sehr vereinzelten 
Stellen der Schicht auf. Dazwischen lagen weite Bereiche, welche bei Be- 7 
liehtung zwar eine Widerstandsverminderung zeigten, aber mit der benutzten 
Meßanordnung keine selbständigen Photospannungen erkennen ließen. Diese 
Beobachtungen deuten nicht mit sehr großer Sicherheit darauf hin, daß die j 
PbS-Schichten aus interkristallinen Sperrschichten aufgebaut sind, wobei = 
noch sichergestellt werden miiBte, ob die beobachteten Photospannungen 
wirklich in interkristallinen Sperrschichten ihre Ursache haben. Die folgenden ; 


Versuchsergebnisse lassen 

aber die Möglichkeit offen, a Au 

daß dieser bei punktförmi- wes 

N 

ger Belichtung entstehende > 

Photoeffekt schon durch 

die Beschaffenheit der 

Schichtoberfläche bedingt 

sein kann. 

Es wurde nämlich fol- 

gendes festgestellt Wäh- 

rend bei senkrechtem Licht- aces 

einfall nur an sehr wenigen 1. zur 

Stellen und bei den einzel- E Einfallsrichtung 

nung zum isotropen Verhaltens der 

weniger zahlreic eine S-Schicht 

Photospannungen auf- 

traten, war es möglich, an allen PbS-Schichten dann photoelektromotorische 

Kräfte bei Belichtung jeder beliebigen Stelle zu erzeugen, wenn man den 

Strahlengang des Lichtes zur Flächennormalen neigte. Diese Photospannung 

nahm mit wachsendem Einfallswinkel zu und zeigte bei einem Lichteinfall 


°») P.Gérlich, Die lichtelektrischen Zellen, ihre Herstellung und Eigenschaften, 
Moskau-Leningrad, 1948, S. 119; sowie H. Pick, Z. Physik 126, 112 (1949). 
!%) G. Schwabe, Z. Naturforschg. 10a, 78 (1955). 
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messer des Lichtpunktes vergrößerte ; lediglich 
spannung nahm unter sonst gleichen Bedingu 


erreicht werden. 
Die PbS-Schichten zeigten hinsichtlich 


an der PbS-Schicht, welche gegeben ist durch 


der Lichtrichtung abhängig. 


Zur Messung der Abhängigkeit der Photo- 


Elektroden 


Abdeckvor - 
richtung 


N 
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Abb. 2. Anordnung zur Messung der Winkelabhängig- 
keit der Photo-EMK. Die PbS-Schichten wurden 
dabei auf Calit-Würfel aufgebracht 


ae 252 ee  Ammalen der Physik. 6. Folge. Band 17. 1956 


unter entgegengesetztem Winkel zur Flächennormalen das umgekehrte Vor. 
zeichen. Daran änderte sich nichts Grundsätzliches, wenn man den Durch. 


die absolute Größe der Photo. 
ngen (Einfallswinkel und Be- 


leuchtungsstärke) zu. Bei vollständiger Ausleuchtung der gesamten Halb- 
leiterfläche konnten auf diese Weise Photospannungen von 0,05 bis 0,1 Volt 


dieser Photo-EMK _ keinerlei 


Anisotropie. Die Photo-EMK entstand vielmehr stets in derjenigen Richtu 


den Schnitt der Einfallsebene 


des Lichtes mit der Schichtoberfläche und zwar so, daß die der Lichtquelle 
abgewandte Seite negatives Vorzeichen annahm!!). Um diesen Tatbestand zu 
erhärten, wurden PbS-Schichten, wie aus Abb.1b ersichtlich, mit kreis- 
förmig angeordneten Goldsektoren bedampft, wovon jeweils zwei sich gegen- 
überliegende als Elektroden benutzt wurden. In allen Fällen war die Photo- 
EMK nach Größe und Vorzeichen ausschließlich in der dargelegten Weise von 


3. Abhängigkeit der Photo-EMK vom Einfallswinkel des Lichtes 


EMK vom Einfallswinkel des 


Lichtes wurden die Zellen auf ein Goniometer mit einer Ablesegenauigkeit von 


0,5 Bogenminuten mon- 
tiert. Die Messung wurde 
bei stets vollständiger Aus- 
leuchtung der PbS-Schicht 
im parallelen Lichte vorge- 
nommen. Um den Einfluß 
von Metall-Halbleiterkon- 
takten auszuschließen, be- 
währte es sich, die PbS- 
Schichten auf Calit-Würfel 
aufzuspiegeln und die seit- 
lich angebrachten Elek- 
troden durch zwei zugleich 
als Halterung wirkende 
Backen vor Lichteinstrah- 
lung zu schützen. (Abb. 2). 
Die Abb. 3 zeigt als Ergeb- 
nis einer solchen Messung 
einen Kurvenverlauf, wie 
er für alle untersuchten 


4) In wenigen Ausnahme- 
fällen wurde bei gleicher Ein- 
strahlungsrichtung auch das 
umgekehrte Vorzeichen der 
Photo-EMK beobachtet; je 
doch war dann deren Größe 
stets weitaus geringer (10? bis 
10-2 Volt). 
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G. Schwabe: Von der Lichtrichtung abhängige Photo-EMK a 

Zellen als typisch gelten kann. Wenn wir zur Beschreibung der Winkel- 
abhängigkeit im einfachsten Falle annehmen, daß die Vorzeichenumkehr an 
eine in Richtung der Elektroden wirkende Komponente des einfallenden 
Liehtes gebunden ist, die 
sieh mit sin e ändert, so 
sollte man, da andererseits 


< 


T 


die Beleuchtungsstärke mit +} 10 vA 
cos e abnimmt, einen durch 
sine cose darstellbaren | P 


Kurvenverlauf erwarten. 
Wir sehen, daß dies nur 
mit gewissen Einschrän- 
kungen zutrifft. Und zwar 
liegt das Maximum der 
Photo-EMK nicht bei 45°, \ 


& -Photo-EMK [10V] - 


Ss 


80 


finfallswinkel €° 


sondern nach größeren u 
Winkeln verschoben, bei vg 10 u 
den einzelnen Proben ver- 
schieden zwischen 50 und u 
60°. Dies könnte darauf A»b.3. Abhängigkeit der Photo-EMK vom Einfalls- . 
zurückzuführen sein, daß winkel des Lichtes = 
bei größeren Einfallswin- a 
keln die Vorwärtsstreuung des Lichtes an der Schichtoberfläche gegenüber ande- u 
ren Streurichtungen bevorzugt ist. Andererseits zeigt die Abb. 3, daß die beiden a 


Maxima von ungleicher Größe sind und daß die Vorzeichenumkehr nicht bei 
genau senkrechtem Lichteinfall erfolgt. Abweichungen bis zu 5° von der 
Senkrechten wurden beobachtet. Die Vermutung, eine unterschiedliche Ober- 
flächenbeschaffenheit der Schichten könnte für dieses Verhalten der Grund 
sein, trifft nicht zu. Jedenfalls war kein deutlicher Zusammenhang mit der 
Oberflächenrauhigkeit festzustellen, etwa in der Weise, daß an rauheren 
Schichten eine größere Abweichung vom senkrechten Lichteinfall bei der 
Vorzeichenumkehr zu verzeichnen wäre. Im Gegenteil, glatte Schichten 
zeigten vielfach die Vorzeichenumkehr bei viel größeren Einfallswinkeln. 
Eines fällt jedoch auf, und dies konnte an allen untersuchten Schichten be- 
stätigt werden: Die schon bei senkrechtem Lichteinfall vorhandene EMK 
hatte stets das gleiche Vorzeichen wie das größere der beiden Maxima. Wir 
werden im Abschnitt 9auf diesen Tatbestand zurückkommen und einen wesent- 
liehen Grund dafür kennen lernen. 


1. Belichtung mit zwei Lichtquellen 


Nach dem Kurvenverlauf in Abb. 3 kann man selbstverständlich bei Be- 
liehtung mit zwei Lichtquellen unter zur Flächennormalen entsprechend ent- 
gegengesetzten Winkeln die von den einzelnen Lichtquellen herrührenden 
Photospannungen zu Null kompensieren. Dabei tritt eine Besonderheit auf: 
Wenn die Zelle mit einer Lichtquelle LZ, senkrecht bestrahlt wird, so daß die 
Photo-EMK jedenfalls Null ist und mit einer zweiten Lichtquelle L, unter 
einem bestimmten von Null verschiedenen Winkel, so ist bei Bestrahlung mit 
beiden Lichtquellen zugleich die Photo-EMK kleiner als bei Bestrahlung mi 
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L, allein. Diese Tatsache kann auf Grund der Sättigungseigenschaften bei 
höheren Lichtintensitäten (Abschnitt 6) befriedigend erklärt werden und ist 
zugleich ein Hinweis auf die Bedeutung der Oberflächenrauhigkeit der Schich- 
ten für die Photo-EMK. 


5. Vergrößerung der Empfindlichkeit durch rauhe Schichtunterlagen 


Für das Zustandekommen einer Photo-EMK mit der beschriebenen Rich- 
tungsabhängigkeit ist offenbar wesentlich, daß bei größeren Einfallswinkeln 
ein möglichst großer Anteil des Lichtes in die Schicht eindringen kann. Bei 
einer ideal ebenen Oberfläche nimmt aber der photoelektrisch unwirksame 
reflektierte Anteil mit wachsendem Einfallswinkel zu. An rauhen Schichten 
werden dagegen wegen der dann vorhandenen, bei größeren Einfallswinkeln 
zum Strahlengange senkrecht stehenden Mikroflächen günstigere Eindring- 
bedingungen herrschen, welche größere Photospannungen zur Folge haben. 
Dies hat sich im vollen Umfange bestätigt. — Es wurden PbS-Schichten auf 
vermittels Sandstrahlgebläse aufgerauhte Glasflächen niedergeschlagen und 
unter den gleichen Bedingungen formiert. An solchen Zellen konnten bei 
einer wirksamen Halbleiterfläche von 25 x 10 mm? und bei einer Beleuch- 
tungsstärke von 15000 Lux Spannungen bis zu 2 Volt erreicht werden. Die 

beträchtliche Größe der 

" Photospannung erklärt sich 
| | andererseits dadurch, daß 
diese als Summe an einer 

N sehr großen Anzahl hinter- 
einander geschalteter Ein- 

| zelelemente entsteht. Als 
rn | solche sind die kleinen 
Zacken der rauhen Ober- 
fläche anzusehen (8. 
Abb. 7b). Damit wird auch 
| die Tatsache verständlich, 
5th | daß eine Vergrößerung der 
Zellenfläche in Richtung 
T der Elektroden die Photo- 

02 0 60 8% spannung erhöht, eine Ver- 


Photospannung in 10?V 


10}— — 


960 6 


| Einfallswinkel in Grad größerung der Zellenfläche 
74 quer zu den Elektroden da- 
| RQ gegen lediglich den Strom. 
| 4 | | Die Photoströme sind 
| wegen des hohen Innen- 
2. m 6640 « widerstandes der Schichten 
1? | klein (Größenordnung pA). 

| n 955 
Intensitätsabhängigkeit 

6. In ensitatsa ngig 
| der Photo-EMK 

4123| Wegen der Schwierigkeit 


Abb. 4. Abhängigkeit der Photo-EMK vom Einfall. im parallelen Lichte eine 
winkel, Lichtintensität als Parameter definierte nicht selektive 
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Liehtschwächung durchzuführen und ausreichend hohe Intensitäten zu er- 
halten, wurde die Messung der Intensitätsabhängigkeit im schwach konver- 
genten Lichte vorgenommen. Es wurde zu diesem Zwecke eine Linsen- 
anordnung benutzt, wobei mit Hilfe einer Aperturblende die Lichtintensität 
variiert werden konnte. Als mittlerer Einfallswinkel wurde der Winkel 
zwischen der Achse des Licht- 

bündels und der Flächen- x» | 
normalen angesehen. Die da- 
durch entstehendeA bweichung 
von der Winkelverteilung der 
Photo-EMK, welche man bei 

einer Messung im parallelen 4% 
Liehte erhalten würde, lag bei 
der im Vergleich zur Länge 
des Strahlenganges kleinen 5 
Abmessung der Zelle unterhalb 
1%. Abb. 4 zeigt die auf diese 
Weise gemessene Winkelver- 
teilung der Photo-EMK mit 
der Lichtintensität als Para- 
meter und Abb. 5 die zuge- 
hörige Intensitätsabhängigkeit der Photo-EMK für verschiedene Einfalls- 
winkel. Die Farbtemperatur der Lichtquelle betrug 2865°K. Die Kurven 
zeigen, daß die Photo-EMK mit wachsender Lichtintensität einem Sättigungs- 
werte zustrebt. 


Photospannung in 10 Volt ——- 
% 
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Abb.5. Abhängigkeit der Photo-EMK von der 
Lichtintensität bei verschiedenen Einfallswinkeln 


7. Spektrale Empfindlichkeit 


Da es sich bei den untersuchten PbS-Schichten um die gleiche Grund- 
substanz handelt wie bei 
den PbS-Widerstandszellen, 
sollte man auch dement- 
sprechend die gleichen spek- 
tralen Eigenschaften erwar- 
ten. Diese sind kennzeich- 
nend für den Anregungs- 
mechanismus, welcher in 4 | \ 
beiden Fällen offenbar der pony 
gleiche ist. Die Messung, | 4 
welche mit einem Spiegel- | 
monochromator im energie- 
gleichen Spektrum durchge- “a 1 2 6 18 2 4 26 28 3 
führt wurde, bestätigt dies sagte 
(Abb. 6). Abb. 6. Spektrale Verteilung der Photo-EMK 
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8. Über die Ursache der richtungsabhingigen Photo-EMK 


Wegen der Richtungsabhängigkeit der Photo-EMK liegt die Vermutung 
nahe, daß es sich dabei um ein Photo-Diffusionspotential handelt, wie es erst- 
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Dr 7 a u Der Unterschied bestünde dann lediglich darin, daßi in unserem Falle eine be; 
272 a = kristalline Schieht quer zur Richtung der Elektroden bei stetiger Änderung 
; des Einfallswinkels belichtet wurde. Außerdem hat neuerdings Moss) 
Br auch an natürlichen PbS-Einkristallen Photo- Diffusionspotentiale festge- 
stellt. Eine quantitative Behandlung dieser Erscheinung, wie sie erstmalig 
von Frenkel4) und in neuester Zeit von Moss, Pincherle und Woodward) 
und gleichzeitig von Aigrain und Bulliard!%) gegeben wurde, geht davon 
aus, daß das Licht beim Eindringen in den Kristall auf seinem Absorptions- 
Be wege ein Konzentrationsgefälle von angeregten Elektronen und Defekt- 

_ elektronen hinterläßt, welches bestrebt ist, sich durch Diffusionsströme aus- 
=. _ zugleichen. Bei diesem Vorgange kommt auf Grund der unterschiedlichen 
Beweglichkeit der positiven und negativen Ladungsträger eine räumliche 
2 ie Ladungstrennung zustande, welche ein elektrisches Feld zur Folge hat. 

Dieses ist so gerichtet, daß es die Ladungsträger mit größerer Beweglichkeit 
bw verzögert und die mit der geringeren beschleunigt. Die auftretenden Ströme 
3 5 _ sind also teilweise Diffusionsströme und teilweise Feldströme. Es gilt danach 
für Elektronen 


je = —9 D, dnidx —qu,n E, 


und für Defektelektronen 
2 
in = D, dhidz + qu,hk, (2) 


und h Elektronendichte bzw. Defektelektronendichte, u, und u, Beweglichkeiten, 
q Elementarladung und D Diffusionskonstante). 


= Im stationären Falle müssen dn/dt und dh/dt verschwinden, wobei dn/dt und 
3 Kontinuitätsgleichung gegeben sind: 


dn/dt = — 1/qdivj, + Q, — R(n + h) 
=— div j, + R (n + h). 


optische Anregung neu und R den Rekombinations- 
faktor. Ebenso wie bei Frenkel ist hier zur Vereinfachung ein monomoleku- 
i dares Rekombinationsgesetz angenommen worden. Auf Grund der exponen- 
_ tiellen zeitlichen Abklingkurven beim PbS, wie sie von Scanlon, Petritz 
& _ und Lummis!) und von Pick?8) gemessen wurden, mag eine solehe Annahme 
berechtigt sein. Unter der weiteren vereinfachenden Annahme, daß Q für 
Bi x > 0 verschwindet, daß also alles Licht in einer extrem dünnen Oberflächen- 
_ schicht absorbiert wird, und unter Hinzunahme der Poissonschen Gleichung 
kann man die Gin. (1) bis (3) lösen und erhält fiir das Photopotential 


i 
> 


_ & —% kT Ny 
12) H. Dember, Physik. Z. 32, 554 u. 856 (1931). 
3) 'T. 8. Moss, Proce. physic. Soc. (B) 66, 993 (1953). Mi 


+ 


M4) J. Frenkel, Physik. Z. Sowjetunion 8, 185 (1935). 
a 8) T. S. Moss, L. Pincherle, A. M. Woodward, Proc. physic. Soc. London (B) 
66, 743 (1953). 

a 16) P. Aigrain u. H. Bulliard, C. R. Acad. Sci. Paris 286, 595 (1953). 
17) W. W. Scanlon, R. L. Petritz, F. L. Lummis, Physic. Rev. 86, 659 (1952). 
18) H. Pick, Ann. Physik 8, 255 (1948). 
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Dabei bedeutet n, die Ladungsträgerkonzentration an der Oberfläche und 
n diejenige im Inneren des Kristalles. 

Diese unter vereinfachten Voraussetzungen gewonnene Beziehung führt 
bei einer ersten überschläglichen Rechnung unter der Annahme, daß die 
Elektronenkonzentration im Inneren der Kristallite auf den e-ten Teil abge- 
sunken ist und im Hinblick darauf, daß sich die an einer Vielzahl von Einzel- 
elementen entstehenden Photospannungen addieren, zu einer den tatsächlichen 
Beobachtungswerten entsprechenden Größenordnung. 


Weitaus begründetere Schlußfolgerungen können aber aus dem Vorzeichen 
der Photospannung gezogen werden, und dies deshalb, weil es möglich ist, mit 
dieser Schlußweise eine Alternativfrage für den Anregungsmechanismus zu 
verbinden. Nehmen wir zunächst an, es handle sich beim Bleisulfid um eine 
jrundgitterabsorption, wobei Paare von Elektronen und Defektelektronen 
gebildet werden. Da die Beweglichkeit der Elektronen im PbS größer ist als 
die der Defektelektronen, müßte in diesem Falle nach Gl. (4) die der Licht- 
quelle abgewandte Seite ein negatives Potential annehmen. Dies würde mit 
dem experimentellen Befunde übereinstimmen. Es erhebt sich aber sofort die 
Frage, welche Rolle der für die Empfindlichkeit der Schichten wesentliche 
Sauerstoff spielt. Fest steht jedenfalls die Tatsache, daß die Schichten nach 
der chemischen Herstellung ohne weitere Formierung keine lichtelektrischen 
Eigenschaften zeigen !®). 

Im Falle einer Störstellenabsorption, wobei der Sauerstoff als Akzeptor 
wirken müßte, wäre nach der obigen Deutung der Diffusionsspannung die bei 
der Ladungstrennung entstehende elektronegative Komponente als ioni- 
siertes Sauerstoffatom im Gitter lokalisiert. Es könnten also nur die allein 
beweglichen Defektelektronen in das Innere des Kristalles hinein diffundieren, 
d.h. die der Lichtquelle abgewandte Seite müßte ein positives Potential an- 
nehmen. Dies ist aber im allgemeinen nicht der Fall. Der Sauerstoff kann 
also nieht einfach die Funktion ins Gitter eingelagerter Störstellen haben 
Es scheint nach diesen Überlegungen ausgeschlossen zu sein, daß wir es bei 
der beschriebenen Photospannung mit einem Diffusionspotential zu tun haben. 

Dagegen scheint folgende Deutung des Effektes mit den Beobachtungs- 
tatsachen gut übereinzustimmen: Es ist ohne weiteres möglich, daß der Sauer- 
stoff während der Formierung nicht vorzugsweise in die Schicht eindringt, 
sondern in Form eines chemisorbierten Belages an der Oberfläche haften bleibt. 
(Wir machen hiermit für die gesamte Schichtoberfläche eine Annahme, wie sie 
für die inneren Grenzflächen der einzelnen Körner von Schwartz, Gibson 
und Smith®)?)8) zur Erklärung der Photoleitfähigkeit diskutiert worden ist.) 
Unter der Voraussetzung, daß die Einwirkung des Sauerstoffes auf die Ober- 
fläche beschränkt bleibt, kann man weiterhin annehmen, daß das Innere der 
Kristallite unter quasi Vakuumbedingungen formiert wurde, also überschuß- 
leitend ist. Der sehr stark elektronenaffine Sauerstoff würde dann mit der 
Halbleitersubstanz derart in Wechselwirkung treten, daß Elektronen unter 
Zurücklassung positiver Störstellenreste im Inneren der Kristallite in die 
Sauerstoffoberflächenzustände wandern. Die Oberfläche wird sich dadurch 
gegenüber dem Inneren negativ aufladen, und es entsteht eine mit einem Rand- 


19) Andererseits konnte festgestellt werden, daß auch im Vakuum getemperte PbS- 
Schichten keine oder nur geringe Photoempfindlichkeit aufweisen. 
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felde siibiniiiiits elektrische Doppelschicht, ähnlich wie bei einem Halbleiter- 
Metallkontakt (siehe Abb. 7). 


Bei Belichtung werden Paare von Elektronen und Defektelektronen ge- 
bildet, welche in dem bestehenden Randfelde getrennt werden und auf diese 
Weise eine Photo-EMK verursachen. Diese Deutung liefert das mit der Beob- 
achtung übereinstimmende Vorzeichen und macht darüber hinaus die 
Bedeutung des Sauer- 
stoffes verständlich; sie 
steht ferner im Einklang 
mit der Feststellung von 
Avery”), daß sich eine 
Sauerstoffbehandlung der 
PbS-Schichten nicht 
wesentlich im Absorp- 
tionsspektrum bemerkbar 
macht; denn wir können 
annehmen, daß die durch 
den Sauerstoff hervor- 
gerufenen Oberflächen- 
Abb. 7. a) Energieschema einer durch Sauerstoff- terme schon im unbelich- 
Oberflächenzustände (SO) hervorgerufenen elektrischen teten Zustande und bei 
Doppelschicht; D Donatorenterme. b) Vereinfachte 
Darstellung einer rauhen PbS-Schicht bei schräger Be- normalen Temperaturen 

lichtung vollständig mit Elek- 
tronen besetzt sind. 

Von Bardeen”) ist darauf hingewiesen worden, daß auch bei Halbleiter- 
Metallkontakten schon ohne die Anwesenheit des Metalles vorhandene Ober- 
flächenzustände eine entscheidende Bedeutung haben können. Ein direkter 
Nachweis ist aber immer dann schwierig, wenn der Halbleiter durch die Me- 
tallelektrode hindurch belichtet werden muß, wie es bei den üblichen Photo- 
elementen der Fall ist. Bei der hier benutzten Beobachtungsart steht es aber 
außer Zweifel, daß die Metallkontakte hinsichtlich der Entstehung der Photo- 
EMK völlig unwesentlich sind. 


Für den Anregungsmechanismus im PbS würde aus dieser Deutung folgen, 
daß wir entgegen anderen Auffassungen®) mit einer Grundgitterabsorption 
entsprechend einem Bandabstand von 0,4eV zu rechnen haben; denn das 
paarweise Auftreten von Elektronen und Defektelektronen bei der Licht- 
absorption ist eine wesentliche Voraussetzung für die Theorie des Sperr- 
schicht-Photoeffektes (siehe Lehovee®)). 


Was den Schichtenaufbau anbetrifft, so geben die hier mitgeteilten Ver- 
suchsergebnisse kein Anzeichen von der Existenz statistisch verteilter ,,p-n- 
Junctions‘. Sie scheinen vielmehr zu zeigen, daß die gesamte Oberfläche der 
PbS-Schicht bei den gewählten Formierbedingungen eine einheitliche Struktur 
erhält. Wenn wir die inneren Oberflächen der einzelnen Kristallite mit ein- 


20) D. G. Avery, Proc. physic. Soc. (B) 66, 134 (1953). 

21) J. Bardeen, Physic. Rev. 71, 717 (1947). Er 
2) L.Genzel, H. Müser, Z. Physik 127, 194 (1950). f 3 
») K. Z. Lehovec, Z. Naturforschg. 1, 258 (1946) u. 2a, 398 (1947. 
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beziehen und an einem Sperrschichtenmodell zur Erklärung der Photoleit- 
fähigkeit festhalten wollen, so haben aus diesem Grunde die von Smith?) 
diskutierten Schottky-Sperrschichten vor den „p-n-junetions‘“ den Vorzug. 


9. Formierung der Schichten in einem Temperaturgradienten 


In Abschnitt 3 wurde schon darauf hingewiesen, daß die Unsymmetrie 
in der Winkelverteilung der Photo-EMK (ungleiche Höhe der beiden Maxima 
und die damit verbundene EMK bei senkrechtem Lichteinfall) nicht mit der 
geometrischen Oberflächenbeschaffenheit der Schichten zusammenhängt. 
Eine besonders gute Symmetrie in der Winkelverteilung der Photo-EMK 
ließ sich dagegen durch das Einhalten möglichst homogener Temperatur- 
bedingungen während der Formierung erreichen. (Die Schichten wurden zu 
diesem Zwecke mit einem Kupferzylinder umgeben, welcher für eine gleich- 
mäßige Wärmeabstrahlung sorgte) Man könnte nach alledem vermuten, 
daß die Unsvmmetrie der Photo-EMK durch die Überlagerung einer weiteren 
Photo-EMK zustande kommt, deren Vorzeichen nicht von der Lichtrichtung 
abhängt, sondern welche an die Struktur der Schichten gebunden ist, wie es 
etwa bei mikrokristallinen ‚‚p-n-junetions“ der Fall ist. 


Um diese Vermutung nachzuprüfen, wurden PbS-Schichten in einem 
Temperaturgradienten formiert. Die Schichtträger (Glasplatten der Größe 
2,7 x 7,5em?) wurden da- 
bei an der Stirnseite des 
Formierofens von einer 
wassergekühlten Halterung 
K und am anderen Ende 
von einem auf der Tem- = 
peratur des Ofens (600° C) 
befindlichen Metallklotz M 
gehalten (Abb. 8). Aus dem 
ersten Drittel der Schicht- 
linge, also aus einem Be- 
reiche, wo bei zugleich ausreichender Höhe der Gesamttemperatur der Tem- 
peraturgradient am größten sein dürfte, wurden Stücke der Größe 
10 x 20 mm? herausgeschnitten und auf die Winkelverteilung der Photo- 
EMK hin untersucht. Dabei ergab sich folgendes: 

Während sich bisher die PbS-Schichten hinsichtlich der Richtungsab- 
hängigkeit der Photo-EMK innerhalb der Meßgenauigkeit vollständig isotrop 
verhielten, hängt der Verlauf der Photo-EMK jetzt von der Lage der Einfalls- 
ebene des Lichtes ab. Greift man die Photospannung in derjenigen Richtung 
an der Schicht ab, in welcher der Temperaturgradient wirkte, und enthält 
die Einfallsebene des Lichtes diese Richtung, so stellt man eine stark unsym- 
metrische Winkelverteilung der Photo-EMK fest (Kurve a in Abb. 9). Es 
wurden sogar häufig Schichten gefunden, bei denen es zu gar keiner Vor- 
zeichenumkehr kam (Kurve a in Abb. 10). Das Vorzeichen der schon bei 
senkrechtem Lichteinfall vorhandenen EMK war stets so, daß die während der 
Formierung höher temperierte Seite ein negatives Potential annahm. Wurde 
die Einfallsebene und entsprechend die Richtung der Elektroden jedoch senk- 
recht zur Richtung des Temperaturgradienten gelegt, so wurde die normale 
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Abb. 8. Anordnung zur Formierung der PbS-Schichten 
in einem Temperaturgradienten 
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Winkelabhängigkeit der Photo-EMK beobachtet mit nur geringfügiger Un- 


symmetrie. 


Wenn wir daran festhalten, daß die durch stark inhomogene Temperatur- 
verteilung entstehende Unsymmetrie der Photo-EMK als eine Überlagerung 


jy) 


oto-EMK [8 x10 


Ph 
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4 Einfallswinkel in Grad _ 


Abb. 9. Winkelverteilung einer unter der Ein- 
wirkung eines Temperaturgradienten formierten PbS- 
Schicht. Die Elektroden wurden in Richtung des 
Temperaturgradienten angebracht. a) Einfallsebene 
enthält die Richtung des Temperaturgradienten. 
b) Einfallsebene senkrecht zur Richtung des Tem- 
peraturgradienten. c) Differenz der Kurven a und b 


zweier photoelektromotori- 
scher Kräfte zu erklären ist, 
so ergibt sich eine einfache 
Möglichkeit, diese beiden 
Photospannungen in ihrer 
Abhängigkeit vom Einfalls- 
winkel getrennt zu erfassen, 
Es wurden zu diesem Zwecke 
die Elektroden in der Rich- 
tung an der PbS-Schicht 
angebracht, in welcher der 
Temperaturgradient wirkte 
und die Einfallsebene einmal 
in die gleiche Richtung ge- 
legt und das andere Mal 
senkrecht dazu. Im ersten 
Falle wird die Summe beider 
Photospannungen gemessen 
(die Kurven a der Abb. 9 
und 10), im zweiten Falle 
dagegen die dem Vorzeichen 


nach nicht riehtungsabhängige Photo-EMK allein (die Kurven b der Abb. 9 
und 10), welche sich nach einem cos-Gesetz mit dem Einfallswinkel ändert. 
Vor allem bei der letzten Messung mußte sorgfältig darauf geachtet werden, 
daß die Einfallsebene ihre zur Richtung des Temperaturgradienten senk- 


Photo-EMK [8x10?v] 


Einfallswinkel in Grad 


Abb. 10. Winkelverteilung einer unter der Einwirkung 
eines Temperaturgradienten formierten PbS-Schicht, 
welche bei keinem Einfallswinkel eine Vorzeichen- 
umkehr aufweist. Die Kurven a, b, ¢ haben die ent- 
sprechende Bedeutung wie in Abb. 9 


rechte Lage bei allen Ein- 
fallswinkeln unverändert 
beibehielt. Anderenfalls 
würde stets eine Kompo- 
nente der dem Vorzeichen 
nach richtungsabhängigen 
Photo-EMK mit erfaßt. 
Ebenso war es nötig, die 
unverspiegelten Seiten- 
kanten der Glasunterlage 
mit schwarzem, licht- 
undurchlässigem Nitrolack 
zu schließen, damit kein 
störendes Nebenlicht auf 
die Rückseite der PbS- 
Schicht treffen konnte. 
Durch Differenzbildung der 
Ordinatenwerte der Kur- 
ven a und b erhält man die 
Kurven e, welche der dem 
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Vorzeichen nach richtungsabhingigen Photo-EMK zuzuordnen sind. Diese 


Un- 
Kurven zeichnen sich in der Tat durch eine recht gute Symmetrie aus. 


atur- Die genannten Messungen wurden an Schichten ausgeführt, deren Glas- 
rung unterlagen nicht extra aufgerauht waren. Der gleiche Sachverhalt ist zwar 
)tori- auch an aufgerauhten Schichten zu beobachten, nur ist dann die entstehende 
n ist, Unsymmetrie in der Winkelverteilung nicht so groß und dies deshalb, weil an 
fache rauhen Schiehten die dem Vorzeichen nach riehtungsabhängige Photo-EMK 
»iden stets viel größer ist als die sich ihr überlagernde Struktur-EMK. 

ihrer Von den PbS-Schichten, welche bei allen Einfallsrichtungen nur ein ein- 
falls- ziges Vorzeichen der Photo-EMK aufwiesen, zeigten die meisten schwache 
ssen, Gleichrichtereigenschaften (Verhältnis der Ströme in Durchlaß- und Sperr- 
ecke Richtung wie 2:1 bei Spannungen bis zu 2 Volt). Dies spricht dafür, daß 
tich- die Struktur-EMK in interkristallinen Sperrschichten ihre Ursache hat. Man 
icht kann sich vorstellen, daß diese durch eine mit unterschiedlicher Beweglichkeit 
der von Pb- und O-Störstellen erfolgende Thermodiffusion im Temperaturgefälle 
irkte während der Formierung zustande kommen. Die schweren Pb-Atome werden 
mal sich dabei an der höher temperierten Seite anreichern und es können auf 
5 Be diese Weise mikrokristalline ‚‚p-n-junetions‘‘ entstehen. Mit dieser Deutung 
Mal wäre die Beobachtung zu erklären, daß die während der Formierung höher 
sten temperierte Seite bei Belichtung ein negatives Potential annimmt, wie es dem 
sider Pb-Überschuß und damit der überschußleitenden Seite der ,,p-n-junction“ 
ssen entsprechen würde. Aus der Tatsache, daß diese Photo- EMK bei senkrechter 
b. 9 Belichtung mit dem gleichen Vorzeichen auch dann vorhanden ist, wenn man 
‘alle nicht die ganze Halbleiterfläche ausleuchtet sondern nur kleine örtlich be- 
*hen grenzte Bereiche, kann man schließen, daß die Störstellendiffusion nicht über 
rb. 9 sehr weite Strecken der Schicht erfolgt, sondern innerhalb kleiner Bezirke etwa 
lert. an den Korngrenzen zum Stillstand kommt. Auf diese Weise können eine 
den, große Anzahl in Serie liegender Mikrosperrschichten entstehen. Es muß jedoch 
enk- hier noch die Frage offen bleiben, ob es sich dabei wirklich um ,,p-n-junctions* 
Ein- handelt. Es ist immerhin denkbar, daß durch Thermodiffusion auch unsym- 
dert metrische Schottky-Sperrschichten an den Korngrenzen entstehen. Zu- 
falls sammenfassend kann man aber auf Grund der Versuchsergebnisse folgendes 
1po- sagen: Bei der dem Vorzeichen nach richtungsabhängigen Photo-EMK haben 
hen wir es mit einem Vorgange zu tun, für welchen die äußeren Oberflächen der 
igen Kristallite entscheidend sind, bei der zusätzlichen Photo-EMK dagegen, 
abt. welche die Unsymmetrie bewirkt, mit einem Vorgange an den inneren Ober- 
= flächen oder im Inneren der Kristallite selbst. ; 

en- 

lage 

ht: 10. Alterungserscheinungen 

ack Nach der Formierung zeigten die Bleisulfidschichten im allgemeinen kein 
ein konstantes Verhalten. In manchen Fallen nahm die Empfindlichkeit im Ver- 


auf laufe von Wochen von etwa 2 Volt auf einen konstanten Endwert einiger 
bS- zehntel Volt bei 15000 Lux ab. Daneben wurden aber auch Zellen gefunden, 
ate. welche kurz nach der Herstellung wenig empfindlich waren und deren Emp- 
der findlichkeit in einem Zeitraume von 4 bis 6 Wochen auf einen konstanten End- 
ur- wert von 1 bis 1,5 Volt bei der gleichen Beleuchtungsstärke anstieg. Darüber 
die hinaus machte sich bei diesen Zellen nach längerer Zeit wieder ein gering- 
lem fügiger Rückgang der Empfindlichkeit bemerkbar. Der Grund für dieses 
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unterschiedliche Verhalten und die Zusammenhänge mit den — 
— sind noch nicht zu übersehen und bedürfen einer näheren Unter- 


Für die Anregung, auf diesem Gebiete zu arbeiten, und für wertvolle 
und fördernde Diskussionen bin ich Herrn Prof. Dr. P. Görlich zu Dank ver- 
pflichtet. Herrn Prof. Dr. W.Schütz danke ich für sein wohlwollendes Inter- 
esse an der Arbeit. 


Jena, Physikalisches Institut der Universitat. _ 
Bei der Redaktion eingegangen am 1. August 1955. 
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& & iwandstabilisierten Argonbogen') 


Von B. Winde 


Mit 5 Abbildungen ne 


Inhaltsübersie 


Aus Spannungsmessungen an langgestreckten Lampen mit ; tee 
lichem Elektrodenabstand wird die Summe von Anoden- und Kathodenfall 
beim Argonbogen zu 13 + 1 Volt gemessen. Sie hängt im Untersuchungs- 
bereich (300... 1100 Torr, 5... 40 Amp.) nicht von Stromstärke und Druck 
ab. Die Lampen werden an der Pumpe (Fülldruck = Brenndruck) betrieben. 
Die Ergebnisse werden mit denen von Bauer und Schulz an Xenon sowie 
von Busz und Finkelnburg am Argon-Hochstrombogen verglichen. 


Über die Bogenentladung in Edelgasen vom wandstabilisierten Typ sind 
bisher relativ wenig experimentelle und theoretische Arbeiten von physikali- 
scher Seite veröffentlicht worden. Gerade dieser Bogentyp bietet aber für 
die theoretische Behandlung gewisse Vorteile, da die Verhältnisse hier noch am 
übersichtlichsten sind und durch eine nicht zu komplizierte Leistungsbilanz 
erfaßt werden können. Versuche in dieser Richtung haben gezeigt, daß man 
bei Edelgasbögen nicht ohne weiteres von der Elenbaas-Hellerschen Dif- 
ferentialgleichung ausgehen kann, in der Form, wie sie beim wandstabili- 
sierten Quecksilberbogen benutzt wird. Man muß, um in ungefähre Über- 
einstimmung selbst mit groben Experimenten?) zu kommen, entweder die 
Säule des Edelgasbogens (die Untersuchungen wurden speziell an Xenon 
durchgeführt) als schwarzen Strahler auffassen*); diese Annahme ist aber 
aus mehreren Gründen unwahrscheinlich und wird vom Experiment nicht 
bestätigt oder man muß gewisse Annahmen über die Temperaturabhängig- 
keit der Wärmeleitfähigkeit, Beweglichkeit usw. machen®). Unsere Kenntnisse 
über diese Abhängigkeit bei Temperaturen in der Größenordnung 10000° 
sind jedoch recht lückenhaft und deshalb scheinen Untersuchungen, die aus 
dem Gradienten und seiner Abhängigkeit von der Leistung Rückschlüsse auf 
diese Größen erlauben, von Bedeutung zu sein. 


!) Die vorliegende Arbeit ist der etwas erweiterte 1. Teil der Dissertationsarbeit des 
Verfassers an der mathem.-naturwissenschaftlichen Fakultät der Humboldt-Universität 
Berlin. (Wiss. Z. d. Humboldt-Universität, im Druck.) 

*) K. Larché, Z. Physik 136, 74 (1953). 

3) H. Schirmer, Z. Physik 136, 87 (1953); Z. angew. Physik 6, 3 (1954); Techn. 
wiss. Abh. d. Osram-Konzerns 6, 20 (1953). 

4) H. Schirmer, Vorträge auf der Tagung der Deutschen Physikal. Ges. in Hamburg 
(1954) und auf dem Symp. of Electr. Disch. in Delft (1955). 
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In dieser Arbeit soll zunächst das Problem des Kathoden- und Anodenfalls 
an Hand experimenteller Untersuchungen an wandstabilisierten Argonent- 
ladungen diskutiert werden, um damit Unterlagen für die Berechnung des 
Gradienten zu gewinnen. 


1. Versuchsanordnung und Meßverfahren 


Die gestreckten Lampen brannten waagerecht. Sie sind nach dem Prinzip 
der XF-Lampen gebaut). Die thorierten Wolframelektroden sind projektil- 
förmig ausgebildet und befinden sich an Wolframstäben, die in Hartglas einge- 
schmolzen sind. Das eigentliche Entladungsrohr aus Quarz ist mit Hilfe von 
Übergangsgläsern an diese 
Durchführungen geschmol- 
zen. Die Messungen wurden 
an der Pumpe durchgeführt. 
Abb. 1 zeigt einige abge- 
schmolzene Lampen. 

Um die Summe von 
Kathoden- und Anodenfall 
zu bestimmen, wurden 
Strom - Spannungscharakte- 
ristiken an Lampen mit 
unterschiedlichem Elek- 

— trodenabstand (20, 40, 60, 

85, 90, 95, 100mm) bei 

gleichem Rohrdurchmesser 

Abb. 1. Lampen mit verschiedenem Elektrodenabstand (11mm) in einem Druck- 

und Durchmesser vor dem Ansetzen an die Pumpe bereich von 300... 1100 Torr 

aufgenommen. Aus dem 

unterschiedlichen Spannungsbedarf bei gleichem Strom und Druck läßt sich 

nach Elenbaas®) durch Extrapolation auf den Abstand Null der Spannungs- 
bedarf für Kathoden- und Anodengebiet bestimmen. 


Die Messungen wurden in einem Stromstärkebereich zwischen 5...40A 
durchgeführt. Das entspricht Leistungen/em Rohrlänge von 40... 120 W/em. 
(Die vorliegenden ähnlichen Messungen von Elenbaas’) sowie von Bauer 
und Schulz‘) wurden bei kleineren Leistungen bzw. Stromstärken durchge- 
führt.) Die Lampen wurden mit Reinstargon bei verschiedenen Kaltdrucken 
gefüllt. Das Reinstargon wurde vom VEB Berliner Glühlampen-Werk in 
Stahlflaschen bezogen. Spektroskopisch konnten — bis auf ganz gering- 
fügige Quecksilberreste, die von der Pumpe herrühren — keine Verunreini- 
gungen festgestellt werden. 


Das Entladungsrohr muß im Betrieb mit Wasser gekühlt werden. Es 
befindet sich zu diesem Zweck in einem speziellen, von Leitungswasser durch- 


5) Herrn A. Ihln vom Berliner Glühlampen-Werk und dem VEB Werk für Fern- 
meldewesen „WF“ bin ich zu Dank verpflichtet für den Bau der Lampen und das Ent- 


gegenkommen bei Sonderwünschen. 
6) W. Elenbaas, Physica 3, 947 (1936). 
7) W. Elenbaas, Phil. Res. Rep. 4, 221 (1949). 
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ulls flossenen Kühltrog aus Glas, der an den Enden mit Gummipfropfen ver- 
Bi schlossen ist. Es zeigt sich, daß die Kühlung innerhalb der Meßgenauigkeit 
les keinerlei Einfluß auf die Entladungsparameter im stationären Zustand hat. 


Fehlt die Kühlung, oder 
wird nur schwach gekühlt, 
so macht jedoch die Zündung 
der Lampe Schwierigkeiten. 


Die Betriebsspannung und 
u die Zündspannung werden 
ne über Klemmschrauben an 
Ber: die Wolframstäbe gelegt. 
é Die fertig montierte Lampe 
w zeigt Abb. 2. Die Lampen 
vr wurden in üblicher Weise 
oi evakuiert, formiert und ge- 
x füllt. Ein Beispiel für die 
ou Meßergebnisse zeigt Abb. 3. 


all Da die Lampe an der 
Pumpe betrieben wurde, Abb. 2. Lampe an der Pumpe betriebsfertig montiert 
stieg der Argondruck in 


Ri der Lampe infolge des relativ großen Volumens der Pumpanlage im Betrieb 
. nieht wesentlich an, in keinem Fall mehr als 5%. Ein zunächst vorgesehenes 
60 Totvolumen für den Druckausgleich konnte später fortgelassen werden. Die 
bei oben angegebenen Elektroden- ww 
abstände sind die vom Her- au 
ser 
a steller angegebenen SollmaBe. 30 
wi Die wirklichen Elektroden- 2 
. . r 
al abstände differieren um einige 
ich mm und wurden durch 50} 
gs- Projektion zusammen mit sot 
einem Maßstab von Fall zu 
Fall ermittelt. “or 
‚A Parallelmessungen an je- „| 
weils zwei gleichartigen Lam- 
ler pen zeigten, daß die Meß- 2 
ge- ergebnisse befriedigend repro- ‚ot 
kon duzierbar waren. Jeder Wert 
wurde mindestens siebenmal 5 0 2 30 35 
gemessen. Aus den gemessenen EIER ‘ 
ni- wende der Mittelnert Abb.3. Beispiel einer gemessenen Strom-Spannungs- 
charakteristik 
gebildet, der im folgenden be- 
Es nutzt wird. Der mittlere Fehler betrug in dem stark fallenden Bereich der 
ch- Charakteristik weniger als 5%, in dem nahezu waagerechten Bereich (bei 
größeren Stromstärken) weniger als 3%. 
int- 2. Die MeBergebnisse 


Es ist mit Recht schon von Engel und Steenbeck°) und kürzlich von 
Bauer und Schulz?) darauf hingewiesen worden, daß der auf die oben ange- 


. ®) A.v. Engel u. M. Steenbeck, Elektrische Gasentladungen Bd. II, Berlin 1934 
153. 
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gebene Weise ermittelte Wert nicht eigentlich die Summe von Anoden- und 
Kathodenfall ist, weil in den — unter Umständen stark — kontrahierten 
Übergangsgebieten vor den Elektroden die Feldstärke viel höher ist als in der 
ungestörten Säule. Es hat sich aber doch wohl eingebürgert, beim Bogen die 


407 


20 8 
°p=300 mmHg 
20 W 60 60 00 le 
Abb. 4. Abhängigkeit der Lampenspannung vom 


Elektrodenabstand bei konstanter Stromstärke; Druck 
als Parameter 
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60 
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Abb. 5. Abhängigkeit der Lampenspannung vom 

Elektrodenabstand bei konstantem Druck; Stromstärke 

als Parameter. Die Abbildung zeigt die Verhältnisse 
bei kleinen Stromstärken (i < 15 Amp.) 


/mm] 
100 Ip 


Definition etwas anders zu 
wählen als bei den analogen 
Verhältnissen inder Glimm- 
entladung. Der voll ausge- 
bildete Bogen besteht aus 
Fallbereichen, Übergangs- 
gebieten und Säule, wobei 
die Kontraktion im Über- 
gangsgebiet unter Um- 
ständen bis weit in die 
Säule wirken kann, ja 
manchmal bei sehr kurzen 
Bögen den gesamten Habi- 
tus des Bogens bestimmt 
(siehe z. B. 10)). 

Hat man einen voll aus- 
gebildeten, wandstabilisier- 
ten Bogen, der das Rohr 
ganz erfüllt, so wird man 
aber zur Gradientenberech- 
nung den zusätzlichen 
Spannungsbedarf der Kon- 
traktionsgebiete ebenfalls 
von der Gesamtspannung 
abziehen müssen und auch 
für technische Probleme er- 
scheint uns die angegebene 
einfache Methode am ge- 
eignetsten zu sein. Wie 
wir im weiteren zeigen 
werden, ist wohl auch die 
von Bauer und Schulz 
geäußerte Meinung, daß 
dem aus der Spannung— 
Elektrodenabstand — Gera- 
den ermittelten Wert 


keinerlei tiefere physikalische Bedeutung zukommt, etwas überspitzt. Die 
Ergebnisse unserer Messungen sind in Abb. 4 und 5 dargestellt. Die Summe 
von Kathoden- und Anodenfall im angegebenen Sinn beträgt für den Argon- 


bogen 13 + 1 Volt. 


Busz und Finkelnburg!!) haben bei stromstarken Bögen in strömendem 
Argon (50...500A) bei Atmosphärendruck durch Extrapolation aus 


10) B. Winde, Ann. Physik (6) 15, 57 (1954). 


1) G. Busz u. W. Finkelnburg, Z. Physik 188, 212 (1954). Br = 
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Messungen bei verschiedenem Elektrodenabstand die Summe von Kathoden- 
und Anodenfall zu 5,5 + 1 Volt bestimmt. Bauer und Schulz‘) fanden 
aus Messungen an sehr kurzen Xe-Bögen einen Anoden- und Kathodenfall 
von 6,5 + 0,75 Volt (Grenzwert bei Stromstärken i > 10 A). 

Diese Ergebnisse sind offenbar einfach dadurch zu erklären, daß bei kleinen 
Elektrodenabständen die Spannung stärker mit abnehmendem Elektroden- 
abstand sinkt als in der Säule, weil in den Kontraktionsgebieten die Feld- 
stärke größer ist als in der Säule. Extrapoliert man in Abb. 4 der Arbeit!!) 
von größeren Elektrodenabständen geradlinig, so kommt man ebenfalls auf 
etwa 13 Volt. Die Busz-Finkelnburgschen Ergebnisse deuten also wohl 
nieht auf einen besonderen ,,Hochstromeffekt**. Es wäre jedoch nicht ausge- 
schlossen, daß bei den sehr kleinen Elektrodenabständen (unterhalb 2—3 mm) 
die Anode in die — bei den von Busz und Finkelnburg benutzten Strom- 
stärken unter Umständen beträchtlich große — kathodische Plasmakugel ein- 
taucht und dadurch ein ganz anderer Entladungstyp entsteht, der die kleinen 
Werte verursacht. In diese Richtung deutet auch die spektroskopisch ermittelte 
sehr hohe Achsentemperatur von etwa 30000°. Aus Potentialsondenmessun- 
gen ermittelten die Autoren am gleichen Objekt für die Elektrodenverluste 
einen Wert von 6 + 1,5 V. 

Vom Druck sind Kathoden- und Anodenfall in dem von uns unter- 
suchten Bereich nicht abhängig (Abb. 4). Wie Abb. 5 zeigt, erhält man brauch- 
bare Werte nur bei großen Stromstärken. Bei Stromstärken /< 15A ist 
die Stromdichte bei den verschiedenen Elektrodenabständen nicht gleich, da 
sich bei kleineren Abständen der Bogen von der Wand löst und die kontra- 
hierten Elektrodengebiete starken Einfluß gewinnen. Da aber auch die Wand 
weiterhin wirksam bleibt, werden die Verhältnisse recht unübersichtlich. 
Aus den vorliegenden Meßergebnissen ist demnach kein klares Bild über eine 
etwaige Stromstärkeabhängigkeit des Kathoden- und Anodenfalls zu gewinnen, 
wie sie von Bauer und Schulz in Xenon gemessen worden ist. Immerhin 
läßt die — allerdings recht grobe — Extrapolation aus den Werten von langen 
Röhren (80—100 mm) bei kleineren Stromstärken (6... 10 A) annehmen, 
daß eine wesentliche Änderung in dem vorliegenden Stromstärkegebiet nicht 
vorhanden ist. Das wird noch durch die Tatsache gestützt, daß die Werte in®) 
erst bei Stromstärken unterhalb 10 A zu wachsen beginnen. 

Zur Klärung der Frage, ob der Kathodenfall von der Stromstärke abhängt, 
wäre es notwendig, die Spannungs-Bogenlängenabhängigkeit bei konstanter 
Stromdichte zu extrapolieren. Zumindest der visuelle Querschnitt müßte 
also in jedem einzelnen Fall gemessen werden. Es genügt aber, zur Berech- 
nung des Gradienten in allen Fällen, wo es sich um lange, wandstabilisierte 
Bögen handelt, 13 Volt von der Elektrodenspannung abzuziehen. Änderungen 
von + 1 Volt geben an unseren Lampen bei der Feldstärkeberechnung einen 
maximalen Fehler von 7%. 

In welchem Verhältnis sich der Elektrodenverlust auf Anode und Kathode 
aufteilt, läßt sich aus den vorliegenden Messungen nicht entscheiden. Man 
wird jedoch keinen großen Fehler machen, wenn man in Analogie zu ®) an- 
nimmt, daß der Spannungsabfall im Kathodengebiet etwa 12 Volt, der im 
Anodengebiet etwa 1 Volt beträgt. 

Schließlich bleibt die Frage, in welcher Weise sich der Kathodenfall auf 
das Raumladungsgebiet und die Kathodische Plasmakugel aufteilt. Wenn 
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unsere Deutung der von Busz und Finkelnburg gemessenen kleinen Werte 
richtig ist, so sollte das Raumladungsgebiet etwa 5—6 Volt oder weniger, der 
Kontraktionsbereich mit der Plasmakugel rund 6—7 Volt benötigen. Diese 
Werte sind in guter Übereinstimmung mit den theoretisch geforderten Werten 
von Ecker!) und den experimentellen Ergebnissen von Bauer und Schulz‘), 

Die Arbeit wurde im Rahmen eines Forschungsauftrages des Zentralamtes 
für Forschung und Technik der Staatlichen Plankommission der Deutschen 
Demokratischen Republik durchgeführt. Herrn Prof. Dr. R. Rompe danke 
ich für die Anregung dieser Arbeit und für zahlreiche fördernde Diskussionen. 
Fräulein E. Kaczmarek danke ich für ihre ausdauernde Unterstützung bei 


den oft langwierigen Messungen. 2 
2) W. Ecker, Tagungsbericht „Gasentladung‘“, Berlin 1954, S. 148. s 


Berlin, Institut fiir Strahlungsquellen der Deutschen Akademie der 
Wissenschaften. 
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Ebene elektrostatische Felder, die eine strenge Berechnung 
der Elektronenbahnen zulassen 


Inhaltsiibersicht 


Es werden ebene elektrostatische Felder gesucht, die eine Bestimmung 
der Trägerbahnen durch Quadraturen zulassen. Dazu ist es nötig, die von 
Staeckel angegebenen Feldtypen weiteren Bedingungen zu unterwerfen. 


Die sich werden explizit 


Felder, die eine strenge Berechnung der Elektronen- oder Ionenbahnen 
zulassen, sind für die Partikeloptik von großer Bedeutung. Sie können als 
Modelle für reale Felder dienen, wobei man an den Modellen die Abhängigkeit 
der optisch wichtigen Größen (Brennweiten, Fehlerkoeffizienten usw.) von 
den Parametern der Bahnen und des Feldes leicht überblicken kann. Bisher 
sind nur wenige derartige Feldmodelle gefunden worden. Im ganzen Raum 
konnte bisher nur der Fall des homogenen Feldes, des (rotationssymmetrischen) 
Hyperbel- und des Parabelfeldes!) beherrscht werden. Darüber hinaus war nur 
ein streng integrierbarer Fall bekannt, der sich allerdings nur auf das paraxiale 
Gebiet bezog?). Es ist dann P. Schiske gelungen, ein im ganzen Raum 
streng behandelbares Einzellinsen-Feld zu finden?). 

Hier soll der Versuch unternommen werden, das Problem allgemeiner und 
systematischer in Angriff zu nehmen. Es handelt sich dabei um folgende Auf- 
gabe: Es sind elektrostatische Felder zu suchen, welche so beschaffen sind, 
daß man die Bahnen der sich in ihnen bewegenden geladenen Partikel streng 
berechnen kann. Der Begriff ‚streng‘ ist allerdings verschwommen und 
soll dahingehend präzisiert werden, daß jene Bahnen als streng berechenbar 
bezeichnet werden sollen, die durch eine einfache Quadratur gewonnen werden 
können (führt die Integration dabei auf bereits tabuliert vorliegende oder gar 
auf elementare Funktionen, so ist das natürlich willkommen). Eine Beschrän- 
kung auf das paraxiale Gebiet soll nicht vorgenommen werden‘). Es ver- 

1) A. Recknagel, Z. Physik 104, 381 (1937). 

) W. Glaser u. H. Robl, ZAMP, Vol. II, 444 (1948). 

3) P. Schiske, Nature 171, 443 (1953). Ausführlicher W. Glaser u. P. Schiske, 
Optik 11, 422 (1954). 

*) Die hier angewendete Methode kann in mehrfacher Weise auf den paraxialen 
Bereich spezialisiert werden. Zunächst kann man die hier gefundenen Felder für die Um- 
gebung einer geraden Feldachse entwickeln. Man kann aber auch in (1) von vorherein 
auf das paraxiale Gebiet spezialisieren. Schließlich kann man von der geraden Feldachse 
Beben und z. B. die um die Linien u = const und v = const entwickelten Felder be- 
Tachten. 
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wee 
steht sich von selbst, daß nur analytisch darstellbare Felder in Betracht 
gezogen werden. 


2. Kriterien für die Separierbarkeit der Hamilton-Jacobischen 
Differentialgleichung (HJDG) 


Der Grundgedanke der vorliegenden Arbeit besteht in folgendem. Wenn 
es gelingt, die HJDG für die Partikelbewegung in einem bestimmten Feld 
durch Summenansatz zu separieren, so ergeben sich die Bahnen in diesem 
Feld automatisch durch Quadraturen und sind damit in dem oben definierten 
Sinn „streng‘‘ berechenbar. Unsere Aufgabe kann also auch dahingehend 
formuliert werden, daß jene elektrostatischen Felder zu suchen sind, die zur 
Separierbarkeit der HJDG führen. Wir sind damit zu einem Sonderfall des 
Problems von Staeckel5) gelangt. Mit dem hier skizzierten Verfahren 
werden wir nicht unbedingt alle Felder des oben definierten Typus erfassen, 
doch wird sich zeigen, daß wir zu einer Reihe von Feldern gelangen, die unter 
Umständen praktische Bedeutung besitzen. Außerdem umfassen wir damit 
systematisch die bisher bekannten Fälle. 
Zunächst sollen in kurzer Form Kriterien für die Separierbarkeit der 
HJDG angegeben werden. Wir legen ein orthogonales — im allgemeinen 
krummliniges — Koordinatensystem (u, v, w) zugrunde, in dem die HJDG 
die Gestalt 

A, S2 + A, S} + A, S2—-2em-y =E (1) 
hat. Die Wirkungsfunktion S (u, v, w, ¢,¢, E) soll die Form Don 
S=U (u) + Ve) + W (w) (2) 
annehmen. Dafür ist notwendig und hinreichend)®), daß sich das Potential 
in der Gestalt 


2em-@p (u,v, w) = A, P (u) + 4,Q (v) + Ay R (v) (3) 
t darstellen läßt, während die A, (u, v, w) die Form 1 be 

> A, => (k = 1, 2, 3) on (4) 
J aufweisen. Dabei ist unter A die Determinante aus den neun Funktionen 


i, zu verstehen, in der die Funktionen der ersten Spalte nur von u, die der 
zweiten Spalte nur von v und die der dritten nur von w abhängen. Sonst sind 
diese Funktionen willkürlich. Unter M, sind die Minoren der ersten Zeile 
von A zu verstehen. 

Die durch (4) charakterisierten Koordinatensysteme wurden durch 
Eisenhart näher untersucht’). Danach gibt es (im euklidischen Raum) 
nur 11 verschiedene derartige Systeme. Sie werden z. T. in dieser und in 
einer folgenden Arbeit herangezogen. 

Die hier gefundenen Feldtypen erlauben übrigens auch die Separation 
der für die Einkörper-Bewegung im elektromagnetischen Feld zuständigen 
Schrödinger-Gleichung®). 


5) P. Staeckel, Habil. Schrift, Halle 1891 und Math. Ann. 42, 546 (1893). 

6) G. Hamel, Theoretische Mechanik, Springer 1949, S. 359. 

’)L.P. Eisenhart, Ann. Math. 35, 284 (1934), ferner Physic Rev. 45, 427 (1934) 
und 74, 87 (1948). 

*) H. P. Robertson, Math. Ann. 98, 749 (1927). 
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Fir die Partikeloptik sind vor allem die ebenen Felder (Elektronen- 
Prismen und -Zylinderlinsen) und die rotationssymmetrischen Felder 
(ot. symm. Elektronenlinsen) wichtig. In dieser Arbeit wird das ebene Feld 
in den Vordergrund gestellt. Die rotationssymmetrischen Felder werden in 


einem zweiten Teil behandelt. 


3. Bedingungen für Feld und Koordinaten 


Das ebene elektrostatische Feld ist von einer Koordinate, z. B. von 
w=z unabhängig. Zweckmäßigerweise wird man allgemeine Zylinder- 
koordinaten verwenden, die sich aus ebenen Orthogonalkoordinaten durch 
Hinzunahme der w-Achse ergeben. Die ebenen Koordinaten u, v werden durch 
eine analytische Funktion f aus kartesischen erzeugt: 


utiv=f(r+iy). 
Das Linienelement sei gegeben durch 


de = 5 (du? + dv?) + dw. 


Wir haben also 

A, =4,=Aluv) und A; = 1. (7a, b) 
Aus (7b) in Verbindung mit (4) ergibt sich A = M,; die Determinante A ist 
also von w unabhängig. Aus (7a) in Verbindung mit (4) folgt M, = M,. 
Da A nach (7) ebenfalls nicht von w abhängt, erkennen wir, daß 93, und 933 
konstant sein müssen. Die Beziehung M, = M, kann also in der Gestalt 


Pea — © Pas = © Piz — Pie (8) 
geschrieben werden. Da links eine Funktion von u allein und rechts eine von 
allein steht, müssen beide einer Konstanten & gleich sein. Mit Hilfe von (8) 
kann man Z. B. 9, und gy, als Funktionen von 95; bzw. gg ausdrücken und 
in A sowie in M, und M, einsetzen. Es zeigt sich, daß der Kehrwert von A als 
Summe einer Funktion p(w) und einer Funktion q (v) darstellbar sein muß: 


= p(w) +4 (0). 


Für das Potential haben wir nach (3) zu fordern: NE ee) re 
2em-g=A[P (u) +Q (v)]. 


Nun haben wir noch sicherzustellen, daß g die Potentialgleichung ae 
A (uu + Pov) = 0 (11) 


befriedigt. Durch Einsetzen von (10) in (11) gelangt man zu der Funktional- 
gleichung 
(Ayu + Ave) (P + Q) + 2(A, P’ + A,Q’) + +Q") =90 (12) 

für P (u) und Q (v). Da nach (10) P nur von u und Q nur von v abhängen darf, 
treten zu (12) die Nebenbedingungen 

Bao (13) 
Haben wir P und Q so bestimmt, daß sie (12) und (13) genügen, so kennen wir 
nach (10) ein Potential, für welches die HJDG separierbar ist. 
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an t. Berechnung der Bahnen 
ie Unter der Voraussetzung, daß das Koordinatensystem so gewählt ist. 
Be“ daß (9) zutrifft und daß das Potential die Form (10) hat, läßt sich (1) sepa- 
er = rieren und wir erhalten die Wirkungsfunktion 
u 
S(u, v, w) =| — a) + | P—xp+Adu+ | VQ—xq—Adr. (14 
Br = Es wurde dabei festgesetzt, daß S im Anfangspunkt (up, v, %,) verschwindet 
FR und das Potential wurde so normiert, daß es dort verschwindet, wo die Ge- 
sehwindigkeit Null wird. 
2B Von (14) gelangen wir mit Hilfe des Jacobischen Satzes zu den Bahn- 
 gleichungen (die Integranden von (14) sollen dabei den erforderlichen Ste- 
tigkeitsbedingungen genügen): N 
u = 
yx J J (19 
ze Die Parameter x und A lassen sich leicht durch Anfangspunkt und Anfangs- 


; richtung im (u, v, w)-System ausdrücken. Implizite Ableitung von (16) und 
al Ableitung von (15) führen zu SER 


dv P—-xp+4 VP-xp+i 
= Die Auflosung ergibt 
(Po + Yo) + (1 + | 
(Po + qo)” (Qo — Py %?) + ws (Vo Po — Po %) (19) 
(Po + %) (Po + Yo) (1 U9") wi? 
. . . . 
a Bei praktischen Rechnungen wird man zweckmäßigerweise den Übergang 
ae vom (u,v, w)-System zum kartesischen möglichst spät vollziehen. Dabei 
= kann man sich vorteilhaft der Konformität der durch f bewirkten Abbildung 
bedienen. 
2 Zur Verdeutlichung der hier entwickelten Methode sollen nun einige ele- 
i mentare Beispiele gegeben werden. Jedes Koordinatensystem wird dabei 
‘c Potentialfelder von mehr oder weniger praktischem Interesse liefern. Die 
hier gefundenen Felder sollen in einer späteren Arbeit im einzelnen untersucht 
werden, soweit sie von Interesse sind. BEN: 


5. Kartesische Koordinaten a. 

Dieser triviale Fall ergibt sich, wenn man in (5) die identische Transfor- 
mation benutzt. Es gilt dann p = 0, q = 1 und die Funktionalgleichung (12) 
lautet 
P” (x) + Q” (y) = 0. (20) 
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Sie kann unmittelbar separiert werden und führt auf das Potential £ 
+6, 2+ Coy + Cg. (21) 


Man kann zwei Fälle unterscheiden. Ist k? = 0, so erhält man nach einer Um- 
formung (wir behalten die Bezeichnung der Variablen auch nach der Drehung 
bei) das homogene elektrostatische Feld: 


2emg=ar+b. (22) 
Ist dagegen k? + 0, so ergibt sich nach Umformung das „Hyperbelfeld‘“ 
2em-g = a? (a? — y*) + b (23) 


das ein Modell fiir eine elektrostatische Zylinderlinse darstellt. 


« 


3 


Wählt man in (5) den Logarithmus, so ergeben sich Kreiszylinder-Koor- 
mit 


= Gos Vv 

Die (12) lautet 
123 [23 | [4 
P r Q } (P Q ) == chi (25) 

Sie kann separiert und durch ne 
Q cos 2 (v — cy) (26) 

2 

gelést werden. 
Wenn man mit 

t=» (27) 


zur iiblichen Schreibung der Polarkoordinaten iibergeht, ergibt sich das 
Potential 


2em:g(r,y) =a+binr+eos2(y —d). (28) 


Dieses Potential ist aus dem gewöhnlichen logarithmischen Potential und 
einem winkelabhängigen Anteil zusammengesetzt. 


bar 


7. Koordinaten des parabolischen Zylinders Wake 


(29) 


In diesem Fall ist für (5) 
zu setzen, woraus weiter 


y = = uv P q eal ) 
ZS 


folgen. Die Funktionalgleichung (12) 
(P” + Q”) +0) —4(P’u + Q’ + 4(P4+ Q) = 0 (31) 


ist hier nicht unmittelbar separierbar, wie in den früheren Fällen. Wir müssen 
von den Nebenbedingungen (13) Gebrauch machen. Wenn man (31) zuerst 
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nach w und dann nach »v ableitet, ergibt sich die Beziehung 
2 v + 2Q” u 0, 
die zu 
1 
P ut cgut Q=- +560 + 6,0 + G (33) 


fiihrt. Durch den DifferentiationsprozeB sind wir allerdings zu allgemeineren 
Lésungen gelangt und miissen daher (33) wieder in (31) einsetzen. Es zeigt 
sich, daß k?,c, und c, beliebig gewählt werden können, daß jedoch die Be- 
ziehungen c, = —c, und c, =c, bestehen. Man gelangt so zum Potential 


cu + dv 


2em-p(u, v) = a(u? — v2) +O 


(34) 


Hier ist nur der letzte Ausdruck von Interesse, weil der erste, wie die 
Einführung kartesischer Koordinaten zeigt, ein homogenes Feld darstellt. 


8. Koordinaten des elliptischen Zylinders 
In diesem Falle hat man in (5) den Arcussinus zu nehmen. Es ist dann 
x=sinu-chv 

und p=-—sijn?u; q=ch*v. (35) 
y= teosu-shv 
Die Funktionalgleichung (12) 
+ (c h? uw — sin? v) (P” + Q”) + P’ sin u cos u ' 
—Q' shvcehv + (ch? v + sin? u — 1) (P + Q) = 


kann wiederum nicht unmittelbar separiert werden. Differenzieren nach u 
und anschließend nach » führen zu der Beziehung 


(P'" +4P) shv-chv = (Q’" —4Q’ sin u - cos u, 
welche die Lésungen 


P =c,3in 24+ c, cos 24 —— 
2 3 


(36) 


u-sin2u-+ C4 


16 


Q=@sh2v+% ch2v+ -sh2v+e, 


besitzt. Das Einsetzen in (36) zeigt, daß k?, c, und c, beliebig gewählt werden 

| men, während die Beziehungen c, = c, und cy = — c, zu berücksichtigen 
I, Auf diese Weise gelangen wir zu ue 

1 


2em-:@ (u,v) = a+ [b-sin2u+c-sh2v 


ch? v — sin? u 
+d(v-sh2v—u-sin 2u)]. 

Das erste Glied dieses Ausdrucks stellt den interessanten Fall eines Ablenk- 


feldes dar, das eine exakte Bahnbestimmung zuläßt, das zweite dagegen eine 


LEE, 
< 
Bei der Redaktion eingegangen am 2. Juli 1955. 
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Ro tationssymmetrische elektrostatische Felder, 
die eine strenge Berechnung der Elektronenbahnen zulassen 


Inhaltsübersicht 


Es werden rotationssymmetrische elektrostatische Felder gesucht, die 
eine Bestimmung der Trägerbahnen durch Quadraturen zulassen. Dazu ist es 
nötig, die von Staeckel angegebenen Feldtypen weiteren Bedingungen zu 
unterwerfen. Die sich ergebenden Felder werden explizit Rs 


1. Einleitung 


In einer früheren Arbeit!) wurden ebene elektrostatische Felder bestimmt, 
die eine strenge Berechnung der in ihnen auftretenden Partikelbahnen zu- 
lassen. Dabei wurde eine Bahn, die durch einfache Quadraturen gewonnen 
werden kann, als „‚streng‘“ berechenbar aufgefaßt. Dieses Problem wurde als 
Sonderfall des Problems von Staeckel?) behandelt, indem es wir folgt for- 
muliert wurde: Es sind jene elektrostatischen Felder zu suchen, die zur Se- 
parierbarkeit der Hamilton-Jacobischen Differentialgleichung (HJDG) 
führen. Die in !) benutzte Methode soll nun auf das praktisch wichtige rota- 
tionssym metrische Feld (Elektronenlinse) übertragen werden. 

Wir erinnern kurz an die Kriterien der Separierbarkeit?). Es wird ein 
orthogonales — im allgemeinen krummliniges — Koordinatensystem (u, v, w) 
zugrunde gelegt, indem die HJDG die Gestalt 


A, S2 + A, S2 + A, S2 = 2em-p+E (1) 


hat. Wenn sich (1) durch einen Summenansatz A en 


(2) 


separieren lassen soll, ist notwendig und hinreichend, daB sich das Potential 
in der Gestalt 


2em-¢ (u,v, w) = A, P (u) + Ag Q (v) + Ag (w) (3) 
darstellen läßt, während die A, die Form er lem 
M. 

A, (k = 1, 2, 3) (4) 


1) H. Ann. Physik 17, 269 (1956). 
2) P. Staeckel, Habil.-Schrift, Halle 1891 u. Math. Ann. 42, 546 (1893), vgl. a. 
G. Hamel, Theoretische Mechanik, Springer 1949, S. 359. 
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2 aufweisen. Dabei ist unter A die Determinante aus den neun Funktionen Pri 
zu verstehen, in der die Funktionen der ersten Spalte nur von u, die der zweiten 
nur von v und die der dritten nur von w abhängen. Sonst sind diese Funk- 
tionen willkürlich. Unter M, sind die Minoren der ersten Zeile von A zu ver- 
stehen. Die erwähnten Bedingungen haben wir nun auf rotationssymmetrische 
Felder zu spezialisieren. ef. 


2. Bedingungen für Koordinatensystem und Feld 


ZweckmaBigerweise werden allgemeine rotationssymmetrische Koordinaten 
verwendet. Sie können aus ebenen Orthogonalkoordinaten, definiert durch 


utiv=f(z+io) (5) 


_erzeugt werden, indem man diese um die z-Achse rotieren läßt. Das Linien- 
element ist dann gegeben Br 


ds? =4 (du? + dv?) + A, dw, (6) 


Aus der erate (4) kann man mit (7) folgern, daB mit 


-A = 0° — Poe Pia) at 
A, = Ag... Poe Pas — 23 = Fis — Pre (9) 


gelten muß. Nach (8) darf A die Variable w nicht enthalten. Da nun A eben- 
falls von w unabhängig sein muß, folgt, daß die Funktionen gg, und g,, kon- 
stant sein müssen. Es ist daher möglich, in (9) etwa 9, und gy, durch 9 
bzw. gg auszudrücken. Wenn man die erhaltenen Beziehungen wieder in (4) 
einsetzt, erkennt man, daß A, und A folgende Gestalt besitzen müssen: 


1 
A, A (10) 
p(u) und g(v) sind dabei willkürliche Funktionen der angeschriebenen Ar- 
gumente. Aus (3) folgt dann sofort die Bedingung für das Feld re 
(ke 
2em g=A[P(u +Q (v)] +5 Rw). 


Da das Feld rotationssymmetrisch sein soll, also vom Azimut w nicht ab- 
hängen darf, ist in (11) die willkürliche Funktion R (w) Null zu setzen. Das 
Potential soll ferner so normiert werden, daß es dort verschwindet, wo die 
Partikelgeschwindigkeit Null wird. Das führt auf E = 0. 

Haben wir nach Wahl eines der Bedingung (10) entsprechenden Koordi- 
natensystems*) passende Funktionen P (u) und Q (v) gefunden, so kann die 
HJDG (1) separiert werden. Das Potential (11) muß aber noch die Potential- 
gleichung befriedigen. Für eine aus leitenden Elektroden im Vakuum be- 


8) Wir folgen in der Bezeichnung E. Madelung, Die math. Hilfsm. des Phys., 
Springer 1953. 


4) Vgl. die Anmerkung *) in der Arbeit }). 
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stehende Linsenanordnung lautet diese Gleichung fiir unser Problem 


du (0 + dv \@ (12) 
Durch Einsetzen von (11) in (12) gelangen wir zu der fundamentalen Funk- 


tionalgleichung 
(P + Q) [o, Ay + Ap +0 (Au + + P’ (Q, A + 20 4,) 


” (13) 
+ (0, A + 20 A,) +0 A (P” + Q”) = 0 
i 
fir P (wu) und Q (v). Dazu treten die Nebenbedingungen 


3. Berechnung der Bahnen 


Unter der Voraussetzung, daß wir das Koordinatensystem nach (10) ge- 
wählt und ein Feld gefunden haben, das den oben formulierten Bedingungen 
entspricht, läßt sich (1) separieren und wir kommen zu der Wirkungsfunktion 


(u, v, w) = (w — wy) + + Adu+ dv. 
(15) 


Sie stimmt formal mit der für das ebene Problem in !) gefundenen überein. 
Auch für die mit Hilfe des Jacobischen Satzes gefundenen Bahngleichun- 
gen und für die Festlegung der Integrationskonstanten gilt das gleiche. Wir 
verweisen diesbezüglich auf !). Zur Illustration des hier dargelegten Ver- 
fahrens sollen die folgenden elementaren Beispiele dienen. Die sich dabei er- 
gebenden Felder werden in einer späteren Arbeit im einzelnen untersucht. 


4. Zylinder-Koordinaten 
Dieser Fall ergibt sich, wenn man in (5—7) 
setzt. Aus (10) folgt 
2=%; 
Die Funktionalgleichung (13) lautet 
Q’ + 0 (P” =0 
Sie kann unmittelbar separiert und durch 
P=--Zw + Cy & + Cg; Q +c, (19) 
gelöst werden. Mit (11) ergibt sich das Potential 
2em-y = a* (o? — 22) +b-Ine+cz+d. (20) 


Dieses Potential ist additiv aus einem ,,Hyperbelfeld‘‘, aus dem gewöhnlichen 
Zylinderpotential und schließlich aus einem homogenen Anteil zusammen- 
gesetzt. 


i 
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5. Kugelkoordinaten 


Wählt man in (5) den Logarithmus, so erhält man Kugelkoordinaten mit 


Nach (10) gilt gene 

Die Funktionalgleichung (13) ae 
2(P+Q)-3P'+Q etgv+ P” +Q” =0 (23) 


kann unmittelbar separiert und durch 


‘ f v 1 
P = c, + ce" + = kh’; Q = cos v + ¢5(1 + cos v-In 


gelöst werden. Mit (11) erhalten wir das Potential 
2em-y =a+be"+c-e cos v + de?" (1 + cos v-In (25) 


das unter anderem das Newtonsche Potential e-“ enthält. 


6. Rotationsparabolische — 


Mit 
u+iv —V2 + io) 


ergeben sich die rotationsparabolischen Koordinaten 


_ 
ur 0 = Uv 
1 
(28) 


Die Funktionalgleichung (13) lautet hier 


um)" —(, 9 
Sie kann nicht unmittelbar separiert werden. Um das zu erreichen, benutzen 


wir die Nebenbedingungen (14) und differenzieren (29) zuerst nach u und dann 
nach v (die Reihenfolge ist übrigens gleichgültig). Es ergibt sich die Beziehung 


2+ 3Pu+ 3Q'v=0 (30) 


die separiert und durch 
P=" inu+ Cg t+ Q =—fino +, (31) 


gelöst werden kann. Wegen des Differentiationsprozesses ist natürlich die 
Lösung von (30) allgemeiner als die der ursprünglichen Funktionalgleichung 
(29). Es muß daher (31) erneut in (29) eingesetzt werden, um die nötigen 
_ Spezialisierungen zu erkennen. Man findet so, daß %®, c, und c, willkürlich 
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bleiben, während andererseits die Beziehungen c, = c, und c, = — c, berück- 
siehtigt werden müssen. Man gelangt so zum Potential 
mit 
= — v2 -In— if 
2em-_ (u,v) =a+b(u In +). (32) 
Unter Benutzung von (27) erkennt man, daß das Glied mit 6 ein homogenes - 
os Feld liefert. in 
7. Koordinaten des abgeplatteten Rotationsellipsoides i 
9% 
25) Wird in (5) der Arcussinus gewählt so ergeben sich die Koordinaten des . 
abgeplatteten Rotationsellipsoides: 5 
24) 0 sinu-chv x? (c h? v.r— sin? u) ( ) 
und weiter 2 
25) Pints.’ (34) 
Die Funktionalgleichung 
2(ch?v + sin? u) (P + Q) + (ch? v- etgu + 3sin u- cos u) P’ a me 
—(3chv-shv+thv-sin® u) Q + (ch? v— sin? u) (P’ + Q”) = (39) 
Y kann ebenfalls nicht unmittelbar separiert werden. Ableiten nach w und 
_ danach nach v führt zu der separierbaren Beziehung 
shi chv(2 P’ — etgu +P ) 
3 (86) 
+ sin u - cos u(2Q Q” thy — = 
8) die durch 
1-+ cos u 16 
(37) 
9) =agshv + byshv- aretg (sh v) +7 v 
- gelöst wird. Auch hier ist erneutes Einsetzen in (35) erforderlich und führt zu 
aie b, = — b,/2, während a,, a, und k* beliebig gewählt werden können. Wir 
ng kommen so zum Ausdruck 
1 
0) ~ zem (u, v) d a2 (c h2 sin? u) 
1 — cosu 
acosu+bshv+ c[ cos u In —2shv-aretg (she 
1 + cos u 
1) 
Setzt man in (38) d = 0 und c = 0. so ergibt sich das von P. Schiske°) 
ie gefundene Feld. 
1g 
an 5) P. Schiske, Nature 171, 443 (1953). Ausführlicher in W. Glaser u. P. Schiske, 
+h Optik 11, 422 (1954). 
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8. Koordinaten des verlängerten Rotationsellipsoids — 
Ganz analog verläuft die Rechnung für die durch 


mt A=- = 
o=xcosu-shv x? (ch? v — sin? u) 
0 


| Koordinaten des verlängerten | 


sich die Feldform: 
1 


c h Ur 1 | 

Damit sind die in unserem Problem anwendbaren Koordinatensysteme®) voll- 
ständig erfaßt. 


VI, Joanelligasse 7. 
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Scheinprobleme der Wissenschaft 
Vierte unveränderte Auflage. 31 Seiten. 1955. DIN A 5. DM 1.50 


Wissenschaftliche Selbstbiographie 
Dritte Auflage.34 Seiten mit einem Bildnis. 1955. DIN A5. DM 1.35 


[Bildet Lebensdarstellungen Deutscher Naturforscher, Nr. 5. Herausgegeben von 
der Deutschen Akademie der Naturforscher (Leopoldina) | 


„Zeitschrift für Naturforschung‘‘: „Plancks Reden gehören wesentlich zu dem 
Gesicht, das die heutige Physik dem geistig interessierten Deutschen zeigt, — sie 
gelten weiterhin als Zeugnis für den geistigen Rang der Naturwissenschaft.‘ 

W. Kossel 


„Kusmos‘“: „In den acht Vorträgen, die der große deutsche Gelehrte, der Schöpfer 
der Quantentheorie und Nobelpreisträger in den Jahren zwischen 1929 und 1946 
meist vor einer breiten Hörerschaft gehalten hat, beschäftigt er sich von hoher 
Warte aus mit Fragen nach den letzten Problemen der Wissenschaft und des 
menschlichen Seins aus dem Gesichtswinkel seines physikalischen Weltbildes, Die 
gedrängte Kürze der Darstellung, bedingt durch die Vortragsform, und der glas- 
klare Stil erleichtern dem naturwissenschaftlichen Laien das Eindringen in die 
schwierigen Probleme.“ D. Herzog 


JOHANN AMBROSIUS BARTH / VERLAG / LEIPZIG 
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Annalen der Physik. 6. Folge, Band 17. Heft 4-5. 1956 


The Physical Society 
London 


Die Mitgliedschaft der Gesellschaft kann jeder erwerben, der sich fiir 

die wissenschaftliche Physik interessiert. Die Mitglieder erhalten alle — 
Veröffentlichungen der Gesellschaft zu wesentlich ermäßigten Sonder- | 
preisen. Ferner wird jedem Mitglied eine Eintrittskarte zur jährlich 
im Frühjahr stattfindenden Londoner Ausstellung der Gesellschaft 
kostenlos überlassen. 


Zu den Veröffentlichungen der Gesellschaft gehören: 


Proceedings of the Physical Society 
Monatlich 2 Hefte und zwar: 


Section A (Physics of Elementary Processes) a 
Section B (Physics of Collective Processes) 


Reports on Progress in Physics 


Jährlich ein Band mit etwa 10 Ubersichtsberichten über den 
gegenwärtigen Stand der Forschungen auf verschiedenen 
physikalischen Gebieten. 


Handbook of Scientific Instruments ER 
and Apparatus 


Erscheint jährlich in Verbindung mit der Ausstellung wissen- 
schaftlicher Instrumente und Apparate. 


= 


Sonderberichte 


Im Jahre 1955 wurden folgende Tagungsberichte herausgegeben: 


Defects in Crystalline Solids (Bristol, Juli 1954) 
The Physies$of the Ionosphere (Cambridge, September 1954) 


Anträge auf Erwerb der Mitgliedschaft der Gesellschaft sind zu richten an: 


THE PHYSICAL SOCIETY, 1 LOWTHER GARDENS, 
PRINCE CONSORT ROAD, LONDON S.W. 7 
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